Contact Dynamics: Legendrian and Lagrangian Submanifolds

Abstract

We are proposing Tulczyjew’s triple for contact dynamics. The most important ingredients of the triple, namely symplectic diffeomorphisms, special symplectic manifolds, and Morse families, are generalized to the contact framework. These geometries permit us to determine so-called generating family (obtained by merging a special contact manifold and a Morse family) for a Legendrian submanifold. Contact Hamiltonian and Lagrangian Dynamics are recast as Legendrian submanifolds of the tangent contact manifold. In this picture, the Legendre transformation is determined to be a passage between two different generators of the same Legendrian submanifold. A variant of contact Tulczyjew’s triple is constructed for evolution contact dynamics.

1. Introduction

Lagrangian dynamics are generated by a Lagrangian function defined on the tangent bundle $T\mathcal{Q}$ of the configuration space of a physical system, whereas Hamiltonian dynamics are governed by a Hamiltonian function on the cotangent bundle $T^{\ast }\mathcal{Q}$, which is canonically symplectic [1,2,3,4]. If a Lagrangian function is regular, that is, if it satisfies the Hessian condition, then the fiber derivative becomes a fibered local diffeomorphism from the tangent bundle to the cotangent bundle. In this case, the fiber derivative turns out to be the Legendre transformation linking the Lagrangian and the Hamiltonian realizations of the physical system.

If a Lagrangian function happens to be degenerate, then the fiber derivative fails to be a local diffeomorphism since its image space turns out only to be, in the best of cases, a proper submanifold of the cotangent bundle $T^{\ast }\mathcal{Q}$. That is, one only arrives at a presymplectic picture determined by some constraint functions. To deal with these constraints, Dirac proposed an algorithm, nowadays called the Dirac–Bergmann algorithm [5,6]. This algorithm proposes a method to arrive at a submanifold (possibly smaller than the image space of the Legendre transformation) of the cotangent bundle where the Hamilton’s equations becomes well-defined. In the final stage of the algorithm, one obtains the so-called Dirac bracket. There also exists a more geometric version of this approach called the Gotay–Nester–Hinds algorithm [7]. Inspired the tools introduced in [7], the Skinner–Rusk unified theory [8] is establishing a unification of Lagrangian and Hamiltonian formalisms on the Whitney sum of tangent and cotangent bundles. In this paper, we shall focus on the Tulczyjew approach for the Legendre transformations of singular Lagrangians.

The Classical Tulczyjew’s Triple. Tulczyjew’s triple is a commutative diagram linking three symplectic bundles, namely, $T{T}^{\ast }\mathcal{Q}$, $T^{\ast }{T}^{\ast }\mathcal{Q}$ and $T^{\ast }T\mathcal{Q}$ via symplectic diffeomorphisms [9]. This geometrization enables one to recast Lagrangian and Hamiltonian dynamical equations as Lagrangian submanifolds of the Tulczyjew symplectic space $T{T}^{\ast }\mathcal{Q}$ [10]. Referring to this geometry, the Legendre transformation is defined as a passage between two different generators of the same Lagrangian submanifold [11,12]. This definition is free from the non-degeneracy requirement that is the Hessian condition. Evidently, this approach is in harmony with the creed by Weinstein “everything is a Lagrangian submanifold” [13]. We reserve Section 2 for a brief summary of the Tulczyjew’s triple and the Legendre transformation in this picture.

Tulczyjew’s triple is modified for many physical systems and it is carried to several geometric frameworks. For higher order Lagrangian dynamics, the triple is upgraded in [14,15]. For physical theories where the configuration space is a Lie group, the triple is examined in a series of works [16,17,18,19]. The triple is examined for principal fiber bundles in [20,21]. It is written for the vector bundle of n-vectors in [22]. For the case of the field theories, we refer to an extensive, but incomplete list [23,24,25,26,27,28,29,30,31,32], and for the higher order field theories see, for example, [33].

Contact Hamiltonian Dynamics. A symplectic manifold must be even dimensional. An odd dimensional generalization of symplectic geometry is contact geometry [34]. Hamiltonian dynamics are available on this generalization as well [35,36,37,38]. In the present work, as a generic model of contact manifold, we consider the extended cotangent bundle

$${\mathcal{T}}^{\ast }\mathcal{Q}=T^{\ast }\mathcal{Q}\times \mathbb{R}$$

(1)

with the contact one-form ${\eta }_{\mathcal{Q}}=dz-{\theta }_{\mathcal{Q}}$. Here, z is the real variable, and ${\theta }_{\mathcal{Q}}$ is the pull-back of the canonical one-form on the cotangent bundle $T^{\ast }\mathcal{Q}$. In a coordinate free formulation, contact Hamiltonian dynamics generated by a Hamiltonian function H are defined as

$${\iota }_{X_H^c}{\eta }_{\mathcal{Q}}=-H,{\iota }_{X_H^c}d{\eta }_{\mathcal{Q}}=dH-\mathcal{R}\left(H\right){\eta }_{\mathcal{Q}},$$

(2)

where $\mathcal{R}$ is the Reeb field associated with the contact form ${\eta }_{\mathcal{Q}}$, and ${\iota }_{X_H^c}$ is the interior product. In Darboux’ coordinates $(q^i,p_i,z)$, the contact Hamilton’s equations are computed to be

$${\dot{q}}^i=\frac{\partial H}{\partial p_i},{\dot{p}}_i=-\frac{\partial H}{\partial q^i}-p_i\frac{\partial H}{\partial z},\dot{z}=p_i\frac{\partial H}{\partial p_i}-H.$$

(3)

Contact Hamiltonian dynamics have different features from classical Hamiltonian dynamics. One interesting characteristic of contact Hamiltonian dynamics is the loss of conservation of Hamiltonian function along the motion [39,40]. Even, the canonical contact volume form is not preserved under the action of the contact Hamiltonian dynamics (see [41]). The dissipative nature of contact Hamiltonian dynamics makes them proper for dissipative dynamical systems. We present here an incomplete list of some recent works along this direction: [37,39,40,42]. Additionally, we can argue that contact framework looks proper for thermodynamics, see, for example, the following incomplete list: [35,43,44,45]. Accordingly, we also wish to cite [46,47,48,49] for the use of contact geometry for qualitative analysis of reversible–irreversible dynamics under GENERIC (an acronym for General Equation for Non-equilibrium Reversible-Irreversible Coupling) formalism. At the end of this paper, we shall provide some applications of our results on thermodynamics. We cite [50,51] for some recent studies on variational aspects of contact Hamiltonian dynamics.

There is evolution contact Hamiltonian formalism on the extended cotangent bundle. In this case, for a Hamiltonian function H, evolution contact Hamilton’s equation is defined to be

$${\iota }_{{\epsilon }_H}{\eta }_{\mathcal{Q}}=0,{\mathcal{L}}_{{\epsilon }_H}{\eta }_{\mathcal{Q}}=dH+\mathcal{R}\left(H\right){\eta }_{\mathcal{Q}}.$$

(4)

Here, ${\mathcal{L}}_{{\epsilon }_H}$ is the Lie derivative. According to the Cartan’s formula, it is computed to be ${\mathcal{L}}_{{\epsilon }_H}=d{\iota }_{{\epsilon }_H}+{\iota }_{{\epsilon }_H}d$. In Darboux’s coordinates, the local picture of the evolution contact Hamilton’s Equation (4) is computed to be

$${\dot{q}}^i=\frac{\partial H}{\partial p_i},{\dot{p}}_i=-\frac{\partial H}{\partial q^i}-p_i\frac{\partial H}{\partial z},\dot{z}=p_i\frac{\partial H}{\partial p_i}.$$

(5)

Note that if 0 is a regular value of H, then the evolutionary vector field ${\epsilon }_H$ is tangent to the hypersurface $H^{-1}\left(0\right)$ and, thus, to every Legendrian submanifold contained in $H^{-1}\left(0\right)$. Therefore, the integral curves of ${\epsilon }_H$ may be interpreted as thermodynamical processes for a system with thermodynamical phase space in the extended cotangent bundle. This fact was used in [52] to discuss the relation between evolution contact dynamics and simple thermodynamical systems with friction.

Contact Lagrangian Dynamics. Corresponding Lagrangian formalism for contact Hamiltonian dynamics is also available in the literature under the name of Herglotz (or generalized Euler–Lagrange) formalism [53]. See also [38]. In this paper, as the geometric framework, we consider the extended tangent bundle

$$\mathcal{T}\mathcal{Q}=T\mathcal{Q}\times \mathbb{R}.$$

(6)

A Lagrangian function $L=L(q,\dot{q},z)$ on $\mathcal{T}\mathcal{Q}$ determines the Herglotz (generalized Euler–Lagrange) equations as

$$\frac{\partial L}{\partial q^i}-\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}^i}\right) + \frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}=0.$$

(7)

Here, z is the real variable appearing in $T\mathcal{Q}\times \mathbb{R}$. Evidently, if the Lagrangian function L is independent of z, then Equation (7) reduces to the classical Euler–Lagrange equations. If the Lagrangian function L is non-degenerate, that is, if the rank of the Hessian matrix

$$\left[\frac{{\partial }^{2}L}{\partial {\dot{q}}^i\partial {\dot{q}}^i}\right]$$

(8)

is maximum, then the fiber derivative

$$\mathbb{F}{L}^c:\mathcal{T}\mathcal{Q}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^i,z)\mapsto (q^i,\frac{\partial L}{\partial {\dot{q}}^i},z)$$

(9)

turns out to be a local diffeomorphism. In this case, a direct calculation shows that the fiber derivative (9) maps the Herglotz equations in (7) to the contact Hamilton’s Equation (3) once the Hamiltonian function is taken to be

$$H(q^i,p_i,z)={\dot{q}}^i{p}_i-L(q,\dot{q},z).$$

(10)

If a Lagrangian function fails to be non-degenerate then, as in the symplectic case, the transformation (9) fails to be a local isomorphism. Then the image space of $\mathbb{F}{L}^c$ can only be, in the best of the cases, a proper submanifold of the extended cotangent bundle. One way to deal with the constraints defining the image space is to employ a version of Dirac algorithm. This is recently studied in [54]. Another way is to generalize the unified formalism for contact dynamics as presented in a recent study [55]. In this paper, our interest is to study the Legendre transformation for contact dynamics following the understanding of Tulczyjew.

Goal of the Present Work. The aim of this work is to define Legendre transformation between the Herglotz Equation (7) and the contact Hamilton’s Equation (3) by properly constructing a Tulczyjew’s triple for the case of contact manifolds. We shall call this as contact Tulczyjew’s triple. Such an attempt involves modifications of the ingredients of the classical Tulczyjew’s triple to the contact geometry. Here, in the contact case, the role of Tulczyjew’s symplectic space $T{T}^{\ast }\mathcal{Q}$ will be played by the extended tangent bundle $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ of the extended cotangent bundle. Accordingly, Tulczyjew’s triple will consist of the iterated extended bundles ${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}$, $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$, ${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}$ as well as contact transformation between them. Maybe, the most vital object in this picture is the introduction of the notion of special contact manifolds. This novel framework permits one to recast both the contact Lagrangian and the contact Hamiltonian dynamics as a Legendrian submanifold of $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$. This enables us to realize that Lagrangian and Hamiltonian functions as generating objects of the same Legendrian submanifold. Thus, by merging special contact geometry with Morse family theory, the Legendre transformation for contact dynamics is defined to be a passage between two different generators of the same Legendrian submanifold. In addition, there is an (evolution) Hamiltonian flow on contact geometry preserving the energy, but not the kernel of the contact form. There is also Lagrangian counterpart of this theory called evolution Herglotz equations. By properly modifying contact Tulczyjew’s triple, the Legendre transformation for the evolution Herglotz equations and the evolution contact Hamilton’s equations are obtained. In this theory, contact manifolds and Legendrian submanifolds are replaced by symplectic manifolds and Legendrian submanifolds, respectively. We call this geometry the evolution contact Tulczyjew’s triple.

The content of this work is as follows. The main body of the paper consists of three sections. In Section 2, for the sake of the completeness of the manuscript and in order to fix the notation, a brief summary of classical Tulczyjew’s triple is given. Section 3 is reserved for the basics on contact dynamics in both Hamiltonian and Lagrangian formulations. Section 4 is the one containing the novel results of the paper where the Tulczyjew’s triple is constructed for the contact and evolution contact dynamics.

2. The Classical Tulczyjew’s Triple

2.1. (Special) Symplectic Manifolds

A manifold $\mathcal{P}$ is said to be symplectic if it is equipped with a non-degenerate closed two-form $\omega$ [1,3,56]. In this case, $\omega$ is called a symplectic two-form. A diffeomorphism between two symplectic manifold is called a symplectic diffeomorphism if it respects the symplectic two-forms.

Submanifolds. Let $(\mathcal{P},\omega )$ be a symplectic manifold, and $\mathcal{S}$ be a submanifold of $\mathcal{P}$. We define the symplectic orthogonal complement of $T\mathcal{S}$ as the vector subbundle of $T\mathcal{P}$

$$T{\mathcal{S}}^{\perp }=\{X\in T\mathcal{P}:\omega (X,Y)=0,\forall Y\in T\mathcal{S}\}.$$

(11)

The rank of the tangent bundle $T\mathcal{P}$ is the sum of the ranks of the tangent bundle $T\mathcal{S}$ and its symplectic orthogonal complement $T{\mathcal{S}}^{\perp }$. We list some of the important cases.

• $\mathcal{S}$ is called an isotropic submanifold if $T\mathcal{S}\subset T{\mathcal{S}}^{\perp }$. In this case, the dimension of $\mathcal{S}$ is less or equal to the half of the dimension of $\mathcal{P}$.
• $\mathcal{S}$ is called a coisotropic submanifold if $T{\mathcal{S}}^{\perp }\subset T\mathcal{S}$. In this case, the dimension of $\mathcal{S}$ is greater or equal to the half of the dimension of $\mathcal{P}$.
• $\mathcal{S}$ is called a Lagrangian submanifold if $T\mathcal{S}=T{\mathcal{S}}^{\perp }$. In this case, the dimension of $\mathcal{S}$ is equal to the half of the dimension of $\mathcal{P}$.

Under a symplectic diffeomorphism, the image of a Lagrangian (isotropic, coisotropic) submanifold is a Lagrangian (resp. isotropic, coisotropic) submanifold.

The Cotangent Bundle. The generic examples of symplectic manifolds are cotangent bundles. To see this, consider a manifold $\mathcal{Q}$, and its cotangent bundle $T^{\ast }\mathcal{Q}$. The canonical (Liouville) one-form ${\theta }_{\mathcal{Q}}$ on $T^{\ast }\mathcal{Q}$ is defined, on a vector X over $T^{\ast }\mathcal{Q}$ as

$${\theta }_{\mathcal{Q}}\left(X\right)=\langle {\tau }_{T^{\ast }\mathcal{Q}}\left(X\right),T{\pi }_{\mathcal{Q}}\left(X\right)\rangle .$$

(12)

Here, ${\tau }_{T^{\ast }\mathcal{Q}}$ is the projection from the tangent bundle $T{T}^{\ast }\mathcal{Q}$ to its base manifold $T^{\ast }\mathcal{Q}$, whereas $T{\pi }_{\mathcal{Q}}$ is the tangent lift of the cotangent projection ${\pi }_{\mathcal{Q}}$. To be more precise, we present the following commutative diagram,

(13)

Minus of the exterior derivative of the canonical one-form ${\theta }_{\mathcal{Q}}$, that is, ${\omega }_{\mathcal{Q}}:=-d{\theta }_{\mathcal{Q}}$, is the canonical symplectic two-form on the cotangent bundle $T^{\ast }\mathcal{Q}$.

Hamiltonian Vector Fields. Let X be a vector field on the symplectic manifold $(T^{\ast }\mathcal{Q},{\omega }_{\mathcal{Q}})$. It is called a local Hamiltonian vector field if it preserves the symplectic two-form

$${\mathfrak{L}}_X{\omega }_{\mathcal{Q}}=0.$$

(14)

Poincaré Lemma assures us the existence of a local function H satisfying the so-called Hamilton’s equation

$${\iota }_{X_H}{\omega }_{\mathcal{Q}}=dH.$$

(15)

Notice that, in this realization, we denote the Hamiltonian vector field by $X_H$. If there exists a global Hamiltonian function H then $X_H$ is called a (global) Hamiltonian vector field.

If the dimension of $\mathcal{Q}$ is n then the cotangent bundle $T^{\ast }\mathcal{Q}$ turns out to be a $2n$ dimensional manifold. In this case, the n-th power ${\omega }_{\mathcal{Q}}^n$ of the symplectic two-form is a non-vanishing top-form on $T^{\ast }\mathcal{Q}$. So that, it determines a volume form, the so-called symplectic volume. The identity (14) gives that a Hamiltonian vector field preserves the symplectic volume. Further, the skew-symmetry of the symplectic two-form manifests the conservation of the Hamiltonian function H. If H is the total energy, then this is called the conservation of energy.

A set of natural bundle coordinates $(q^i,p_i)$ on $T^{\ast }\mathcal{Q}$ are Darboux coordinates, meaning that the canonical one-form ${\theta }_{\mathcal{Q}}$ and the symplectic two-form ${\omega }_{\mathcal{Q}}$ are written as

$${\theta }_{\mathcal{Q}}=p_{i}d{q}^i,{\omega }_{\mathcal{Q}}=d{q}^i\wedge d{p}_i,$$

(16)

respectively. In this realization, the Hamilton’s Equation (15) is computed to be

$$\frac{d{q}^i}{dt}=\frac{\partial H}{\partial p_i},\frac{d{p}_i}{dt}=-\frac{\partial H}{\partial q^i}.$$

(17)

Special Symplectic Structures. Let $\mathcal{P}$ be a symplectic manifold carrying an exact symplectic two-form $\omega =-d\theta$, with $\theta$ is being a potential one-form. Assume also that, $\mathcal{P}$ is the total space of a fiber bundle $(\mathcal{P},\pi ,\mathcal{Q})$. Here, $\mathcal{Q}$ is the base space, and $\pi$ is the projection. A special symplectic structure is a quintuple

$$(\mathcal{P},\pi ,\mathcal{Q},\theta ,\varphi )$$

(18)

where $\varphi$ is a fiber preserving symplectic diffeomorphism from $\mathcal{P}$ to the cotangent bundle $T^{\ast }\mathcal{Q}$ (see [57,58]). The symplectic diffeomorphism $\varphi$ relating canonical symplectic two-form ${\omega }_{\mathcal{Q}}=-d{\theta }_{\mathcal{Q}}$ on $T^{\ast }\mathcal{Q}$ and the symplectic two-form $\omega =-d\theta$ satisfies both

$${\varphi }^{\ast }{\theta }_{\mathcal{Q}}=\theta ,{\varphi }^{\ast }{\omega }_{\mathcal{Q}}=\omega .$$

(19)

Here, ${\theta }_{\mathcal{Q}}$ is the canonical one-form defined in (12). Accordingly, the symplectic diffeomorphism $\varphi$ can be characterized by the following pairing

$$〈\varphi \left(x\right),T\pi \circ X\left(x\right)〉=〈\theta \left(x\right),X\left(x\right)〉$$

(20)

for a vector field X on $\mathcal{P}$, for any point x in $\mathcal{P}$. To see this definition, simply evaluate X with the pull-back one form ${\varphi }^{\ast }{\theta }_{\mathcal{Q}}$ that is,

$$\begin{array}{cc}\hfill \langle {\varphi }^{\ast }{\theta }_{\mathcal{Q}},X\rangle \left(x\right)& =\langle {\theta }_{\mathcal{Q}}\left(\varphi \left(x\right)\right),T\varphi \circ X\left(x\right)\rangle \hfill \\ & =\langle {\tau }_{T^{\ast }\mathcal{Q}}\circ T\varphi \circ X\left(x\right),T{\pi }_{\mathcal{Q}}\circ T\varphi \circ X\left(x\right)\rangle \hfill \\ & =\langle \varphi \left(x\right),T\pi \circ X\left(x\right)\rangle ,\hfill \end{array}$$

(21)

where we have employed the definition (12) of the canonical one-form in the second line whereas the identities ${\tau }_{T^{\ast }\mathcal{Q}}\circ T\varphi \circ X=\varphi$ and ${\pi }_{\mathcal{Q}}\circ \varphi =\pi$ in the third line. We exhibit a special symplectic manifold $(\mathcal{P},\pi ,\mathcal{Q},\theta ,\varphi )$ with the following commutative diagram

(22)

The two-tuple $(\mathcal{P},\omega )$ is called underlying symplectic manifold of the special symplectic structure.

2.2. Morse Families

First of all, we will review the definition of a special kind of Lagrangian submanifolds of a cotangent bundle endowed with the canonical symplectic structure (for more details, see [4,59]). Let $\mathcal{N}$ be a submanifold of a smooth manifold $\mathcal{Q}$ and ${\mathsf{\Theta }|}_{\mathcal{N}}$ be a closed one-form defined on $\mathcal{N}$. Then, we can consider a Lagrangian submanifold ${\mathcal{S}}_{{\mathsf{\Theta }|}_{\mathcal{N}}}$ of $T^{\ast }\mathcal{Q}$ given by

$${\mathcal{S}}_{{\mathsf{\Theta }|}_{\mathcal{N}}}=\{\gamma \in T^{\ast }\mathcal{Q}|{\pi }_{\mathcal{Q}}\left(\gamma \right)\in \mathcal{N},\left(T_{{\pi }_{\mathcal{Q}}\left(\gamma \right)}^{\ast }\right)\iota \left(\gamma \right)=\mathsf{\Theta }{|}_{\mathcal{N}}\left({\pi }_{\mathcal{Q}}\left(\gamma \right)\right)\}$$

(23)

where $\iota :\mathcal{N}\hookrightarrow \mathcal{Q}$ is the canonical inclusion. Note that for every point q in $\mathcal{N}$, there exists an open subset $U\subseteq \mathcal{N}$ including q, and a real valued smooth function $\mathsf{\Delta }:U\mapsto \mathbb{R}$ such that ${\mathsf{\Theta }|}_{\mathcal{N}}=d\mathsf{\Delta }$. In this case, $\mathsf{\Delta }$ is called the (local) generating function of ${\mathcal{S}}_{{\mathsf{\Theta }|}_{\mathcal{N}}}$. Next, we will review the definition of Morse families and the Lagrangian submanifolds associated with them.

Consider a fiber bundle $(\mathcal{W},\tau ,\mathcal{Q})$ where $\mathcal{W}$ is the total space of dimension $n+K$, whereas $\mathcal{Q}$ is the base manifold of dimension n. Here, $\tau$ is the bundle projection. The vertical bundle $V\tau$ over the manifold $\mathcal{W}$ is a vector sub-bundle of $T\mathcal{W}$ containing vectors that belong to the kernel of the tangent mapping $T\tau$ that is,

$$V\tau =\{X\in T\mathcal{W}:T\tau \circ X=0\}.$$

(24)

The co-normal bundle $V^0\tau$ is the space of co-vectors in $T^{\ast }\mathcal{W}$ annihilating vectors in the vertical bundle $V\tau$. A real-valued function E on the total space of a fiber bundle $(\mathcal{W},\tau ,\mathcal{Q})$ is called Morse family (or generating family) [4,58] if

$$T_{\gamma }\mathit{im}\left(dE\right)+T_{\gamma }{V}^0\tau =T_{\gamma }{T}^{\ast }\mathcal{W},$$

(25)

for all $\gamma$ in the intersection $\mathit{im}\left(dE\right)\cap V^0\tau$.

Generating Lagrangian Submanifolds. A Morse family E defined on $(\mathcal{W},\tau ,\mathcal{Q})$ generates an immersed Lagrangian submanifold of the cotangent bundle $T^{\ast }\mathcal{Q}$ as

$$\mathcal{S}= \left\{z\in T^{\ast }\mathcal{Q}:T^{\ast }\tau \left(z\right)=dE\left(w\right),\forall w\in \mathcal{W},\tau \left(w\right)={\pi }_{\mathcal{Q}}\left(z\right)\right\}.$$

(26)

The inverse of this statement is also true. Thus, we state the following theorem, which is generalizing the well-known Poincaré Lemma for non-horizontal Lagrangian submanifolds. Accordingly, it is called as generalized Poincaré Lemma or Maslov-Hörmander Theorem [4,60,61].

Theorem 1.

For a Lagrangian submanifold of a cotangent bundle $T^{\ast }\mathcal{Q}$, there always exists, at least locally, a Morse family E generating $\mathcal{S}$.

We picture the Lagrangian submanifold $\mathcal{S}$ generated by a Morse family E on $(\mathcal{W},\tau ,\mathcal{Q})$ as follows.

(27)

Given a Lagrangian submanifold, its Morse family generator is far from being unique. For example, one may find a Morse family with less number of fiber variables generating the same Lagrangian submanifold. This procedure is called reduction of Morse family. See [60] for further discussions on this subject.

We merge a Morse family E defined on $(\mathcal{W},\tau ,\mathcal{Q})$ and a special symplectic structure $(\mathcal{P},\pi ,\mathcal{Q},\theta ,\varphi )$ depicted diagrammaticality in (22). This permits us to define a Lagrangian submanifold of the symplectic manifold $(\mathcal{P},\omega )$ [9]. Here is the diagram,

(28)

Here, $\mathcal{S}$ is the Lagrangian submanifold of $T^{\ast }\mathcal{Q}$ given in (26), and the inverse of the symplectic diffeomorphism $\varphi$ maps $\mathcal{S}$ to a Lagrangian submanifold ${\mathcal{S}}_E$ of $\mathcal{P}$. That is, $\varphi \left({\mathcal{S}}_E\right)=\mathcal{S}$.

Local Picture. Now, we present local realization of these discussions. Let $\left(q^i\right)$ be local coordinates on $\mathcal{Q}$ and we consider the induced local coordinates $(q^i,{ϵ}^a)$ on the the total space $\mathcal{W}$. In this image, a function E is called a Morse family if the rank of the following matrix

$$\left(\frac{{\partial }^{2}E}{\partial q^i\partial {ϵ}^a},\frac{{\partial }^{2}E}{\partial {ϵ}^a\partial {ϵ}^b}\right)$$

(29)

is equal to K. In this case, the Lagrangian submanifold (26) generated by E is viewed locally as

$$\mathcal{S}= \left\{\left(q^i,\frac{\partial E}{\partial q^j}(x,ϵ)\right) \in T^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {ϵ}^a}(x,ϵ)=0\right\}.$$

(30)

Note that the dimension of the submanifold $\mathcal{S}$ is half of the dimension of the cotangent bundle $T^{\ast }\mathcal{Q}$ and the canonical symplectic two-form ${\omega }_{\mathcal{Q}}$ vanishes on $\mathcal{S}$.

2.3. Merging Two Special Symplectic Structures

Let $(\mathcal{P},\omega )$ be an exact symplectic manifold. Assume that $\mathcal{P}$ admits two different fiber bundle structures denoted by $(\mathcal{P},\pi ,\mathcal{Q})$ and $(\mathcal{P},{\pi }^{\prime },{\mathcal{Q}}^{\prime })$. Further, let these fibrations lead to two special symplectic symplectic structures denoted by $(\mathcal{P},\pi ,\mathcal{Q},\theta ,\varphi )$ and $(\mathcal{P},{\pi }^{\prime },{\mathcal{Q}}^{\prime },{\theta }^{\prime },{\varphi }^{\prime })$, respectively. We merge these two special symplectic structures in one diagram as follows:

(31)

This is the most abstract realization of classical Tulczyjew’s triple. Notice that, in this geometry, the symplectic two-form $\omega$ on $\mathcal{P}$ admits two different potential one forms ${\theta }^{\prime }$ and $\theta$ so that

$$\omega =-d{\theta }^{\prime }=-d\theta .$$

(32)

Further, by employing the canonical symplectic two-forms ${\omega }_{\mathcal{Q}}=-d{\theta }_{\mathcal{Q}}$ and ${\omega }_{{\mathcal{Q}}^{\prime }}=-d{\theta }_{{\mathcal{Q}}^{\prime }}$ over the corresponding cotangent bundles $T^{\ast }\mathcal{Q}$ and $T^{\ast }{\mathcal{Q}}^{\prime }$, the following properties hold

$${\varphi }^{\ast }\left({\theta }_{\mathcal{Q}}\right)=\theta ,{\varphi }^{\ast }\left({\omega }_{\mathcal{Q}}\right)=\omega ,{\varphi }^{\prime \ast }\left({\theta }_{{\mathcal{Q}}^{\prime }}\right)={\theta }^{\prime },{\varphi }^{\prime \ast }\left({\omega }_{{\mathcal{Q}}^{\prime }}\right)=\omega .$$

(33)

Consider the product manifold $\mathcal{Q}\times {\mathcal{Q}}^{\prime }$ and let $\pi$ and ${\pi }^{\prime }$ be the maps defined in (31). Assume that the image space of

$$\chi :\mathcal{P}⟶\mathcal{Q}\times {\mathcal{Q}}^{\prime },p\mapsto (\pi \left(p\right),{\pi }^{\prime }\left(p\right))$$

(34)

is an embedded submanifold $\mathit{im}\chi =\mathcal{N}$ of $\mathcal{Q}\times {\mathcal{Q}}^{\prime }$. Further, we assume that the map $\chi :\mathcal{P}\mapsto \mathcal{N}$ is a surjective submersion. The equalities (32) show that the difference ${\theta }^{\prime }-\theta$ is a closed one-form. Using Poincaré Lemma, we have that for every point p in $\mathcal{P}$ there exists an open subset $U_p\subseteq \mathcal{P}$ containing p and a smooth function ${\mathsf{\Delta }}_p$ on $U_p$ such that

$$d{\mathsf{\Delta }}_p=({\theta }^{\prime }-\theta ){|}_{U_p}.$$

(35)

The vertical bundle with respect to the fibration $\chi$ is precisely the intersection of the vertical bundles with respect to $\pi$ and ${\pi }^{\prime }$, that is, $V\chi =V\pi \cap V{\pi }^{\prime }$. Notice that ${\theta }^{\prime }-\theta$ is a $\chi$-basic one-form. In fact, due to the fiber preserving character of the symplectic diffeomorphisms $\varphi$ and ${\varphi }^{\prime }$, the potential one-forms $\theta$ and ${\theta }^{\prime }$ are taking values in the co-normal bundles $V^0\pi$ and $V^0{\pi }^{\prime }$, respectively. We deduce that the difference of potential one-forms takes values in the co-normal bundle $V^0\chi$ since

$${\theta }^{\prime }-\theta \in \mathsf{\Gamma }\left(V^0{\pi }^{\prime }\right)+\mathsf{\Gamma }\left(V^0\pi \right)=\mathsf{\Gamma }(V^0{\pi }^{\prime }+V^0\pi )=\mathsf{\Gamma }{(V{\pi }^{\prime }\cap V\pi )}^0=\mathsf{\Gamma }\left(V^0\chi \right),$$

(36)

where $\chi$ is the fibration in (34). So, there exists a unique closed one-form ${\mathsf{\Theta }|}_{\mathcal{N}}$ such that ${\chi }^{\ast }{\mathsf{\Theta }|}_{\mathcal{N}}={\theta }^{\prime }-\theta$. In fact, ${\mathsf{\Delta }}_p$ given in (35) is a $\chi$-basic function and, therefore, there exists a smooth function $\mathsf{\Delta }$ on $\chi \left(U_p\right)$ satisfying $d\mathsf{\Delta }=\theta {|}_{U_p}$. The following diagram illustrates the above situation.

(37)

Next, we consider the following Lagrangian submanifold of the canonical symplectic manifold $(T^{\ast }({\mathcal{Q}}^{\prime }\times \mathcal{Q}),{\omega }_{{\mathcal{Q}}^{\prime }\times \mathcal{Q}})$ given by

$$\overline{\mathcal{S}}= \left\{\overline{\gamma }\in T^{\ast }({\mathcal{Q}}^{\prime }\times \mathcal{Q}):{\pi }_{{\mathcal{Q}}^{\prime }\times \mathcal{Q}}\left(\overline{\gamma }\right)=(q^{\prime },q)\in \mathcal{N},\overline{\gamma }{|}_{T_{(q^{\prime },q)}\mathcal{N}}{ =\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)\right\}.$$

(38)

Consider the symplectic diffeomorphism

$$\mathsf{\Psi }:\left(T^{\ast }({\mathcal{Q}}^{\prime }\times \mathcal{Q}),{\omega }_{{\mathcal{Q}}^{\prime }\times \mathcal{Q}}\right) ⟶ \left(T^{\ast }{\mathcal{Q}}^{\prime }\times T^{\ast }\mathcal{Q},{\omega }_{{\mathcal{Q}}^{\prime }}\ominus {\omega }_{\mathcal{Q}}\right),({\mu }_{q^{\prime }},{\mu }_q)\mapsto ({\mu }_{q^{\prime }},-{\mu }_q)$$

(39)

where ${\omega }_{{\mathcal{Q}}^{\prime }\times \mathcal{Q}}$ is the canonical form on the cotangent bundle, whereas

$${\omega }_{{\mathcal{Q}}^{\prime }}\ominus {\omega }_{\mathcal{Q}}={\mathit{pr}}_2^{\ast }\left({\omega }_{{\mathcal{Q}}^{\prime }}\right)-{\mathit{pr}}_1^{\ast }\left({\omega }_{\mathcal{Q}}\right).$$

Referring to (39), we map the Lagrangian submanifold $\overline{\mathcal{S}}$ in (38) to a Lagrangian submanifold $\mathcal{S}=\mathsf{\Psi }\left(\overline{\mathcal{S}}\right)$ of $T^{\ast }{\mathcal{Q}}^{\prime }\times T^{\ast }\mathcal{Q}$. Then, we will prove the following result.

Theorem 2.

The $\mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})$ of the symplectic diffeomorphism

$$\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1}:T^{\ast }{\mathcal{Q}}^{\prime }\mapsto T^{\ast }\mathcal{Q}$$

(40)

pictured in (31) is an open subset of the Lagrangian submanifold $\mathcal{S}=\mathsf{\Psi }\left(\overline{\mathcal{S}}\right)$.

Proof. 

Evidently, $\mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})$ is a Lagrangian submanifold of $T^{\ast }{\mathcal{Q}}^{\prime }\times T^{\ast }\mathcal{Q}$ so that it is enough to show that $\mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})$ is a subset of $\mathcal{S}$. In fact, we will see that

$$\begin{array}{cc}\hfill \mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})=& \{({\mu }_{q^{\prime }},-{\mu }_q)\in T^{\ast }{\mathcal{Q}}^{\prime }\times T^{\ast }\mathcal{Q}:q=\pi \left({\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)\right),\hfill \\ & ({\mu }_{q^{\prime }},-{\mu }_q){|}_{T_{(q^{\prime },q)}\mathcal{N}}{ =\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)\}\hfill \end{array}$$

which implies the result. Suppose that $({\mu }_{q^{\prime }},-{\mu }_q)$ be an element of $\mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})$, then we have that $\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)$ is precisely ${\mu }_q$. Assume that ${\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)=p$ in $\mathcal{P}$ so that $\varphi \left(p\right)=-{\mu }_q$. This gives

$${\pi }_{{\mathcal{Q}}^{\prime }}\left(\varphi \left(p\right)\right)=q\Rightarrow \pi \left(p\right)=q\Rightarrow \pi \circ {\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)=q.$$

(41)

Next, we will see that $({\mu }_{q^{\prime }},-{\mu }_q){|}_{T_{(q^{\prime },q)}\mathcal{N}}{ =\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)$. Using (20), it follows that

$$\begin{array}{cc}\hfill {\varphi }^{\prime }\left(p\right)& ={\mu }_{q^{\prime }}\Rightarrow {\theta }^{\prime }\left(p\right)=T_p^{\ast }{\pi }^{\prime }\left({\mu }_{q^{\prime }}\right)\hfill \\ \hfill \varphi \left(p\right)& =-{\mu }_q\Rightarrow -\theta \left(p\right)=T_p^{\ast }\pi \left({\mu }_q\right)\hfill \end{array}$$

(42)

therefore, one has that

$${\theta }^{\prime }\left(p\right)-\theta \left(p\right)=T_p^{\ast }\chi \circ T_{(q^{\prime },q)}^{\ast }\iota ({\mu }_{q^{\prime }},{\mu }_q).$$

(43)

However, since

$$T_p^{\ast }{\chi \circ \mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)={\theta }^{\prime }\left(p\right)-\theta \left(p\right)$$

(44)

we conclude that

$$\begin{array}{cc}\hfill T_p^{\ast }\chi \circ \mathsf{\Theta }{|}_{\mathcal{N}}(q^{\prime },q))& =T_p^{\ast }\chi \circ T_{(q^{\prime },q)}^{\ast }\iota ({\mu }_{q^{\prime }},{\mu }_q),\hfill \\ \hfill {\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)& =T_{(y^{\prime },y)}^{\ast }\iota ({\mu }_{q^{\prime }},{\mu }_q)=({\mu }_{q^{\prime }},{\mu }_q){|}_{T_{(q^{\prime },q)}\mathcal{N}}.\hfill \end{array}$$

(45)

Conversely, suppose that $\gamma =({\mu }_{q^{\prime }},-{\mu }_q)$ in $T^{\ast }{\mathcal{Q}}^{\prime }\times T^{\ast }\mathcal{Q}$ and $q=\pi \left({\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)\right)$, $({\mu }_{q^{\prime }},-{\mu }_q){|}_{T_{(q^{\prime },q)}\mathcal{N}}{ =\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)$. If $p={\left({\varphi }^{\prime }\right)}^{-1}\left({\mu }_{q^{\prime }}\right)$, we will see that $\varphi \left(p\right)=-{\mu }_q$, which implies that $({\mu }_{q^{\prime }},-{\mu }_q)$ in $\mathrm{graph}(\varphi \circ {\left({\varphi }^{\prime }\right)}^{-1})$. We see that ${\varphi }^{\prime }\left(p\right)={\mu }_{q^{\prime }}$ which gives ${\theta }^{\prime }\left(p\right)=T_p^{\ast }{\pi }^{\prime }\left({\mu }_{q^{\prime }}\right)$. On the other hand,

$$\begin{array}{cc}\hfill T_{(q^{\prime },q)}^{\ast }\iota ({\mu }_{q^{\prime }},{\mu }_q)& {=\mathsf{\Theta }|}_{\mathcal{N}}(q^{\prime },q)\hfill \\ \hfill T_p^{\ast }\chi \circ T_{(q^{\prime },q)}^{\ast }\iota ({\mu }_{q^{\prime }},{\mu }_q)& ={\theta }^{\prime }\left(p\right)-\theta \left(p\right)\hfill \\ \hfill T^{\ast }(\iota \circ \chi )({\mu }_{q^{\prime }},{\mu }_q)& ={\theta }^{\prime }\left(p\right)-\theta \left(p\right)\hfill \\ \hfill T_p^{\ast }{\pi }^{\prime }\left({\mu }_{q^{\prime }}\right)+T_p^{\ast }\pi \left({\mu }_q\right)={\theta }^{\prime }\left(p\right)-\theta \left(p\right).\end{array}$$

(46)

Thus, as $T_p^{\ast }{\pi }^{\prime }\left({\mu }_{q^{\prime }}\right)={\theta }^{\prime }\left(p\right)$, we conclude that

$$T_p^{\ast }\pi (-{\mu }_q)=\theta \left(p\right)\Rightarrow \pi \left(p\right)=-{\mu }_q.$$

(47)

□

To illustrate the Legendre transformation in this geometry, we present the following discussion. Assume the existence of the triple in (31). Start with the left wing of the triple by considering a real valued function F on $\mathcal{Q}$. This determines a Lagrangian submanifold $\mathit{im}\left(dF\right)$ of $T^{\ast }\mathcal{Q}$. By employing the inverse of the symplectic diffeomorphism $\varphi$ to $\mathit{im}\left(dF\right)$, one arrives at a Lagrangian submanifold ${\mathcal{S}}_F$ of the symplectic manifold $\mathcal{P}$. This submanifold can also be defined in terms of the potential one-form $\theta$ as

$${\mathcal{S}}_F=\{p\in \mathcal{P}:d(F\circ \pi )\left(p\right)=\theta \left(p\right)\}.$$

The Legendre transformation in terms of Tulczyjew [9,12] is to determine a generating family for ${\mathcal{S}}_{\mathcal{F}}$ referring to the right wing of the triple (31). To have that, first map ${\mathcal{S}}_F$ to a Lagrangian submanifold of $T^{\ast }{\mathcal{Q}}^{\prime }$ by means of the symplectic diffeomorphism ${\varphi }^{\prime }$. As a manifestation of the generalized Poincaré lemma, there exits a Morse family $E^{\prime }$ on a fiber bundle $(\mathcal{W},\tau ,\mathcal{Q})$ generating $\varphi \left({\mathcal{S}}_F\right)$.

2.4. The Classical Tulczyjew’s Triple

In this subsection, we draw Tulczyjew’s triple for classical dynamics assuming a configuration manifold $\mathcal{Q}$. This is to construct the triple (31) by replacing ${\mathcal{Q}}^{\prime }$ with the tangent bundle $T\mathcal{Q}$, whereas replacing $\mathcal{Q}$ with the cotangent bundle $T^{\ast }\mathcal{Q}$. On the upper level of (31), this results with the iterated bundles $T^{\ast }T\mathcal{Q}$, $T{T}^{\ast }\mathcal{Q}$ and $T^{\ast }{T}^{\ast }\mathcal{Q}$ in order. See that, being cotangent bundles, both $T^{\ast }T\mathcal{Q}$ and $T^{\ast }{T}^{\ast }\mathcal{Q}$ are symplectic. We now establish the symplectic structure on $T{T}^{\ast }\mathcal{Q}$ admitting two potential one forms.

Symplectic structure on $T{T}^{\ast }\mathcal{Q}$. Consider the canonical symplectic manifold $T^{\ast }Q$ equipped with the exact symplectic two-form ${\omega }_{\mathcal{Q}}=-d{\theta }_{\mathcal{Q}}$. The derivation $i_T$ takes the symplectic two-form ${\omega }_{\mathcal{Q}}$ on $T^{\ast }\mathcal{Q}$ to a one-form on $T{T}^{\ast }\mathcal{Q}$ as

$$i_T{\omega }_{\mathcal{Q}}\left(Y\right)={\omega }_{\mathcal{Q}}({\tau }_{T{T}^{\ast }\mathcal{Q}}\left(Y\right),T{\tau }_{T^{\ast }Q}\left(Y\right)),$$

for any tangent vector Y on $T{T}^{\ast }\mathcal{Q}$. Here, ${\tau }_{T{T}^{\ast }\mathcal{Q}}$ is the tangent bundle projection $TT{T}^{\ast }\mathcal{Q}$ to $T{T}^{\ast }\mathcal{Q}$ whereas $T{\tau }_{T^{\ast }\mathcal{Q}}$ is the tangent mapping of the bundle projection ${\tau }_{T^{\ast }\mathcal{Q}}:T{T}^{\ast }\mathcal{Q}\mapsto T^{\ast }\mathcal{Q}$. We define two one-forms on $T{T}^{\ast }\mathcal{Q}$ as

$${\vartheta }_1=-i_T{\omega }_{\mathcal{Q}},{\vartheta }_2=d_T{\theta }_{\mathcal{Q}}=i_{T}d{\theta }_Q+d{i}_T{\theta }_{\mathcal{Q}}.$$

(48)

where the derivation $d_T$ is the commutator $i_{T}d+d{i}_T$. The exterior derivatives of these one-forms results with a symplectic two-form on $T{T}^{\ast }\mathcal{Q}$ defined to be

$$d_T{\omega }_{\mathcal{Q}}=-d{\vartheta }_1=-d{\vartheta }_2.$$

(49)

We record this in the following theorem [12].

Theorem 3.

The tangent bundle $T{T}^{\ast }\mathcal{Q}$ is a symplectic manifold with the symplectic two-form $d_T{\omega }_{\mathcal{Q}}$ in (49) admitting two potential one-forms given in (48).

Consider the Darboux’s coordinates $(q^i,p_i)$ on the cotangent bundle $T^{\ast }\mathcal{Q}$. In terms of the induced local coordinate chart $(q^i,p_i;{\dot{q}}^i,{\dot{p}}_i)$ on the tangent bundle $T{T}^{\ast }\mathcal{Q}$, the potential one-forms in (48) are computed to be

$${\vartheta }_1=-i_T{\omega }_{\mathcal{Q}}={\dot{p}}_{i}d{q}^i-{\dot{q}}^{i}d{p}_i,{\vartheta }_2=d_T{\theta }_{\mathcal{Q}}={\dot{p}}_{i}d{q}^i+p_{i}d{\dot{q}}^i.$$

(50)

Notice that, in this case, the symplectic two-form turns out to be

$$d_T{\omega }_{\mathcal{Q}}=d{q}^i\wedge d{\dot{p}}_i+d{\dot{q}}^i\wedge d{p}_i.$$

(51)

Note that, the value ${\vartheta }_2-{\vartheta }_1$ is an exact one-form. Actually, it is the exterior derivative of coupling function $i_T{\theta }_{\mathcal{Q}}:T{T}^{\ast }\mathcal{Q}\mapsto \mathbb{R}$.

One can arrive at the tangent bundle symplectic two-form $d_T{\omega }_{\mathcal{Q}}$ on $T{T}^{\ast }\mathcal{Q}$ as the complete lift of the canonical symplectic two-form ${\omega }_{\mathcal{Q}}$ on $T^{\ast }\mathcal{Q}$. To have this, from [3,62] we recall the definition of the complete lift of a differential form. The complete lift of the canonical one-form ${\theta }_{\mathcal{Q}}=p_{i}d{q}^i$ is computed to be

$${\theta }_{\mathcal{Q}}^C={\dot{p}}_{i}d{q}^i+p_{i}d{\dot{q}}^i.$$

(52)

Complete lift of forms commutes with exterior derivative. Thus, we compute

$${\omega }_{\mathcal{Q}}^C={(-d{\theta }_{\mathcal{Q}})}^C=d{q}^i\wedge d{\dot{p}}_i+d{\dot{q}}^i\wedge d{p}_i=d_T{\omega }_{\mathcal{Q}}.$$

(53)

As manifested in the display, we conclude that ${\omega }_{\mathcal{Q}}^C=d_T{\omega }_{\mathcal{Q}}$.

By recalling Diagram (13), notice that $T{T}^{\ast }\mathcal{Q}$ admits two bundle structures. It is a vector bundle over $T\mathcal{Q}$ with respect to the vector bundle projection $T{\pi }_{\mathcal{Q}}$, and it is a vector bundle over $T^{\ast }\mathcal{Q}$ with vector bundle projection ${\tau }_{T^{\ast }\mathcal{Q}}$. Thus, we can consider the projection

$$\sigma :T{T}^{\ast }\mathcal{Q}⟶\mathfrak{N}=T\mathcal{Q} {\times }_{\mathcal{Q}}{T}^{\ast }\mathcal{Q},Z\mapsto (T{\pi }_{\mathcal{Q}}\left(Z\right),{\tau }_{T^{\ast }\mathcal{Q}}\left(Z\right))$$

(54)

from $T{T}^{\ast }\mathcal{Q}$ to the Whitney sum $\mathfrak{N}=T\mathcal{Q} {\times }_{\mathcal{Q}}{T}^{\ast }\mathcal{Q}$. The function $i_T{\theta }_{\mathcal{Q}}$ is $\sigma$-basic and it induces a smooth function $\mathsf{\Delta }$ on $\mathfrak{N}$, which is just the coupling function

$$\mathsf{\Delta }={\dot{q}}^i{p}_i.$$

(55)

In addition, $\mathfrak{N}$ is a submanifold of the product manifold $T\mathcal{Q}\times T^{\ast }\mathcal{Q}$. Thus, we can consider the Lagrangian submanifold ${\mathcal{S}}_{{\theta }_{\mathfrak{N}}}$ induced by the exact one-form ${\theta }_{\mathfrak{N}}=d\mathsf{\Delta }$ as in Section 2.2. Moreover, using the construction presented in Section 2.3, we deduce that ${\mathcal{S}}_{{\theta }_{\mathfrak{N}}}$ is the graph of a symplectic diffeomorphism $\psi$ between $T^{\ast }T\mathcal{Q}$ and $T^{\ast }{T}^{\ast }\mathcal{Q}$. In terms of the Darboux’ coordinates $(q^i,{\dot{q}}^i,a_i,{\dot{a}}_i)$ on $T^{\ast }T\mathcal{Q}$, it is computed to be

$$\psi :T^{\ast }T\mathcal{Q}⟶T^{\ast }{T}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^i,a_i,{\dot{a}}_i)\mapsto (q^i,{\dot{a}}_i,a_i,-{\dot{q}}^i).$$

(56)

The canonical involution on $TT\mathcal{Q}$. The iterated tangent bundle $TT\mathcal{Q}$ is a tangent bundle with the base manifold $T\mathcal{Q}$ along with the tangent bundle projection ${\tau }_{T\mathcal{Q}}$. It is possible to show that $TT\mathcal{Q}$ can be written as a vector bundle over $T\mathcal{Q}$ apart from the canonical tangent bundle fibration. This fibration is achieved by the tangent mapping $T{\tau }_{\mathcal{Q}}$ of the projection ${\tau }_{\mathcal{Q}}$. To manifest this, we plot here the double vector bundle structure of the iterated tangent bundle $TT\mathcal{Q}$ in the following commutative diagram

(57)

Referring to this double bundle structure [1,63,64], consider a differential mapping $\mathsf{\Gamma }=\mathsf{\Gamma }(t,s)$ from an open domain of ${\mathbb{R}}^2$ to $\mathcal{Q}$. The differential of $\mathsf{\Gamma }$ with respect to t, at $t=0$, results with a curve $\gamma \left(s\right)$ lying in the tangent bundle $T\mathcal{Q}$ depending on the free variable s. If one further takes the differential of the curve $\gamma \left(s\right)\subset T\mathcal{Q}$ with respect to s, at $s=0$, then arrives at a vector in $TT\mathcal{Q}$. Accordingly, the canonical involution ${\kappa }_{\mathcal{Q}}$ on $TT\mathcal{Q}$ is defined by changing the order of differentiations

$$\kappa :TT\mathcal{Q}⟶TT\mathcal{Q}:{\left.\frac{d}{dt}\right|}_{t=0}{\left.\frac{d}{ds}\right|}_{s=0}\mathsf{\Gamma }(t,s)⟶{\left.\frac{d}{ds}\right|}_{s=0}{\left.\frac{d}{dt}\right|}_{t=0}\mathsf{\Gamma }(t,s).$$

(58)

An observation gives that the involution $\kappa$ is changing the order of the fibrations

$${\tau }_{T\mathcal{Q}}\circ \kappa =T{\tau }_{\mathcal{Q}},T{\tau }_{\mathcal{Q}}\circ \kappa ={\tau }_{T\mathcal{Q}}.$$

(59)

Pairing between $T{T}^{\ast }\mathcal{Q}$ and $TT\mathcal{Q}$. We now establish a pairing between an element Z in $T{T}^{\ast }\mathcal{Q}$ and an element W in $TT\mathcal{Q}$ such that $T{\pi }_{\mathcal{Q}}\left(Z\right)=T{\tau }_{\mathcal{Q}}\left(W\right)$ and that ${\pi }_{\mathcal{Q}}\circ z={\tau }_{\mathcal{Q}}\circ v$. Recall that, there is a curve $z\left(t\right)$ in $T^{\ast }\mathcal{Q}$ so that $Z=\dot{z}\left(0\right)$ and, there is a curve $v\left(t\right)$ in $T\mathcal{Q}$ so that $W=\dot{v}\left(0\right)$. In this framework, the pairing is defined as

$${\langle •,•\rangle }^{\tilde{}}:T{T}^{\ast }\mathcal{Q}\times TT\mathcal{Q}⟶\mathbb{R},{\langle Z,W\rangle }^{\tilde{}}={\left.\frac{d}{dt}\langle z\left(t\right),v\left(t\right)\rangle \right|}_{t=0}.$$

(60)

Here, the pairing on the right hand side is the one between $T^{\ast }\mathcal{Q}$ and $T\mathcal{Q}$.

In the induced coordinates $(q^i,{\dot{q}}^i,q^{\prime i},{\dot{q}}^{\prime i})$ on the iterated tangent bundle $TT\mathcal{Q}$, the fibration in (57) read

$${\tau }_{T\mathcal{Q}}(q^i,{\dot{q}}^i,q^{\prime i},{\dot{q}}^{\prime i})=(q^i,{\dot{q}}^i),T{\tau }_{\mathcal{Q}}(q^i,{\dot{q}}^i,q^{\prime i},{\dot{q}}^{\prime i})=(q^i,q^{\prime i}),$$

(61)

whereas the canonical involution $\kappa$ in (59) is computed to be

$$\kappa (q^i,{\dot{q}}^i,q^{\prime i},{\dot{q}}^{\prime i})=(q^i,q^{\prime i},{\dot{q}}^i,{\dot{q}}^{\prime i}).$$

Coordinate expression of the pairing (60) is as follows. Let us choose coordinates on $T{T}^{\ast }\mathcal{Q}$ as $Z=(q^i,p_i,{\dot{q}}^i,{\dot{p}}_i)$, and coordinates $W=(q^i,q^{\prime i},{\dot{q}}^i,{\dot{q}}^{\prime i})$ on $TT\mathcal{Q}$ then,

$${\langle Z,W\rangle }^{\tilde{}}={〈(q^i,p_i,{\dot{q}}^i,{\dot{p}}_i),(q^i,q^{\prime i},{\dot{q}}^i,{\dot{q}}^{\prime i})〉}^{\tilde{}}=p_i{\dot{q}}^{\prime i}+q^{\prime i}{\dot{p}}_i.$$

Now, we define two special symplectic structures for the symplectic manifold $(T{T}^{\ast }\mathcal{Q},d_T{\omega }_{\mathcal{Q}})$. One is to $T^{\ast }T\mathcal{Q}$ and the other is to $T^{\ast }{T}^{\ast }\mathcal{Q}$. We assume that, being cotangent bundles, $T^{\ast }T\mathcal{Q}$ and $T^{\ast }{T}^{\ast }\mathcal{Q}$ are equipped with the canonical symplectic forms ${\omega }_{T\mathcal{Q}}=-d{\theta }_{T\mathcal{Q}}$ and ${\omega }_{T^{\ast }\mathcal{Q}}=-d{\theta }_{T^{\ast }\mathcal{Q}}$, respectively.

Left Wing of the Tulczyjew’s triple. Start with defining the vector fibration morphism

$$\alpha :T{T}^{\ast }\mathcal{Q}⟶T^{\ast }T\mathcal{Q},\langle \alpha \left(Z\right),W\rangle ={\langle Z,\kappa \left(W\right)\rangle }^{\tilde{}},$$

(62)

where $\kappa$ is the canonical involution defined in (59), whereas the pairing on the right hand side is the one given in (60). Here, the pairing of the left hand side is the canonical pairing between $T^{\ast }T\mathcal{Q}$ and $TT\mathcal{Q}$. Locally, one has that

$$\alpha (q^i,p_i,{\dot{q}}^i,{\dot{p}}_i)=(q^i,{\dot{q}}^i,{\dot{p}}_i,p_i).$$

(63)

See that, $\alpha$ is a symplectic diffeomorphism by satisfying ${\alpha }^{\ast }{\omega }_{T\mathcal{Q}}=d_T{\omega }_{\mathcal{Q}}$. Here, ${\omega }_{T\mathcal{Q}}$ is the canonical symplectic two-form on $T^{\ast }T\mathcal{Q}$, and $d_T{\omega }_{\mathcal{Q}}$ is the lifted symplectic two-form (49) on $T{T}^{\ast }\mathcal{Q}$. So that we arrive at a special symplectic structure

$$(T{T}^{\ast }Q,T{\pi }_{\mathcal{Q}},T\mathcal{Q},{\vartheta }_2,\alpha )$$

(64)

where ${\vartheta }_2$ is the one-form in (50), and $T{\pi }_{\mathcal{Q}}$ is the tangent lift of the cotangent bundle projection ${\pi }_{\mathcal{Q}}$. We include this in the following diagram for future reference

(65)

Right Wing of the Tulczyjew’s triple. The nondegeneracy of the canonical symplectic two-form ${\omega }_{\mathcal{Q}}$ leads to the existence of a (musical) diffeomorphism

$$\beta :T{T}^{\ast }\mathcal{Q}⟶T^{\ast }{T}^{\ast }\mathcal{Q},Z\mapsto -{\omega }_{\mathcal{Q}}^{♭}\left(Z\right)=-{\omega }_{\mathcal{Q}}(Z,•).$$

(66)

Locally, one has that

$$\beta (q^i,p_i,{\dot{q}}^i,{\dot{p}}_i)=(q^i,p_i,{\dot{p}}_i,-{\dot{q}}^i).$$

(67)

Note that, $\beta$ is a symplectomorphism by satisfying ${\beta }^{\ast }{\omega }_{T^{\ast }\mathcal{Q}}=d_T{\omega }_{\mathcal{Q}}$. Here, ${\omega }_{T^{\ast }\mathcal{Q}}$ is the canonical symplectic two-form on $T^{\ast }{T}^{\ast }\mathcal{Q}$. So that we arrive at the following special symplectic structure

$$(T{T}^{\ast }Q,{\tau }_{T^{\ast }Q},T^{\ast }Q,{\vartheta }_1,\beta ).$$

(68)

Here, ${\vartheta }_1$ is the potential one-form in (48) and that ${\tau }_{T^{\ast }Q}$ is the tangent bundle projection. The diagram is for the future reference.

(69)

The Tulczyjew’s triple for classical mechanics is the following commutative diagram merging the two special symplectic structures (65) and (69) in one picture

(70)

where $\beta$ and $\alpha$ are symplectomorphisms such that

$${\beta }^{\ast }{\omega }_{T^{\ast }\mathcal{Q}}={\omega }_{\mathcal{Q}}^C=-d{\vartheta }_1,{\alpha }^{\ast }{\omega }_{T\mathcal{Q}}={\omega }_{\mathcal{Q}}^C=-d{\vartheta }_2.$$

By employing the composition of two symplectic diffeomorphisms ${\alpha }^{-1}$ and $\beta$, we arrive at a symplectic diffeomorphism $\psi =\beta \circ {\alpha }^{-1}$. This is exactly the one in (56).

2.5. The Legendre Transformation

Let us now realized Euler–Lagrange equations generated by a Lagrangian function L on the tangent $T\mathcal{Q}$ as a Lagrangian submanifold of $T{T}^{\ast }\mathcal{Q}$ using the left side of Tulczyjew’s triple (70). The image space of the exterior derivative $dL$ is a Lagrangian submanifold of the manifold $T^{\ast }T\mathcal{Q}$. By employing the inverse of $\alpha$, we transfer this Lagrangian submanifold to a Lagrangian submanifold ${\mathcal{S}}_L$ of the manifold $T{T}^{\ast }\mathcal{Q}$. In terms of the local coordinates $(q^i,p_i,{\dot{q}}^i,{\dot{p}}_i)$ on $T{T}^{\ast }\mathcal{Q}$, $\mathcal{S}$ is computed to be

$${\mathcal{S}}_L=\alpha \left(\mathit{im}dL\right)= \left\{\left(q^i,\frac{\partial L}{\partial {\dot{q}}^i},{\dot{q}}^i,\frac{\partial L}{\partial q^i}\right) \in T{T}^{\ast }\mathcal{Q}:L=L(q^i,{\dot{q}}^i)\right\} \subset T{T}^{\ast }\mathcal{Q}.$$

(71)

Dynamics determined by the Lagrangian submanifold ${\mathcal{S}}_L$ is computed simply by taking the time derivatives of the base coordinates $\left(q^i,\partial L/\partial {\dot{q}}^i\right)$ and equate them to the fiber coordinates given by $({\dot{q}}^i,\partial L/\partial q^i)$, respectively. By this, we arrive at the Euler–Lagrange equations

$$\frac{d{q}^i}{dt}={\dot{q}}^i,\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}^i}\right) =\frac{\partial L}{\partial q^i}.$$

(72)

Notice that, for this realization, we have not asked any regularity conditions for the Lagrangian function. That is, this theory is valid for singular Lagrangians as well.

The Legendre Transformation of Euler–Lagrange Equations. We now generate the Lagrangian submanifold ${\mathcal{S}}_L$ given in (71) referring to the right wing of the Tulczyjew’s triple (70). Let us remark once more that this is the Legendre transformation in the understanding of Tulczyjew. To have that, at first consider the Morse family over the Pontryagin bundle as

$$E:T\mathcal{Q}{\times }_{\mathcal{Q}}{T}^{\ast }\mathcal{Q}⟶\mathbb{R},(u,\varsigma )\mapsto \langle \varsigma ,u\rangle -L\left(u\right)$$

(73)

where $\langle •,•\rangle$ is the canonical pairing between $T^{\ast }\mathcal{Q}$ and $T\mathcal{Q}$. In coordinates, we write the Morse family as

$$E(q^i,{\dot{q}}^i,p_i)=p_i{\dot{q}}^i-L(q,\dot{q}).$$

(74)

We remark that the matrix

$$\left(\frac{{\partial }^{2}E}{\partial {\dot{q}}^i\partial q^j},\frac{{\partial }^{2}E}{\partial {\dot{q}}^i\partial p_j},\frac{{\partial }^{2}E}{\partial {\dot{q}}^i\partial {\dot{q}}^j}\right)$$

(75)

has maximal rank. For the choice of sign conventions, we consider the minus of the Morse family, that is, $-E$. So, according to the definition in (30), the minus of the Morse family $-E$ generates a Lagrangian submanifold of the cotangent bundle $T^{\ast }{T}^{\ast }\mathcal{Q}$ as

$$\mathcal{S}= \left\{\left(q^i,p_i,-\frac{\partial E}{\partial q^i},-\frac{\partial E}{\partial p_i}\right) \in T^{\ast }{T}^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {\dot{q}}^i}=0\right\} \subset T^{\ast }{T}^{\ast }\mathcal{Q}.$$

(76)

In order to transfer this Lagrangian submanifold to $T{T}^{\ast }\mathcal{Q}$, we merge the Pontryagin bundle with the right wing of the Tulczyjew’s triple (70) as follows

(77)

This permits us to arrive at the following Lagrangian submanifold of $T{T}^{\ast }\mathcal{Q}$

$${\mathcal{S}}_{-E}={\beta }^{-1}\left(\mathcal{S}\right)= \left\{\left(q^i,p_i,\frac{\partial E}{\partial p_i},-\frac{\partial E}{\partial q^i}\right) \in T{T}^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {\dot{q}}^i}=0\right\} \subset T{T}^{\ast }\mathcal{Q}.$$

(78)

A direct computation proves that the Lagrangian submanifold ${\mathcal{S}}_{-E}$ in (78) and the Lagrangian submanifold ${\mathcal{S}}_L$ in (71) are the same. So that the Legendre transformation is achieved. If the Lagrangian function is non-degenerate, then from the equation

$$\frac{\partial E}{\partial {\dot{q}}^i}(q,\dot{q},p)=p_i-\frac{\partial L}{\partial {\dot{q}}^i}(q,\dot{q})=0$$

(79)

one can explicitly determine the velocity ${\dot{q}}^i$ in terms of the momenta $(q^i,p_i)$. In other words, for a non-degenerate Lagrangian function $L=L(q,\dot{q})$ the fiber derivative

$$\mathbb{F}L:T\mathcal{Q}⟶T^{\ast }\mathcal{Q},(q^i,{\dot{q}}^j)⟶ \left(q^i,\frac{\partial L}{\partial {\dot{q}}^i}(q,\dot{q})\right)$$

(80)

is a local diffeomorphism. In this case, the Morse family E can be reduced to a well-defined Hamiltonian function

$$H(q^i,p_i)=p_i{\dot{q}}^i(q,p)-L\left(q,\dot{q}(q,p)\right)$$

(81)

on $T^{\ast }\mathcal{Q}$.

Inverse Legendre Transformation. The inverse Legendre transformation is also possible in a similar way. This time, one starts with a Hamiltonian system $(T^{\ast }\mathcal{Q},{\omega }_{\mathcal{Q}},H)$ where H is a Hamiltonian function. See that, in this notation Hamiltonian vector field $X_H$ defined in (15) is determined through

$$\beta \circ X_H=-dH.$$

(82)

Notice that the Lagrangian submanifold determined by the equality (82) is written in coordinates as

$${\mathcal{S}}_{-H}= \left\{\left(q^i,p_i,\frac{\partial H}{\partial p_i},-\frac{\partial H}{\partial q^i}\right) \in T{T}^{\ast }\mathcal{Q}\right\} \subset T{T}^{\ast }\mathcal{Q}.$$

(83)

Evidently, this Lagrangian submanifold is precisely determining the Hamilton’s Equation (17). In the present picture, the inverse Legendre transformation is to generate the Lagrangian submanifold (83) by referring to the right wing of the triple. If the Hamiltonian function is not regular then one needs to employ a Morse family

$$F:T\mathcal{Q}{\times }_{\mathcal{Q}}{T}^{\ast }\mathcal{Q}⟶\mathbb{R},(u,\varsigma )\mapsto \langle \varsigma ,u\rangle -H\left(\varsigma \right).$$

(84)

So, if we consider the Pontryagin bundle over $T\mathcal{Q}$ and we proceed as in the previous subsection, we will obtain the inverse Legendre transformation.

3. Contact Dynamics

3.1. Contact Manifolds

A $(2n+1)$—dimensional manifold $\mathcal{M}$ is called contact manifold if it is equipped with a contact one-form $\eta$ satisfying $d{\eta }^n\wedge \eta \ne 0$ [4,34]. We denote a contact manifold by a two-tuple $(\mathcal{M},\eta )$. The Reeb vector field $\mathcal{R}$ is the unique vector field satisfying

$${\iota }_{\mathcal{R}}\eta =1,{\iota }_{\mathcal{R}}d\eta =0.$$

(85)

At each point of the manifold $\mathcal{M}$, the kernel of the contact form $\eta$ determines the contact structure $H\mathcal{M}$. The complement of this structure, denoted by $V\mathcal{M}$, is determined by the kernel of the exact two-form $d\eta$. These give the following decomposition of the tangent bundle

$$T\mathcal{M}=H\mathcal{M}\oplus V\mathcal{M},H\mathcal{M}=ker\eta ,V\mathcal{M}=kerd\eta .$$

(86)

Here, $H\mathcal{M}$ is a vector sub-bundle of rank $2n$. The restriction of $d\eta$ to $H\mathcal{M}$ is non-degenerate so that $(H\mathcal{M},d\eta )$ is a symplectic vector bundle over $\mathcal{M}$. The rank of $V\mathcal{M}$ is 1, and it is generated by the Reeb field $\mathcal{R}$.

Contactization. It is possible to arrive at a contact manifold starting from a symplectic manifold. To have this, consider a symplectic manifold $\mathcal{P}$ admitting an integer symplectic two form $\omega$. Introduce the principal circle (quantization) bundle

$$S^1⇝(\mathcal{M},\eta )\stackrel{\mathit{pr}}{⟶}(\mathcal{P},\omega ).$$

(87)

The contact one-form on $\mathcal{M}$ is the connection one-form associated with a principal connection on the principal $S^1$-bundle $pr:\mathcal{M}\mapsto \mathcal{P}$ with curvature $\omega$. This procedure is called contactization.

Another example of a contact manifold can be obtained from an exact symplectic manifold as follows. Consider a trivial line bundle over a manifold given by $\mathcal{Q}\times \mathbb{R}\mapsto \mathcal{Q}$. The first jet bundle, denoted by ${\mathcal{T}}^{\ast }\mathcal{Q}$ is diffeomorphic to the product space $T^{\ast }\mathcal{Q}\times \mathbb{R}$ that is,

$${\mathcal{T}}^{\ast }\mathcal{Q}=T^{\ast }\mathcal{Q}\times \mathbb{R}.$$

(88)

We call this space the extended cotangent bundle. There exist two projections

$$\begin{array}{cc}\hfill {\pi }_{\mathcal{Q}}^1& :{\mathcal{T}}^{\ast }\mathcal{Q}=T^{\ast }\mathcal{Q}\times \mathbb{R}⟶T^{\ast }\mathcal{Q},(\zeta ,z)\mapsto \zeta \hfill \\ \hfill {\pi }_{\mathcal{Q}}^0& :{\mathcal{T}}^{\ast }\mathcal{Q}=T^{\ast }\mathcal{Q}\times \mathbb{R}⟶\mathcal{Q},(\zeta ,z)\mapsto {\pi }_{\mathcal{Q}}\left(\zeta \right),\hfill \end{array}$$

(89)

where ${\pi }_{\mathcal{Q}}$ is the cotangent bundle projection whereas z is the standard coordinate on $\mathbb{R}$. Referring to the fibration defined by ${\pi }_{\mathcal{Q}}^1$, we have the following globally trivial contactization of the canonical symplectic manifold

$$\mathbb{R}⇝({\mathcal{T}}^{\ast }\mathcal{Q},{\eta }_{\mathcal{Q}})\stackrel{{\pi }_{\mathcal{Q}}^1}{⟶}(T^{\ast }\mathcal{Q},{\omega }_{\mathcal{Q}}).$$

(90)

Here, the contact one-form on the jet bundle ${\mathcal{T}}^{\ast }\mathcal{Q}$ is defined to be

$${\eta }_{\mathcal{Q}}:=dz-{\theta }_{\mathcal{Q}}$$

(91)

where ${\theta }_{\mathcal{Q}}$ is the canonical one-form (12) on the cotangent bundle $T^{\ast }\mathcal{Q}$. Notice that, we have employed abuse of notation by identifying z and ${\theta }_{\mathcal{Q}}$ with their pull-backs on the total space ${\mathcal{T}}^{\ast }\mathcal{Q}$. The previous construction also works if we replace $T^{\ast }\mathcal{Q}$ by an arbitrary exact symplectic manifold $\mathcal{P}$ and, in such a case, we obtain a contact structure on the product manifold $\mathcal{P}\times \mathbb{R}$. There exist Darboux’ coordinates $(q^i,p_i,z)$ on ${\mathcal{T}}^{\ast }\mathcal{Q}$, where i is running from 1 to n. In these coordinates, the contact one-form and the Reeb vector field are computed to be

$${\eta }_{\mathcal{Q}}=dz-p_{i}d{q}^i,\mathcal{R}=\frac{\partial }{\partial z},$$

(92)

respectively. Notice that, in this realization, the horizontal bundle is generated by the vector fields

$$H{\mathcal{T}}^{\ast }\mathcal{Q}=span\{{\xi }_i,{\xi }^i\},{\xi }_i= \frac{\partial }{\partial q^i}+ p_i\frac{\partial }{\partial z},{\xi }^i=\frac{\partial }{\partial p_i}.$$

(93)

It is important to note that these generators are not closed under the Jacobi–Lie bracket that is,

$$[{\xi }^i,{\xi }_j]={\delta }_j^i\mathcal{R},$$

(94)

where ${\delta }_j^i$ stands for the Kronecker delta. The Darboux’s theorem manifests that local picture presented in this subsection is generic for all contact manifolds of dimension $2n+1$.

Musical Mappings. For a contact manifold $(\mathcal{M},\eta )$, there is a musical isomorphism ♭ from the tangent bundle $T\mathcal{M}$ to the cotangent bundle $T^{\ast }\mathcal{M}$ defined to be

$$♭:T\mathcal{M}⟶T^{\ast }\mathcal{M},v\mapsto {\iota }_{v}d\eta +\eta \left(v\right)\eta .$$

(95)

This mapping takes the Reeb field $\mathcal{R}$ to the contact one-form $\eta$. We denote the inverse of this mapping by ♯. Referring to this, we define a bi-vector field $\mathsf{\Lambda }$ on $\mathcal{M}$ as

$$\mathsf{\Lambda }(\alpha ,\beta )=-d\eta (♯\alpha ,♯\beta ).$$

(96)

The couple $(\mathsf{\Lambda },-\mathcal{R})$ induces a Jacobi structure [59,65,66]. This is a manifestation of the equalities

$$[\mathsf{\Lambda },\mathsf{\Lambda }]=-2\mathcal{R}\wedge \mathsf{\Lambda },[\mathcal{R},\mathsf{\Lambda }]=0,$$

(97)

where the bracket is the Schouten–Nijenhuis bracket. We cite [4,66,67,68] for more details on the Jacobi structure associated with a contact one-form. Referring to the bi-vector field $\mathsf{\Lambda }$ we introduce the following musical mapping

$${♯}_{\mathsf{\Lambda }}:T^{\ast }\mathcal{M}⟶T\mathcal{M},\alpha \mapsto \mathsf{\Lambda }(\alpha ,•)=♯\alpha -\alpha \left(\mathcal{R}\right)\mathcal{R}.$$

(98)

Evidently, the mapping ${♯}_{\mathsf{\Lambda }}$ fails to be an isomorphism. Notice that the kernel is spanned by the contact one-form $\eta$. So that, the image space of ${♯}_{\mathsf{\Lambda }}$ is precisely the horizontal bundle $H\mathcal{M}$ exhibited in (93).

In terms of the Darboux coordinates $(q^i,p_i,z)$, we compute the image of a one-form in $T^{\ast }\mathcal{M}$ by ${♯}_{\mathsf{\Lambda }}$ as

$${♯}_{\mathsf{\Lambda }}:{\alpha }_{i}d{q}^i+{\alpha }^{i}d{p}^i+udz\mapsto {\alpha }^i\frac{\partial }{\partial q^i}-({\alpha }_i+p_{i}u)\frac{\partial }{\partial p_i}+{\alpha }^i{p}_i\frac{\partial }{\partial z}.$$

(99)

Symplectization. The symplectization of a contact manifold $(\mathcal{M},\eta )$ is the symplectic manifold $(\mathcal{M}\times \mathbb{R},d\left(e^t\eta \right))$, where t denotes the standard coordinate on $\mathbb{R}$ factor. In this case, $\mathcal{M}\times \mathbb{R}$ is said to be the symplectification of $\mathcal{M}$. The inverse of this assertion is also true. That is, if $(\mathcal{M}\times \mathbb{R},d\left(e^t\eta \right))$ is a symplectic manifold, then $(\mathcal{M},\eta )$ turns out to be contact.

3.2. Submanifolds of Contact Manifolds

Let $(\mathcal{M},\eta )$ be a contact manifold. Recall the associated bi-vector field $\mathsf{\Lambda }$ defined in (96). Consider a linear subbundle $\mathsf{\Xi }$ of the tangent bundle $T\mathcal{M}$ (that is, a distribution on $\mathcal{M}$). We define the contact complement of $\mathsf{\Xi }$ as

$${\mathsf{\Xi }}^{\perp }:={♯}_{\mathsf{\Lambda }}\left({\mathsf{\Xi }}^o\right),$$

(100)

where the sharp map on the right hand side is the one in (98) and ${\mathsf{\Xi }}^o$ is the annihilator of $\mathsf{\Xi }$. Let $\mathcal{N}$ be a submanifold of $\mathcal{M}$. We say that $\mathcal{N}$ is:

• Isotropic if $T\mathcal{N}\subseteq {T\mathcal{N}}^{\perp }$.
• Coisotropic if $T\mathcal{N}\supseteq {T\mathcal{N}}^{\perp }$.
• Legendrian if $T\mathcal{N}={T\mathcal{N}}^{\perp }$.

Assume that a submanifold $\mathcal{N}$ of a contact manifold $\mathcal{M}$ is defined to be the zero level set of k real smooth functions ${\varphi }_a:U\to \mathbb{R}$. We determine k vector fields $Z_a={♯}_{\mathsf{\Lambda }}\left(d{\varphi }_a\right)$. The image space of these vector fields are spanning the contact complement

$${T\mathcal{N}}^{\perp }=\mathrm{span}\left\{Z_a\right|a=1,\dots ,k\}.$$

(101)

In this geometry, $\mathcal{N}$ is coisotropic if and only if, $Z_a\left({\varphi }_b\right)=0$ for all $a,b$. According to the local computation in (99), we have that $\mathcal{N}$ is co-isotropic if and only if

$$\frac{\partial {\varphi }_a}{\partial p_i}\left(\frac{\partial {\varphi }_b}{\partial q^i}+p_i\frac{\partial {\varphi }_b}{\partial z}\right)-(\frac{\partial {\varphi }_a}{\partial q^i}+p_i\frac{\partial {\varphi }_a}{\partial z})\frac{\partial {\varphi }_b}{\partial p_i}=0.$$

(102)

Referring to this local observation, one can easily prove that a submanifold $\mathcal{N}$ of $\mathcal{M}$ is Legendrian if and only if it is a maximal integral manifold of $ker\eta$. In this case, the dimension of $\mathcal{N}$ must be n (see [4,37,69]).

Generating Functions. Consider the first order jet bundle ${\mathcal{T}}^{\ast }\mathcal{Q}$ endowed with the contact structure given in (91). Let F be a real valued function on the base manifold $\mathcal{Q}$. Its first prolongation is a section of the bundle ${\pi }_{\mathcal{Q}}^0$ displayed in (89) that is,

$${\mathcal{T}}^{\ast }F:\mathcal{Q}⟶{\mathcal{T}}^{\ast }\mathcal{Q}=T^{\ast }\mathcal{Q}\times \mathbb{R},q\mapsto (dF\left(q\right),F\left(q\right)).$$

(103)

The image space of the first prolongation ${\mathcal{T}}^{\ast }F$ is a Legendrian submanifold of ${\mathcal{T}}^{\ast }\mathcal{Q}$. The converse of this assertion is also true, that is, if the image space of a section $\sigma$ of ${\pi }_{\mathcal{Q}}^0$ is a Legendrian submanifold, then it is the first prolongation of a function F. Evidently, this is not the only way to obtain a Legendrian submanifold.

Consider, for example, a Morse family E defined on a smooth bundle $(\mathcal{W},\tau ,\mathcal{Q})$ according to Section 2.2. Then, referring to the definition (26), we define a Lagrangian submanifold $\mathcal{S}$ of the cotangent bundle $T^{\ast }\mathcal{Q}$. In the light of the first jet prolongation in (103), we lift this Lagrangian submanifold to a Legendrian submanifold of the contact manifold ${\mathcal{T}}^{\ast }\mathcal{Q}$. To see this, consider a local system of coordinates $\left(q^i\right)$ on the base manifold $\mathcal{Q}$, and the induced coordinates $(q^i,{ϵ}^a)$ on the total space $\mathcal{W}$. Then, referring to the Darboux’s coordinates on ${\mathcal{T}}^{\ast }\mathcal{Q}$, the Legendrian submanifold $\mathcal{N}$ generated by a Morse family $E=E(q,ϵ)$ is computed to be

$$\mathcal{N}= \left\{\left(q^i,\frac{\partial E}{\partial q^i}(q,ϵ),E(q,ϵ)\right) \in {\mathcal{T}}^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {ϵ}^a}=0\right\} \subset {\mathcal{T}}^{\ast }\mathcal{Q}.$$

(104)

On the other hand, the lift of a Legendrian submanifold to the symplectification is a Lagrangian submanifold. In fact, $\mathcal{N}$ is a Legendrian submanifold of a contact manifold $(\mathcal{M},\eta )$ if and only if $\mathcal{N}\times \mathbb{R}$ is a Lagrangian submanifold of the symplectic manifold $(\mathcal{M}\times \mathbb{R},d\left(e^t\eta \right))$, see [70].

3.3. Contact Diffeomorphisms and Contact Hamiltonian Systems

Let $({\mathcal{M}}_1,{\eta }_1)$ and $({\mathcal{M}}_2,{\eta }_2)$ be two contact manifolds. A diffeomorphism $\phi$ from ${\mathcal{M}}_1$ to ${\mathcal{M}}_2$ is said to be a contact diffeomorphism (or contactomorphism) if it preserves the contact structures that is, $T\phi (ker{\eta }_1)=ker{\eta }_2$. In terms of the contact forms, a contact diffeomorphism $\phi$ is the one satisfying

$${\phi }^{\ast }{\eta }_2=\mu {\eta }_1.$$

(105)

where $\mu$ is a non-zero conformal factor. To manifest the existence of this conformal factor, a mapping $\phi$ satisfying (105) is also called as conformal contact diffeomorphism. In this understanding, the contact mapping is denoted by a two-tuple $(\phi ,\mu )$.

For a contact manifold $(\mathcal{M},\eta )$, we denote the group of contact diffeomorphisms [71] by

$${\mathrm{Diff}}_{con}\left(\mathcal{M}\right)= \left\{\phi \in \mathrm{Diff}\left(\mathcal{M}\right):{\phi }^{\ast }\eta =\mu \eta ,\mu \in \mathcal{F}\left(\mathcal{M}\right)\right\}.$$

(106)

Here, $\mathrm{Diff}\left(\mathcal{M}\right)$ is standing for the group of all diffeomorphism on $\mathcal{M}$. A vector field on the contact manifold $\left(\mathcal{M},\eta \right)$ is a contact vector field (called also as infinitesimal conformal contactomorphism) if it generates one-parameter group of contact diffeomorphisms. We have that the space of contact vector fields is given by

$${\mathfrak{X}}_{con}\left(\mathcal{M}\right)= \left\{X\in \mathfrak{X}\left(\mathcal{M}\right):{\mathcal{L}}_X\eta =-\lambda \eta ,\lambda \in \mathcal{F}\left(\mathcal{M}\right)\right\}.$$

(107)

Sometimes a contact vector field is denoted by a two-tuple $(X,\lambda )$ in order to manifest the existence of the conformal factor $\lambda$. In order to discuss the geometry of $\lambda$, we perform the following observation. This permits us to introduce Hamiltonian dynamics on the present framework as well (for more details, see [35,36,37,72]).

(Contact) Hamiltonian Vector Fields. For a real valued function H on a contact manifold $(\mathcal{M},\eta )$, there is a corresponding contact vector field $X_H^c$, defined as follows:

$${\iota }_{X_H^c}\eta =-H,{\iota }_{X_H^c}d\eta =dH-\mathcal{R}\left(H\right)\eta ,$$

(108)

where $\mathcal{R}$ is the Reeb vector field. Here, H is called the (contact) Hamiltonian function and $X_H^c$ is called the (contact) Hamiltonian vector field. We denote a contact Hamiltonian system as a three-tuple $(\mathcal{M},\eta ,H)$ where $(\mathcal{M},\eta )$ is a contact manifold and H is a smooth real function on M. A direct computation determines the conformal factor for a given Hamiltonian vector fields as

$${\mathcal{L}}_{X_H^c}\eta =d{\iota }_{X_H^c}\eta + {\iota }_{X_H^c}d\eta =-\mathcal{R}\left(H\right)\eta .$$

(109)

That is, $\lambda =\mathcal{R}\left(H\right)$.

In this realization, the contact Jacobi bracket of two smooth functions on $\mathcal{M}$ is defined by

$${\{F,H\}}^c={\iota }_{[X_F^c,X_H^c]}\eta ,$$

(110)

where $X_F$ and $X_H$ are Hamiltonian vectors fields determined through (108). Here, $\left[•,•\right]$ is the Lie bracket of vector fields. Then, the identity

$$- \left[X_K^c,X_H^c\right] =X_{{\left\{K,H\right\}}^c}^c$$

(111)

establishes the isomorphism

$$\left({\mathfrak{X}}_{con}\left(\mathcal{M}\right),-\left[•,•\right]\right) ⟷ \left(\mathcal{F}\left(\mathcal{M}\right),{\left\{•,•\right\}}^c\right)$$

(112)

between the Lie algebras of real smooth functions and contact vector fields.

According to (109), the flow of a contact Hamiltonian system preserves the contact structure, but it does not preserve neither the contact one-form nor the Hamiltonian function. Instead, we obtain

$${\mathcal{L}}_{X_H^c}H=-\mathcal{R}\left(H\right)H.$$

(113)

Being a non-vanishing top-form we can consider $d{\eta }^n\wedge \eta$ as a volume form on $\mathcal{M}$. Hamiltonian motion does not preserve the volume form since

$${\mathcal{L}}_{X_H^c}(d{\eta }^n\wedge \eta )=-(n+1)\mathcal{R}\left(H\right)d{\eta }^n\wedge \eta .$$

(114)

However, it is immediate to see that, for a nowhere vanishing Hamiltonian function H, the quantity $H^{-(n+1)}{\left(d\eta \right)}^n\wedge \eta$ is preserved along the motion (see [41]).

Referring to the Darboux’s coordinates $(q^i,p_i,z)$, for a Hamiltonian function H, the Hamiltonian vector field, determined in (108), is computed to be

$$X_H^c=\frac{\partial H}{\partial p_i}\frac{\partial }{\partial q^i}- \left(\frac{\partial H}{\partial q^i}+\frac{\partial H}{\partial z}{p}_i\right)\frac{\partial }{\partial p_i}+(p_i\frac{\partial H}{\partial p_i}-H)\frac{\partial }{\partial z},$$

(115)

whereas the contact Jacobi bracket (110) is

$${\{F,H\}}^c=\frac{\partial F}{\partial q^i}\frac{\partial H}{\partial p_i}-\frac{\partial F}{\partial p_i}\frac{\partial H}{\partial q^i}+ \left(F-p_i\frac{\partial F}{\partial p_i}\right)\frac{\partial H}{\partial z}- \left(H-p_i\frac{\partial H}{\partial p_i}\right)\frac{\partial F}{\partial z}.$$

(116)

Thus, we obtain that the Hamilton’s equations for H as

$${\dot{q}}^i=\frac{\partial H}{\partial p_i},{\dot{p}}_i=-\frac{\partial H}{\partial q^i}-p_i\frac{\partial H}{\partial z},\dot{z}=p_i\frac{\partial H}{\partial p_i}-H.$$

(117)

Evolution vector fields Another vector field can be defined from a Hamiltonian function H on a contact manifold $(M,\equiv )$: the evolution vector field of H [52], denoted as ${\epsilon }_H$, which is the one that satisfies

$${\mathcal{L}}_{{\epsilon }_H}\eta =dH-\mathcal{R}\left(H\right)\eta ,\eta \left({\epsilon }_H\right)=0.$$

(118)

In local coordinates, it is given by

$${\epsilon }_H=\frac{\partial H}{\partial p_i}\frac{\partial }{\partial q^i}- \left(\frac{\partial H}{\partial q^i}+\frac{\partial H}{\partial z}{p}_i\right)\frac{\partial }{\partial p_i}+p_i\frac{\partial H}{\partial p_i}\frac{\partial }{\partial z},$$

(119)

so that the integral curves satisfy the evolution equations

$${\dot{q}}^i=\frac{\partial H}{\partial p_i},{\dot{p}}_i=-\frac{\partial H}{\partial q^i}-p_i\frac{\partial H}{\partial z},\dot{z}=p_i\frac{\partial H}{\partial p_i}.$$

(120)

The evolution and Hamiltonian vector fields are related by

$${\epsilon }_H=X_H^c+H\mathcal{R}.$$

(121)

Quantomorphisms. By asking the conformal factor $\mu$ in the definition (105) to be the unity, one arrives the conservation of the contact forms

$${\phi }^{\ast }{\eta }_2={\eta }_1.$$

(122)

We call such a mapping as a strict contact diffeomorphism (or quantomorphism). For a contact manifold $(\mathcal{M},\eta )$ we denote the space of all strict contact transformations as

$${\mathrm{Diff}}_{con}^{st}\left(\mathcal{M}\right)= \left\{\phi \in \mathrm{Diff}\left(\mathcal{M}\right):{\phi }^{\ast }\eta =\eta \right\}\subset {\mathrm{Diff}}_{con}\left(\mathcal{M}\right).$$

(123)

The Lie algebra of this group is consisting of the infinitesimal quantomorphisms

$${\mathfrak{X}}_{con}^{st}\left(\mathcal{M}\right)= \left\{X\in {\mathfrak{X}}_{con}\left(\mathcal{M}\right):{\mathcal{L}}_{X_H}\eta =0\right\}.$$

(124)

If the contact vector field is determined through a smooth function H as in (108), then $X_H$ falls into the subspace ${\mathfrak{X}}_{con}^{st}\left(\mathcal{M}\right)$ if and only if $\lambda =-dH\left(\mathcal{R}\right)=0$. This reads that, to generate an infinitesimal quantomorphism, a function H must not depend on the fiber variable z.

Now, consider the canonical contact manifold $({\mathcal{T}}^{\ast }\mathcal{Q},{\eta }_{\mathcal{Q}})$. For two functions, those that are not dependent on the fiber variable z, the contact Jacobi bracket ${\{•,•\}}^c$ in (116) locally turns out to be equal to the canonical Poisson bracket on $T^{\ast }Q$, therefore we find that

$$\left[X_H^{st},X_F^{st}\right] =-X_{\left\{H,F\right\}}^{st},$$

(125)

where $X_H^{st}$ is the infinitesimal quantomorphism generated by H. Accordingly, one arrives at an isomorphism

$${\mathfrak{X}}_{ham}\left(T^{\ast }\mathcal{Q}\right)⟶{\mathfrak{X}}_{con}^{st}\left({\mathcal{T}}^{\ast }\mathcal{Q}\right),X_H\mapsto X_H-H\frac{\partial }{\partial z}$$

(126)

from the Lie algebra of Hamiltonian vector fields ${\mathfrak{X}}_{ham}\left(T^{\ast }\mathcal{Q}\right)$ into the canonical cotangent bundle $T^{\ast }\mathcal{Q}$ to the Lie algebra of infinitesimal quantomorphisms ${\mathfrak{X}}_{con}^{st}\left({\mathcal{T}}^{\ast }\mathcal{Q}\right)$ on the extended cotangent bundle.

3.4. Contact Lagrangian Dynamics

Once more, we consider the extended configuration space $\mathcal{Q}\times \mathbb{R}$, but this time take it as the total space of the standard fibration from $\mathcal{Q}\times \mathbb{R}$ to $\mathbb{R}$. In this bundle structure, the base manifold is $\mathbb{R}$ and the fibration is simply the projection to the second factor. The first jet manifold is diffeomorphic to

$$\mathcal{T}\mathcal{Q}=T\mathcal{Q}\times \mathbb{R}.$$

(127)

We call this space the extended tangent bundle.

Suppose that $L:\mathcal{T}\mathcal{Q}\mapsto \mathbb{R}$ is a Lagrangian function. In order to arrive at the dynamical equations governed by such a Lagrangian function, one needs to employ the Herglotz principle, which is defined by an action functional [38,53,73].

The value of the functional attains its extremum if $q\left(t\right)$ is a solution of the Herglotz equations (also known as the generalized Euler–Lagrange equations):

$${\dot{q}}^i=\frac{d{q}^i}{dt},\frac{\partial L}{\partial q^i}-\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}^i}\right) + \frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}=0,$$

(128)

and z is a solution of the Cauchy problem

$$\dot{z}=L(t,q^i,{\dot{q}}^i,z),0\le t\le \tau .$$

(129)

It is important to notice that (128) represents a family of differential equations, since for each function $q\left(t\right)$, a different differential equation arises, hence $z\left(t\right)$ depends on $q\left(t\right)$. Without the explicit dependence of z, this problem reduces to a problem of the classical calculus of variations. If the functional z defined in (129) is invariant with respect to translation in time, then the quantity

$$I=exp\left(-{\int }^t\frac{\partial L}{\partial z}d\theta \right)\left(L(q,\dot{q},z)-\frac{\partial L}{\partial {\dot{q}}^i}{\dot{q}}^i\right)$$

(130)

is conserved on solutions of the Herglotz equations for regular Lagrangians.

For a regular Lagrangian function L, the fiber derivative determines a diffeomorphism from the extended tangent bundle $\mathcal{T}\mathcal{Q}$ to the extended cotangent bundle ${\mathcal{T}}^{\ast }\mathcal{Q}$ as

$$\mathbb{F}{L}^c:\mathcal{T}\mathcal{Q}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^i,z)\mapsto (q^i,\frac{\partial L}{\partial {\dot{q}}^j},z)$$

(131)

A direct calculation shows that the Legendre transformation (131) maps the Herglotz equations in (128) to the contact Hamilton’s equations (117) if the Hamiltonian function is defined to be

$$H(q^i,p_i,z)={\dot{q}}^i{p}_i-L(q,\dot{q},z).$$

(132)

The evolution-Herglotz equations. One can obtain the Lagrangian formalism for the evolution vector field by using a nonlinear nonholonomic action principle, as shown in [74]. The resulting equations are the evolution Herglotz equations

$$\begin{array}{c}\frac{\partial L}{\partial q^i}-\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}^i}\right) + \frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}=0,\\ \dot{z}={\dot{q}}^i\frac{\partial L}{\partial {\dot{q}}^i}.\end{array}$$

(133)

The Legendre transformation (131) maps the evolution Herglotz equations in (133) to the evolution contact Hamiltonian dynamics in (120).

4. Tulczyjew’s triple for Contact Geometry

4.1. Special Contact Structures

Let $(\mathcal{M},\eta )$ be a contact manifold and the total space of a fibre bundle $(\mathcal{M},\rho ,\mathcal{Q})$. We introduce a special contact structure as a quintuple

$$(\mathcal{M},\rho ,\mathcal{Q},\eta ,\mathsf{\Phi }),$$

(134)

where $\mathsf{\Phi }$ is a fiber preserving contact diffeomorphism from $\mathcal{M}$ to the canonical contact manifold $({\mathcal{T}}^{\ast }\mathcal{Q},{\eta }_{\mathcal{Q}})$. The two-tuple $(\mathcal{M},\eta )$ is said to be the underlying contact manifold of the special contact structure. Here, we have a diagram exhibiting a special contact structure in a pictorial way

(135)

where ${\pi }_{\mathcal{Q}}^0$ is the fibration given in (89).

It is possible to define a special contact space starting with a special symplectic space defined in Section 2.1. For this, suppose that $(\mathcal{P},\omega =-d\theta )$ is an exact symplectic manifold and in the product manifold $\mathcal{P}\times \mathbb{R}$ we consider the standard contact structure $\eta =dz-\theta$ (that is, the contactization of the exact symplectic structure $\omega =-d\theta$). Assume also that $\mathcal{P}$ admits a special symplectic structure $(\mathcal{P},\pi ,\mathcal{Q},\theta ,\varphi )$ as pictured in (22). Then, $\mathcal{M}$ admits a special contact structure $(\mathcal{M},\rho ,\mathcal{Q},\eta ,\mathsf{\Phi })$, where $\mathsf{\Phi }(p,z)=\left(\varphi \right(p),z)$, and the following diagram

(136)

is commutative, with ${\pi }_{\mathcal{Q}}^1$ the fibration given in (89). Here, it is considered that $\rho =\pi \circ \mathit{pr}$. This construction is the contactization of the special symplectic structure.

We now merge a Morse family E defined on a fiber bundle $(\mathcal{W},\tau ,\mathcal{Q})$ and a special contact space $(\mathcal{M},\rho ,\mathcal{Q},\eta ,\mathsf{\Phi })$ in order to arrive at a Legendrian submanifold of $(\mathcal{M},\eta )$. For this, consider the following commutative diagram

(137)

Referring to the definition in (104), we obtain a Legendrian submanifold $\mathcal{N}$ of the jet bundle ${\mathcal{T}}^{\ast }\mathcal{Q}$. Then, by employing the inverse of the contact diffeomorphism $\mathsf{\Phi }$, we arrive at a Legendrian submanifold ${\mathcal{N}}_E$ of $\mathcal{M}$. Referring to this realization, we shall exhibit both the contact Hamiltonian and contact Lagrangian dynamics as Legendrian submanifolds of the same contact manifold in the following subsection.

4.2. Tangent Contact Manifold

We start by lifting a contact structure $\eta$ on a contact manifold $\mathcal{M}$ a contact structure on the extended tangent bundle $\mathcal{T}\mathcal{M}$. This lifting is in introduced in [70] to characterize the contact vector fields on $\mathcal{M}$ (in particular, the Hamiltonian vector fields in $\mathcal{M}$) in terms of Legendrian submanifolds of the contact manifold $\mathcal{T}\mathcal{M}$. In fact, more later, in this direction, we shall use some other results those available in [70].

Theorem 4.

For a contact manifold $(\mathcal{M},\eta )$, the extended tangen bundle $\mathcal{T}\mathcal{M}\simeq T\mathcal{M}\times \mathbb{R}$ is a contact manifold by admitting a contact one-form

$${\eta }^{\mathcal{T}}:=u{\eta }^V+{\eta }^C$$

(138)

where u is coordinate on $\mathbb{R}$ whereas ${\eta }^C$ and ${\eta }^V$ are the complete and vertical lifts of η, respectively.

The one-form ${\eta }^{\mathcal{T}}$ is said to be the tangent contact structure and we will denote the tangent contact manifold as a two-tuple

$$(\mathcal{T}\mathcal{M},{\eta }^{\mathcal{T}})=(T\mathcal{M}\times \mathbb{R},u{\eta }^V+{\eta }^C).$$

(139)

Contact Hamiltonian Dynamics as a Legendrian Submanifold. Let $(\mathcal{M},\eta )$ be a contact manifold. Consider a vector field X, a real valued function $\lambda$ on $\mathcal{M}$, hence a section

$$(X,\lambda ):\mathcal{M}⟶\mathcal{T}\mathcal{M}=T\mathcal{M}\times \mathbb{R},m\mapsto \left(X\right(m),\lambda (m\left)\right).$$

(140)

of the fibration ${\tau }_{\mathcal{M}}^0:\mathcal{T}\mathcal{M}\mapsto \mathcal{M}$. We plot the following commutative diagram to see this

(141)

Using Theorem 3.13 in [70] and the comments at the beginning of this subsection, we deduce that the pair $(X,\lambda )$ is a contact vector field (an infinitesimal conformal contactomorphism), that is, an element of ${\mathfrak{X}}_{con}\left(\mathcal{M}\right)$ in (107) if and only if the image space of $(X,\lambda )$ is a Legendrian submanifold of the tangent contact manifold $(\mathcal{T}\mathcal{M},{\eta }^{\mathcal{T}})$. This result states that the image of a contact Hamiltonian vector field $X_H^c$, after suitably included in the contactified tangent bundle, turns out to be a Legendrian submanifold. As discussed in the previous section, the conformal factor $\lambda$ in the present case is $\mathcal{R}\left(H\right)$. So that, the image of the mapping

$$(X_H^c,\mathcal{R}\left(H\right)):\mathcal{M}⟶\mathcal{T}\mathcal{M}=T\mathcal{M}\times \mathbb{R},m\mapsto (X_H^c\left(m\right),\mathcal{R}\left(H\right)\left(m\right))$$

(142)

is a Legendrian submanifold of the tangent contact manifold $\mathcal{T}\mathcal{M}$. To see this more clearly, let us discuss this geometry in the realm of special contact spaces.

Consider a contact manifold $(\mathcal{M},\eta )$. Its extended tangent bundle $\mathcal{T}\mathcal{M}$ is a contact manifold endowed with the contact structure given by (139) and ${\mathcal{T}}^{\ast }\mathcal{M}$, as an extended cotangent bundle, is also a contact manifold. Between these two extended spaces, we introduce a fiber preserving contact diffeomorphism ${\beta }^c$ as follows

$${\beta }^c:\mathcal{T}\mathcal{M}⟶{\mathcal{T}}^{\ast }\mathcal{M},(V,u)\mapsto \left(-{\iota }_{V}d\eta -u\eta ,\eta \left(V\right)\right)$$

(143)

where $\eta$ is the contact one-form on $\mathcal{M}$. It is a direct computation to see that

$${\left({\beta }^c\right)}^{\ast }{\eta }_{{\mathcal{T}}^{\ast }\mathcal{M}}={\eta }^{\mathcal{T}}$$

(144)

where ${\eta }_{{\mathcal{T}}^{\ast }\mathcal{M}}$ is the contact one-form on the canonical contact manifold ${\mathcal{T}}^{\ast }\mathcal{M}$ and ${\eta }^{\mathcal{T}}$ is the lifted contact one-form on $\mathcal{T}\mathcal{M}$ given in (138). This observation permits us to determine a special contact space and following the order given in (134), we write this special contact manifold as

$$(\mathcal{T}\mathcal{M},{\tau }_{\mathcal{M}}^0,\mathcal{M},{\eta }^{\mathcal{T}},{\beta }^c).$$

(145)

Accordingly, following the picture in (135), we plot the following diagram

(146)

where we have employed the projections ${\pi }_{\mathcal{M}}^0:{\mathcal{T}}^{\ast }\mathcal{M}\mapsto \mathcal{M}$ and ${\tau }_{\mathcal{M}}^0:\mathcal{T}\mathcal{M}\mapsto \mathcal{M}$. This diagram also manifests how one can transfer the Legendrian submanifolds one onto the other. If H is a Hamiltonian function on the contact manifold $\mathcal{M}$, then the image space of $-{\mathcal{T}}^{\ast }H$, as defined in (103), is a Legendrian submanifold of ${\mathcal{T}}^{\ast }\mathcal{M}$. Thus, since ${\beta }^c$ is a contact diffeomorphism, by pulling back the image space of $-{\mathcal{T}}^{\ast }H$, we arrive at a Legendrian submanifold of the tangent contact manifold $\mathcal{T}\mathcal{M}$. However, using (143), we have that

$${\beta }^c\circ (X_H^c,\mathcal{R}\left(H\right))=-{\mathcal{T}}^{\ast }H=-(dH,H).$$

(147)

This is the contact version of the identity (82). Notice that, this observation can be considered as an indirect proof of the assertion that the image space of a contact vector field is a Legendrian submanifold.

Local Picture. Consider Darboux’s coordinates $(q^i,p_i,z)$ on $\mathcal{M}$ then we assume the induced coordinates on $\mathcal{T}\mathcal{M}\simeq T\mathcal{M}\times \mathbb{R}$ as $(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,\dot{z},u)$. In this local realization, the lifted contact one-form ${\eta }^{\mathcal{T}}$ defined in (138) is computed to be

$$\begin{array}{cc}\hfill {\eta }^{\mathcal{T}}& =u{\eta }^V+{\eta }^C=u(dz-p_{i}d{q}^i)+(d\dot{z}-{\dot{p}}_{i}d{q}^i-p_{i}d{\dot{q}}^i)\hfill \\ \hfill & =d\dot{z}+udz-({\dot{p}}_i+u{p}_i)d{q}^i-p_{i}d{\dot{q}}^i,\hfill \end{array}$$

(148)

and the Reeb vector field is ${\mathcal{R}}^{\mathcal{T}}=\partial /\partial \dot{z}$. In this realization, the contact mapping ${\beta }^c$ in (143) turns out to be

$${\beta }^c(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,\dot{z},u)=(q^i,p_i,z,u{p}_i+{\dot{p}}_i,-{\dot{q}}^i,-u,\dot{z}-p_i{\dot{q}}^i).$$

(149)

The fact that ${\beta }^c$ is a contactomorphism, that is the identity (144), follows from a direct calculation in coordinates. Observe that,

$$\begin{array}{cc}\hfill {\left({\beta }^c\right)}^{\ast }\left({\eta }_{{\mathcal{T}}^{\ast }\mathcal{M}}\right)& =d(\dot{z}-p_i{\dot{q}}^i)-(u{p}_i+{\dot{p}}_i)q^i+{\dot{q}}^{i}d{p}_i+udz\hfill \\ & =u(dz-p_{i}d{q}^i)+(d\dot{z}-{\dot{p}}_{i}d{q}^i-p_{i}d{\dot{q}}^i)={\eta }^{\mathcal{T}}.\hfill \end{array}$$

(150)

Consider now a Hamiltonian function H on the contact manifold $\mathcal{M}$. Minus of its first prolongation defines a Legendrian submanifold of ${\mathcal{T}}^{\ast }\mathcal{M}$ given by

$$\mathit{im}(-{\mathcal{T}}^{\ast }H)= \left\{(q^i,p_i,z,-\frac{\partial H}{\partial q^i},-\frac{\partial H}{\partial p_i},-\frac{\partial H}{\partial z},-H)\in {\mathcal{T}}^{\ast }\mathcal{M}:H=H(q,p,z)\right\}.$$

(151)

Referring to the inverse of the contact diffeomorphism ${\beta }^c$ we compute a Legendrian submanifold of $\mathcal{T}\mathcal{M}$ as

$$\begin{array}{cc}\hfill {\mathcal{N}}_{-H}& ={\left({\beta }^c\right)}^{-1}\left(\mathit{im}(-{\mathcal{T}}^{\ast }H)\right)\hfill \\ \hfill & = \left\{(q^i,p_i,z,\frac{\partial H}{\partial p_i},-\frac{\partial H}{\partial z}{p}_i-\frac{\partial H}{\partial q^i},p_i\frac{\partial H}{\partial p_i}-H,\frac{\partial H}{\partial z})\in \mathcal{T}\mathcal{M}:H=H(q,p,z)\right\}.\hfill \end{array}$$

(152)

It is immediate to see that the Legendrian submanifold can be alternatively obtained by

$${\mathcal{N}}_{-H}=im(X_H^c,\mathcal{R}\left(H\right)).$$

(153)

So, the contact Hamiltonian dynamics are determined by the Legendrian submanifold ${\mathcal{N}}_{-H}$ as follows. Let $\sigma :I\mapsto \mathcal{M}$ be a smooth curve on $\mathcal{M}$ and consider the lift of $\sigma$ to $\mathcal{T}\mathcal{M}$ defined by

$${\sigma }^{\mathcal{T}}= (\dot{\sigma }\left(t\right),\mathcal{R}\left(H\right)\left(\sigma \left(t\right)\right),\forall t\in I,$$

(154)

where $\dot{\sigma }$ is the tangent lift to $T\mathcal{M}$. Then, $\sigma$ is a solution of the contact Hamilton’s equations for H if and only if its lift ${\sigma }^{\mathcal{T}}$ to $\mathcal{T}\mathcal{M}$ is contained in the Legendrian submanifold ${\mathcal{N}}_{-H}$. In fact, if the local expression of $\sigma$ is

$$\sigma \left(t\right)=(q^i\left(t\right),p_i\left(t\right),z\left(t\right))$$

(155)

then, from (152), it follows that ${\sigma }^{\mathcal{T}}\in {\mathcal{N}}_{-H}$ for every $t\in I$, if and only if $\sigma$ satisfies the following equations

$$\frac{d{q}^i}{dt}=\frac{\partial H}{\partial p_i},\frac{d{p}_i}{dt}=-\frac{\partial H}{\partial z}{p}_i-\frac{\partial H}{\partial q^i},\frac{dz}{dt}=p_i\frac{\partial H}{\partial p_i}-H.$$

(156)

They are precisely the contact Hamilton’s equations in (117).

4.3. Contact Tulczyjew’s Triple

Classical Tulczyjew’s triple is obtained by properly merging two special symplectic structures. Following the same understanding, we introduce Tulczyjew’s triple for contact dynamics by properly merging two special contact structures.

Contact Lift of Diffeomorphism $\psi$ in (56). Recall the symplectic diffeomorphism $\psi$, defined in (56), from the cotangent bundle $T^{\ast }T\mathcal{Q}$ to the cotangent bundle $T^{\ast }{T}^{\ast }\mathcal{Q}$. We extend this mapping to the level of contact manifolds

$${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}\simeq T^{\ast }(T^{\ast }\mathcal{Q}\times \mathbb{R})\times \mathbb{R},{\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}\simeq T^{\ast }(T\mathcal{Q}\times \mathbb{R})\times \mathbb{R}$$

(157)

by assuming that extension is the identity on $T\mathbb{R}$ and it vanishes on the zero section. Accordingly, we recall the mapping $\psi$ in (56), consider a local chart $(q^i,{\dot{q}}^i,z,a_i,{\dot{a}}_i,v,u)$ on ${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}$ then we have an arbitrary extension

$${\psi }^c:{\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}⟶{\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^i,z,a_i,{\dot{a}}_i,v,u)\mapsto (q^i,{\dot{a}}_i,z,a_i,-{\dot{q}}^i,v,w).$$

To determine w, we compute

$${\left({\psi }^c\right)}^{\ast }{\eta }_{{\mathcal{T}}^{\ast }\mathcal{Q}}=dw-{\psi }^{\ast }{\theta }_{T^{\ast }\mathcal{Q}}-vs.dz=dw-du+d\left({\dot{q}}^i{\dot{a}}_i\right)+{\eta }_{\mathcal{T}\mathcal{Q}},$$

(158)

where ${\eta }_{{\mathcal{T}}^{\ast }\mathcal{Q}}$ and ${\eta }_{\mathcal{T}\mathcal{Q}}$ are the canonical contact one-forms on ${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}$ and ${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}$, respectively. Hence, ${\psi }^c$ is a contact mapping if and only if $dw=d(u-{\dot{q}}^i{\dot{a}}_i)$. Since we ask w to vanish on the zero section, necessarily $w=u-{\dot{q}}^i{\dot{a}}_i$. Therefore, the unique extension is, locally, computed to be

$${\psi }^c:{\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}⟶{\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^i,z,a_i,{\dot{a}}_i,v,u)\mapsto (q^i,{\dot{a}}_i,z,a_i,-{\dot{q}}^i,v,u-{\dot{a}}_i{\dot{q}}^i).$$

(159)

Next, we will present an intrinsic definition of ${\psi }^c$. For this purpose, we will use the following identifications

$${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}\cong (T^{\ast }T\mathcal{Q}\times T^{\ast }\mathbb{R})\times \mathbb{R},{\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}\cong (T^{\ast }{T}^{\ast }\mathcal{Q}\times T^{\ast }\mathbb{R})\times \mathbb{R}.$$

(160)

Then, we have that

$${\psi }^c((\gamma ,vdz),u)= \left(\psi \left(\gamma \right),vdz,u-\langle ({\pi }_{T^{\ast }\mathcal{Q}}\left(\psi \left(\gamma \right)\right),{\pi }_{T\mathcal{Q}}\left(\gamma \right)\rangle \right)$$

(161)

where $\langle •,•\rangle$ is the canonical pairing between $T^{\ast }\mathcal{Q}$ and $T\mathcal{Q}$.

The Left Wing of the triple. Now, using the contact mapping ${\psi }^c$ in (159), we define a contact diffeomorphism

$$\begin{array}{cc}\hfill {\alpha }^c:\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}⟶{\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q},V& \mapsto {\left({\psi }^c\right)}^{-1}\circ {\beta }^c\left(V\right)\hfill \\ \hfill (q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,\dot{z},u)& \mapsto (q^i,{\dot{q}}^i,z,u{p}_i+{\dot{p}}_i,p_i,-u,\dot{z}).\hfill \end{array}$$

(162)

Indeed, it is immediate to show that

$${\left({\alpha }^c\right)}^{\ast }{\eta }_{\mathcal{T}\mathcal{Q}}={\eta }_{\mathcal{Q}}^{\mathcal{T}},$$

(163)

where ${\eta }_{\mathcal{Q}}^{\mathcal{T}}$ is the lifted contact one-form on $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ given in (148). Then, we define the following special contact structure

$$(\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q},\mathcal{T}{\pi }_{\mathcal{Q}}^0,\mathcal{T}\mathcal{Q},{\eta }_{\mathcal{Q}}^{\mathcal{T}},{\alpha }^c)$$

(164)

which is diagrammatically given by

(165)

with the projections

$$\begin{array}{cc}\hfill {\pi }_{\mathcal{T}\mathcal{Q}}^0& :{\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}\simeq T^{\ast }\mathcal{T}\mathcal{Q}\times \mathbb{R}⟶\mathcal{T}\mathcal{Q},(W,u)\mapsto {\pi }_{\mathcal{T}\mathcal{Q}}\left(W\right)\hfill \\ \hfill \mathcal{T}{\pi }_{\mathcal{Q}}^0& :\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}\simeq T{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶\mathcal{T}\mathcal{Q},(V,s)\mapsto (T{\pi }_{\mathcal{Q}}\left(U\right),z).\hfill \end{array}$$

(166)

Here, we have employed the following global trivialization $V=(U,z,\dot{z})$ in $T{T}^{\ast }\mathcal{Q}\times T\mathbb{R}$.

Generalized Euler–Lagrange Equations as a Lagrangian Submanifold. Consider a Lagrangian function L on $\mathcal{T}\mathcal{Q}=T\mathcal{Q}\times \mathbb{R}$. The image space of its first prolongation, that is, $\mathit{im}\left({\mathcal{T}}^{\ast }L\right)$, is a Legendrian submanifold of ${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}$ computed to be

$$\mathit{im}\left({\mathcal{T}}^{\ast }L\right)= \left\{(q^i,{\dot{q}}^i,z,\frac{\partial L}{\partial q^i},\frac{\partial L}{\partial {\dot{q}}^i},\frac{\partial L}{\partial z},L)\in {\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}:L=L(q,\dot{q},z)\right\} \subset {\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}$$

(167)

Referring to the left wing of the contact triple (176), that is by applying the inverse of the mapping ${\alpha }^c$, we arrive at a Legendrian submanifold of $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ as

$$\begin{array}{cc}\hfill {\mathcal{N}}_L& ={\left({\alpha }^c\right)}^{-1}\left(\mathit{im}\left({\mathcal{T}}^{\ast }L\right)\right)\hfill \\ \hfill & = \left\{(q^i,\frac{\partial L}{\partial {\dot{q}}^i},z,{\dot{q}}^i,\frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}+\frac{\partial L}{\partial q^i},L,-\frac{\partial L}{\partial z})\in \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}:L=L(q,\dot{q},z)\right\} \subset \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}.\hfill \end{array}$$

(168)

The Lagrangian dynamics are determined by the Legendrian submanifold ${\mathcal{N}}_L$ as follows. First of all, we consider the transformation $\mathbb{F}L:\mathcal{T}\mathcal{Q}\mapsto {\mathcal{T}}^{\ast }\mathcal{Q}$ induced by L given by

$$\mathbb{F}L(u,z)=(\mathbb{F}{L}_1(u,z),z),$$

(169)

with $\mathbb{F}{L}_1(u,z)$ defined by

$$\langle \mathbb{F}{L}_1(u,z),u^{\prime }\rangle =\frac{d}{dt}{|}_{t=0}L(u+t{u}^{\prime },z).$$

(170)

Then, a smooth curve $c:I\mapsto \mathcal{T}\mathcal{Q}$ is a solution of the dynamics if and only if the lift ${(\mathbb{F}L\circ c)}^{\mathcal{T}}:I\mapsto \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ of the curve $\mathbb{F}L\circ c:I\mapsto {\mathcal{T}}^{\ast }\mathcal{Q}$ given by

$${(\mathbb{F}L\circ c)}^{\mathcal{T}}\left(t\right)= \left(\dot{\mathbb{F}L\circ c}\left(t\right),\frac{\partial L}{\partial z}\left(c\left(t\right)\right)\right),\forall t\in I$$

(171)

is contained in the Legendrian submanifold ${\mathcal{N}}_L$. In fact, if the local expression of c is

$$c\left(t\right)=(q^i\left(t\right),{\dot{q}}^i\left(t\right),z\left(t\right))$$

(172)

then, using (168), it follows that ${(\mathbb{F}L\circ c)}^{\mathcal{T}}\left(t\right)\in {\mathcal{N}}_L$, $\forall t\in I$, if and only if c satisfies

$$\frac{d{q}^i}{dt}={\dot{q}}^i,\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{q}}^i}\right) =\frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}+\frac{\partial L}{\partial q^i},\frac{dz}{dt}=L(q,\dot{q},z)$$

(173)

where the first two equations are the Herglotz equations in (128) and the latter is the Herglotz’s differential principle in (129). We remark that this theory does not require any regularity conditions on the Lagrangian function.

The Right Wing of the triple. To have the right wing, we simply replace the arbitrary contact manifold $(\mathcal{M},\eta )$ with the canonical contact manifold $({\mathcal{T}}^{\ast }\mathcal{Q},{\eta }_{\mathcal{Q}})$ in the special contact manifold (145). In this case, the total spaces become ${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}$ and $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ with the fibrations

$$\begin{array}{cc}\hfill {\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}^0& :\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}\simeq T{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(V,u)\mapsto {\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}\left(V\right),\hfill \\ \hfill {\pi }_{{\mathcal{T}}^{\ast }\mathcal{Q}}^0& :{\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}\simeq T^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(Z,u)\mapsto {\pi }_{{\mathcal{T}}^{\ast }\mathcal{Q}}\left(Z\right).\hfill \end{array}$$

(174)

We denote this special contact manifold as

$$(\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q},{\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}^0,{\mathcal{T}}^{\ast }\mathcal{Q},{\eta }_{\mathcal{Q}}^{\mathcal{T}},{\beta }^c).$$

(175)

Since Darboux’s coordinates are employed in Section 4.2, all the local formulations available in that subsection hold also in the present case.

Contact triple. We now merge the special contact spaces (164) and (175) in order to construct a Tulzyjew’s triple for the contact geometry. Accordingly, we couple the commutative diagrams in (146) and (165) and arrive at the following diagram

(176)

where the contact diffeomorphisms ${\beta }^c$ and ${\alpha }^c$ are those defined in (143) and (162), respectively. Here, the projections ${\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}^0$ and ${\pi }_{{\mathcal{T}}^{\ast }\mathcal{Q}}^0$ (respectively ${\pi }_{\mathcal{T}\mathcal{Q}}^0$ and $\mathcal{T}{\pi }_{\mathcal{Q}}^0$) are given by (174) (respectively, (166)).

4.4. Evolution Contact Tulczyjew’s Triple

We will analyze how the Tulczyjew’s triple can be used to understand the evolution dynamics. First of all, we notice that there is not, a natural way, to describe an evolution vector field on a contact manifold $\mathcal{M}$ as a Legendrian submanifold of $\mathcal{T}\mathcal{M}$. However, we will see that such a vector field may be considered as a Lagrangian submanifold of an exact symplectic submanifold of $\mathcal{T}\mathcal{M}$.

Consider a contact manifold $(\mathcal{M},\eta )$. In (86), we have realized the kernel of the contact one-form as a symplectic vector bundle $H\left(\mathcal{M}\right)$ of the tangent bundle $T\mathcal{M}$. Consider now the inclusion mapping

$$j=(i,{\mathrm{Id}}_{\mathbb{R}}):H\mathcal{M}\times \mathbb{R}\hookrightarrow \mathcal{T}\mathcal{M}=T\mathcal{M}\times \mathbb{R}$$

(177)

where i is the inclusion of $H\mathcal{M}$ into $T\mathcal{M}$. Consider now the one-form ${\theta }_{\eta }=j^{\ast }{\eta }^{\mathcal{T}}$ on $H\mathcal{M}\times \mathbb{R}$, and the two-form ${\omega }_{\eta }=d{\theta }_{\eta }$. In the following proposition, we state that ${\omega }_{\eta }$ is a symplectic two-form on $H\left(\mathcal{M}\right)\times \mathbb{R}$

Proposition 1.

The exact two-form ${\omega }_{\eta }=d{\theta }_{\eta }$ induces a symplectic structure on $H\mathcal{M}\times \mathbb{R}$.

Proof. 

Let $(q^i,p_i,z)$ be the Darboux’s coordinates on $\mathcal{M}$, so that $(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,\dot{z},u)$ are the induced coordinates for $\mathcal{T}\mathcal{M}$. In the light of the inclusion (177), we can employ $(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,u)$ as a local coordinate chart on $H\mathcal{M}\times \mathbb{R}$ since

$$j(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,u)=(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,p_i{\dot{q}}^i,u).$$

(178)

In other words,

$$H\mathcal{M}\times \mathbb{R}= \left\{(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,\dot{z},u)\in \mathcal{T}\mathcal{M}\times \mathbb{R}:\dot{z}-p_i{\dot{q}}^i=0\right\}.$$

(179)

Thus, the Reeb vector field ${\mathcal{R}}^{\mathcal{T}}=\partial /\partial \dot{z}$ of the contact manifold $(\mathcal{T}\mathcal{M},{\eta }^{\mathcal{T}})$ is transverse to the manifold of codimension one $H\mathcal{M}\times \mathbb{R}$. This implies that the exact two-form ${\omega }_{\eta }=d{\theta }_{\eta }=d\left(j^{\ast }{\eta }^{\mathcal{T}}\right)$ induces a symplectic structure on $H\mathcal{M}\times \mathbb{R}$.

Note that, from (148) and (178), it follows that

$$\begin{array}{cc}\hfill {\theta }_{\eta }& =udz-({\dot{p}}_i+u{p}_i)d{q}^i+{\dot{q}}^{i}d{p}_i,\hfill \\ \hfill {\omega }_{\eta }& =du\wedge dz-d{\dot{p}}_i\wedge d{q}^i-p_{i}du\wedge d{q}^i-ud{p}_i\wedge d{q}^i+d{\dot{q}}^i\wedge d{p}_i.\hfill \end{array}$$

(180)

□

Lagrangian submanifolds and evolution vector fields. Given a vector field X on $\mathcal{M}$ and a real smooth function $f:\mathcal{M}\to \mathbb{R}$, one can construct a section $(X,f):\mathcal{M}\to \mathcal{T}\mathcal{M}$ of the extended tangent bundle $\mathcal{T}\mathcal{M}$.

Theorem 5.

Let $(\mathcal{M},\eta )$ be a contact manifold. The map $(X,f):\mathcal{M}\mapsto \mathcal{T}\mathcal{M}$ defines a Lagrangian submanifold of the exact symplectic manifold $(H\mathcal{M}\times \mathbb{R},{\omega }_{\eta })$ if and only if $\eta \left(X\right)=0$ and ${\mathcal{L}}_X\eta +f\eta$ is closed.

Proof. 

First of all, note that the image of X lies on $H\mathcal{M}\times \mathbb{R}$ if and only if $\eta \left(X\right)=0$. For the second condition, we use well-known properties of complete and vertical lifts (see [62]);

$${(X,f)}^{\ast }{\eta }^{\mathcal{T}}={(X,f)}^{\ast }(u{\eta }^V+{\eta }^C)=f{X}^{\ast }{\eta }^V+X^{\ast }{\eta }^C=f\eta +{\mathcal{L}}_X\eta .$$

(181)

Thus, the image of $(X,f)$ is Lagrangian if and only if

$${(X,f)}^{\ast }{\omega }_{\eta }=d{(X,f)}^{\ast }{\eta }^{\mathcal{T}}=d(f\eta +{\mathcal{L}}_X\left(\eta \right))=0.$$

(182)

□

Looking at the definition of the evolution vector field (118), we obtain the following result.

Corollary 1.

The map $(X,f):\mathcal{M}\mapsto \mathcal{T}\mathcal{M}$ defines a Lagrangian submanifold of the exact symplectic manifold $(H\mathcal{M}\times \mathbb{R},{\omega }_{\eta })$ if and only if, locally, $X={\epsilon }_H$ and $f=\mathcal{R}\left(H\right)$ for a (local) smooth function $H:\mathcal{M}\to \mathcal{R}$.

Proof. 

Since the evolution vector field satisfies

$${\mathcal{L}}_{{\epsilon }_H}\eta =dH-\mathcal{R}\left(H\right)\eta ,\eta \left({\epsilon }_H\right)=0,$$

(183)

one arrives that ${\mathcal{L}}_X\eta +f\eta =dH$ is closed. Conversely, if $\eta \left(X\right)=0$ and ${\mathcal{L}}_X\eta +f\eta$ is closed, in a local chart U, one has that

$${\left.{\mathcal{L}}_X\eta +f\eta \right|}_U={\left.{\iota }_{X}d\eta +f\eta \right|}_U=dH$$

(184)

for a local smooth function H. By contracting the equality with the Reeb vector field, we obtain that $f=\mathcal{R}\left(H\right)$. Thus, $X={\epsilon }_H$. □

We start by presenting the following commutative diagram:

(185)

Here, ${\beta }^c$ is the contact mapping in (143), the mapping $({\mathrm{Id}}_{T\mathcal{M}},0)$ is the canonical inclusion of the cotangent bundle into the extended cotangent bundle as a zero section and the smooth map ${\beta }^0$ is given by

$${\beta }^0:H\mathcal{M}\times \mathbb{R}⟶T^{\ast }\mathcal{M},(V,u)\mapsto -{\iota }_{V}d\eta -u\eta .$$

(186)

In coordinates,

$${\beta }^0(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,u)=(q^i,p_i,z,u{p}_i+{\dot{p}}_i,-{\dot{q}}^i,-u).$$

(187)

Now, we are ready to plot the right wing of Tulczyjew’s triple for the evolution contact dynamics. For this, we replace the contact manifold $\mathcal{M}$ with the extended cotangent bundle ${\mathcal{T}}^{\ast }\mathcal{Q}$. Note that, in this case, we have

(188)

where ${\widehat{\tau }}_{{\mathcal{T}}^{\ast }\mathcal{Q}}$ is the projection taking a two-tuple $(V,u)$ mapping to ${\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}\left(V\right)$, that is,

$${\widehat{\tau }}_{{\mathcal{T}}^{\ast }\mathcal{Q}}:H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(V,u)\mapsto {\tau }_{{\mathcal{T}}^{\ast }\mathcal{Q}}\left(V\right)$$

(189)

using that $H{\mathcal{T}}^{\ast }\mathcal{Q}$ is a vector subbundle of $T{\mathcal{T}}^{\ast }\mathcal{Q}$. In order to see that the triangle commutes, we compute that

$${\beta }^0\circ ({\epsilon }_H,\mathcal{R}\left(H\right))=-{\iota }_{{\epsilon }_H}d{\eta }_{\mathcal{Q}}-\mathcal{R}\left(H\right){\eta }_{\mathcal{Q}}=-{\mathcal{L}}_{{\epsilon }_H}{\eta }_{\mathcal{Q}}-\mathcal{R}\left(H\right){\eta }_{\mathcal{Q}}=-dH,$$

(190)

where we have used (118) and the Cartan’s formula. Following notation presented in (18), we write that the quintuple

$$(H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R},{\widehat{\tau }}_{{\mathcal{T}}^{\ast }\mathcal{Q}},{\mathcal{T}}^{\ast }\mathcal{Q},-{\theta }_{{\eta }_{\mathcal{Q}}},{\beta }^0)$$

(191)

determines a special symplectic structure. Since ${\beta }^0$ is a symplectic diffeomorphism, we can realize one more time that the image space $im({\epsilon }_H,\mathcal{R}\left(H\right))$ as a Lagrangian submanifold of $H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}$. Moreover, if ${\sigma }^e:I\mapsto {\mathcal{T}}^{\ast }\mathcal{Q}$ is a smooth curve on ${\mathcal{T}}^{\ast }\mathcal{Q}$ and we consider the lift ${\left({\sigma }^e\right)}^{\mathcal{T}}:I\mapsto \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ of the curve ${\sigma }^e$ given by (154), then one may prove that ${\sigma }^e$ is a solution of the contact evolution equations for H if and only if ${\left({\sigma }^e\right)}^{\mathcal{T}}\left(t\right)\in im({\epsilon }_H,\mathcal{R}\left(H\right))$, for every $t\in I$.

Recall the mapping ${\alpha }^c$ given in (162). Consider the following commutative diagram

(192)

where j is the inclusion mapping in (178), $pr$ is the natural projection from the extended cotangent bundle ${\mathcal{T}}^{\ast }\mathcal{T}\mathcal{Q}=T^{\ast }\mathcal{T}\mathcal{Q}\times \mathbb{R}$ to the first factor that is the cotangent bundle $T^{\ast }\mathcal{T}\mathcal{Q}$. Using (162) and (178), we deduce that the local expression of ${\alpha }^0$ is

$${\alpha }^0:H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶T^{\ast }\mathcal{T}\mathcal{Q},(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,u)\mapsto (q^i,{\dot{q}}^i,z,u{p}_i+{\dot{p}}_i,p_i,-u).$$

(193)

A direct calculation proves that ${\left({\alpha }^0\right)}^{\ast }{\omega }_{\mathcal{T}\mathcal{Q}}={\omega }_{{\eta }_{\mathcal{Q}}}$ that is ${\alpha }^0$ is a symplectic diffeomorphism. Further, we can determine another potential one-form ${\theta }_{{\eta }_{\mathcal{Q}}}^{\prime }$ for the symplectic two-form ${\omega }_{{\eta }_{\mathcal{Q}}}$ given in (180) as follows

$${\theta }_{{\eta }_{\mathcal{Q}}}^{\prime }:={\left({\alpha }^0\right)}^{\ast }{\theta }_{\mathcal{T}\mathcal{Q}}=u{p}_{i}d{q}^i+{\dot{p}}_{i}d{\dot{q}}^i-udz,$$

(194)

where ${\theta }_{\mathcal{T}\mathcal{Q}}$ is the canonical one-form on the cotangent bundle $T^{\ast }\mathcal{T}\mathcal{Q}$. We now construct the left wing of evolution contact Tulczyjew’s triple as follows

(195)

where $\widetilde{\mathcal{T}{\pi }_{\mathcal{Q}}^0}$ is the projection defined by $\widetilde{\mathcal{T}{\pi }_{\mathcal{Q}}^0}=\mathcal{T}{\pi }_{\mathcal{Q}}^0\circ j={\pi }_{\mathcal{T}\mathcal{Q}}\circ {\alpha }^0$. The local expression of $\widetilde{\mathcal{T}{\pi }_{\mathcal{Q}}^0}$ is

$$\widetilde{\mathcal{T}{\pi }_{\mathcal{Q}}^0}:H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}⟶\mathcal{T}\mathcal{Q},(q^i,p_i,z,{\dot{q}}^i,{\dot{p}}_i,u)\mapsto (q^i,{\dot{q}}^i,z).$$

(196)

Note that, following the notation in (18), we have that the quintuple

$$(H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R},\widetilde{\mathcal{T}{\pi }_{\mathcal{Q}}^0},\mathcal{T}\mathcal{Q},-{\theta }_{{\eta }_{\mathcal{Q}}}^{\prime },{\alpha }^0)$$

(197)

is a special symplectic structure.

For a given Lagrangian function $L:\mathcal{T}\mathcal{Q}\mapsto \mathbb{R}$, ${\alpha }^0$ pulls the Lagrangian submanifold $im\left(dL\right)$ to a Lagrangian submanifold ${\left({\alpha }^0\right)}^{-1}\left(im\left(dL\right)\right)$ of the symplectic manifold $H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}$. This Lagrangian submanifold is the realization of the evolution Herglotz equations. In fact,

$${\left({\alpha }^0\right)}^{-1}\left(im\left(dL\right)\right)= \left\{(q^i,\frac{\partial L}{\partial {\dot{q}}^i},z,{\dot{q}}^i,\frac{\partial L}{\partial z}\frac{\partial L}{\partial {\dot{q}}^i}+\frac{\partial L}{\partial q^i},\dot{z},-\frac{\partial L}{\partial z})\in \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}:\dot{z}={\dot{q}}^i\frac{\partial L}{\partial {\dot{q}}^i}\right\}.$$

(198)

Thus, a curve $c^e:I\mapsto \mathcal{T}\mathcal{Q}$ is a solution of the evolution Herglotz equations for L if and only if ${(\mathbb{F}L\circ c^e)}^{\mathcal{T}}\left(t\right)\in {\left({\alpha }^0\right)}^{-1}\left(im\left(dL\right)\right)$, for every $t\in I$, where ${(\mathbb{F}L\circ c^e)}^{\mathcal{T}}:I\mapsto \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ is the lift of the curve $\mathbb{F}L\circ c^e:I\mapsto {\mathcal{T}}^{\ast }\mathcal{Q}$ given by (154).

We merge the right and left wings of triple given in (188) and (195), respectively. Thus, we have the evolution contact Tulczyjew’s triple:

(199)

Note that this triple is consisting of two special symplectic structures.

4.5. The Legendre Transformation

At first we introduce the following Whitney product

$$\mathcal{T}\mathcal{Q}{\times }_{\mathcal{Q}\times \mathbb{R}}{\mathcal{T}}^{\ast }\mathcal{Q}=\{(u,\zeta ,z)\in T\mathcal{Q}\times T^{\ast }\mathcal{Q}\times \mathbb{R}:{\tau }_{\mathcal{Q}}\left(u\right)={\pi }_{\mathcal{Q}}\left(\zeta \right)\}$$

(200)

over the extended manifold $\mathcal{Q}\times \mathbb{R}$. This is an $3n+1$ dimensional manifold which we call extended Pontryagin bundle. We introduce the local coordinates $(q^i,{\dot{q}}^i,p_i,z)$ on $\mathcal{T}\mathcal{Q}{\times }_{\mathcal{Q}\times \mathbb{R}}{\mathcal{T}}^{\ast }\mathcal{Q}$. We now generate the Legendrian submanifold ${\mathcal{N}}_L$ on $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ given in (168) referring to the right wing of the contact Tulczyjew’s triple (176). This is the Legendre transformation in the understanding of Tulczyjew. To have that, define the following Morse family

$$E:\mathcal{T}\mathcal{Q}{\times }_{\mathcal{Q}\times \mathbb{R}}{\mathcal{T}}^{\ast }\mathcal{Q}⟶\mathbb{R},(q^i,{\dot{q}}^i,p_i,z)\mapsto p_i{\dot{q}}^i-L(q,\dot{q},z).$$

(201)

To represent the minus of the Morse family $-E$, we are drawing the right side of the Tulczyjew’s triple (176) by equipping it with the Whitney product (200) as follows

(202)

This geometry fits well in the general picture in (137). Recalling (104), a direct computation determines the Legendrian submanifold of ${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}$ generated by minus of the Morse family $-E$ as

$$\begin{array}{cc}\hfill {\mathcal{N}}_{-E}& = \left\{(q^i,p_i,z,-\frac{\partial E}{\partial q^i},-\frac{\partial E}{\partial p_i},-\frac{\partial E}{\partial z},-E)\in {\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {\dot{q}}^i}=0\right\}\hfill \\ & = \left\{(q^i,p_i,z,\frac{\partial L}{\partial q^i},-{\dot{q}}^i,\frac{\partial L}{\partial z},-p_i{\dot{q}}^i+L)\in {\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}:p_i-\frac{\partial L}{\partial {\dot{q}}^i}=0\right\}.\hfill \end{array}$$

(203)

Using the inverse of the contact diffeomorphism ${\beta }^c$, we map the Legendrian submanifold ${\mathcal{N}}_{-E}$ to a Legendrian submanifold of the contact manifold $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ as

$${\left({\beta }^c\right)}^{-1}\left({\mathcal{N}}_{-E}\right)= \left\{\left(q^i,p_i,z,{\dot{q}}^i,p_i\frac{\partial L}{\partial z}+\frac{\partial L}{\partial q^i},L,-\frac{\partial L}{\partial z}\right) \in \mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}:p_i-\frac{\partial L}{\partial {\dot{q}}^i}=0\right\}$$

(204)

See that, this Legendrian submanifold is exactly the Legendrian submanifold ${\mathcal{N}}_L$ in (168) realizing the Herglotz Equations (128) and (129). This completes the Legendre transformation of the Herglotz equations to the contact Hamiltonian formalism.

If, further, the Lagrangian function $L=L(q,\dot{q},z)$ is non-degenerate, then from the equation

$$\frac{\partial E}{\partial {\dot{q}}^i}(q,\dot{q},p,z)=p_i-\frac{\partial L}{\partial {\dot{q}}^i}(q,\dot{q},z)=0$$

(205)

one can explicitly determine the velocity ${\dot{q}}^i$ in terms of $(q^i,p_i,z)$. This gives

$$\mathbb{F}L:\mathcal{T}\mathcal{Q}⟶{\mathcal{T}}^{\ast }\mathcal{Q},(q^i,{\dot{q}}^j,z)⟶ \left(q^i,\frac{\partial L}{\partial {\dot{q}}^j}(q,\dot{q},z),z\right)$$

(206)

as a local diffeomorphism relating $\mathcal{T}\mathcal{Q}$ and ${\mathcal{T}}^{\ast }\mathcal{Q}$. In this case, the Morse family E can be reduced to a well-defined contact Hamiltonian function

$$H(q,p,z)=p_i{\dot{q}}^i(q,p,z)-L\left(q,\dot{q}(q,p,z),z\right)$$

(207)

on ${\mathcal{T}}^{\ast }\mathcal{Q}$.

It is possible to perform the inverse Legendre transformation of the contact Hamiltonian dynamics as well. This time, one needs to generate the Legendrian submanifold ${\mathcal{N}}_{-H}$ in (152) realizing the contact Hamilton’s equations (156) referring to the left wing of the contact triple (176).

The Legendre Transformation for Evolution Dynamics. Recall the Tulczyjew’s triple (199) exhibited for the case of evolution contact dynamics. We consider the Lagrangian submanifold ${\left({\alpha }^0\right)}^{\ast }\left(im\left(dL\right)\right)$ of $H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}$ generated by a Lagrangian $L=L(q,\dot{q},z)$ on $\mathcal{T}\mathcal{Q}$ referring to the left wing (195) of the triple. Once more, consider the total space given in (200), and the energy function E given in (201). In this evolution case, we plot the following diagram merging the right wing (188) of the evolution Tulczyjew’s triple and the Morse family determined by $-E$ that is

(208)

From (30), we deduce that the Lagrangian submanifold ${\mathcal{S}}_{-E}$ of the cotangent bundle $T^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}$ generated by $-E$ is computed to be

$$\begin{array}{cc}\hfill {\mathcal{S}}_{-E}& = \left\{(q^i,p_i,z,-\frac{\partial E}{\partial q^i},-\frac{\partial E}{\partial p_i},-\frac{\partial E}{\partial z})\in T^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}:\frac{\partial E}{\partial {\dot{q}}^i}=0\right\}\hfill \\ & = \left\{(q^i,p_i,z,\frac{\partial L}{\partial q^i},-{\dot{q}}^i,\frac{\partial L}{\partial z})\in {\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }\mathcal{Q}:p_i-\frac{\partial L}{\partial {\dot{q}}^i}=0\right\}.\hfill \end{array}$$

(209)

Using the inverse of the symplectic diffeomorphism ${\beta }^0$, we transfer the Lagrangian submanifold ${\mathcal{S}}_{-E}$ to a Lagrangian submanifold of $H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}$ as follows

$$\begin{array}{cc}\hfill {\left({\beta }^0\right)}^{-1}\left({\mathcal{S}}_{-E}\right)=& \{\left(q^i,p_i,z,{\dot{q}}^i,p_i\frac{\partial L}{\partial z}+\frac{\partial L}{\partial q^i},\dot{z},-\frac{\partial L}{\partial z}\right) \in H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}:\hfill \\ & p_i-\frac{\partial L}{\partial {\dot{q}}^i}=0,\dot{z}-{\dot{q}}^i\frac{\partial L}{\partial {\dot{q}}^i}=0\}.\hfill \end{array}$$

This is exactly the Lagrangian submanifold ${\left({\alpha }^0\right)}^{-1}\left(im\left(dL\right)\right)$ realizing the evolution Herglotz equations. In a similar way, one may obtain the inverse Legendre transformation of the contact evolution dynamics for a Hamiltonian function $H:{\mathcal{T}}^{\ast }\mathcal{Q}\mapsto \mathbb{R}$.

5. Example: The Ideal Gas

5.1. A Quantomorphism on the Euclidean Space

Thermodynamics have been studied extensively in the framework of contact geometry. For some recent work directly related with the present discussions, we cite [35,43,44,45]. In this section, we shall be applying the theoretical results obtained in the previous sections to some thermodynamical models.

We start this subsection by providing the following theorem realizing a strict contact diffeomorphism (quantomorphism) on the extended cotangent bundle ${\mathcal{T}}^{\ast }{\mathbb{R}}^m$ of the Euclidean space, see also [43]. The proof follows by a direct calculation.

Theorem 6.

Consider a disjoint partition $I\cup J$ of the set of indices $\{1,\dots ,m\}$ so that the coordinates on ${\mathbb{R}}^m$ is given as $(x^a,x^{\rho })$, where $a\in I$ and $\rho \in J$. Then the following mapping

$$\varphi :{\mathcal{T}}^{\ast }{\mathbb{R}}^m⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^m,(x^a,x^{\rho },y_a,y_{\rho },u)\mapsto (x^a,y_{\rho },y_a,-x^{\rho },u-x^{\rho }{y}_{\rho })$$

(210)

preserves the canonical contact one-form ${\eta }_{{\mathbb{R}}^m}=du-y_{a}d{x}^a-y_{\rho }d{x}^{\rho }$. Here, $(x^a,x^{\rho },y_a,y_{\rho },u)$ are the Darboux’s coordinates on the extended cotangent bundle ${\mathcal{T}}^{\ast }{\mathbb{R}}^m$.

In Section 3.2, we have stated that the image of the first prolongation of a smooth function on the base manifold is a Legendrian submanifold of the extended cotangent bundle. Accordingly, consider a smooth function $U=U(x^a,x^{\rho })$ on ${\mathbb{R}}^m$ so that its first prolongation ${\mathcal{T}}^{\ast }U$ to the extended cotangent bundle turns out to be a Legendrian submanifold of ${\mathcal{T}}^{\ast }{\mathbb{R}}^m$ as given in (103). Under the quantomorphism $\varphi$ in (210), we have a Legendrian submanifold on the image space as

$$\varphi (x^a,x^{\rho },\frac{\partial U}{\partial x^a},\frac{\partial U}{\partial x^{\rho }},U)=(x^a,\frac{\partial U}{\partial x^{\rho }},\frac{\partial U}{\partial x^a},-x^{\rho },U-x^{\rho }\frac{\partial U}{\partial x^{\rho }}).$$

(211)

This alternative realization of the Legendrian submanifold is important for geometric characterization of reversible thermodynamics.

Remark 1.

If $\mathcal{Q}={\mathbb{R}}^n$ then the extended cotangent and the extended tangent bundles turn out to be isomorphic that is ${\mathcal{T}}^{\ast }\mathcal{Q}\simeq \mathcal{T}\mathcal{Q}\simeq {\mathbb{R}}^{2n+1}$. In this particular instance, by assuming $m=2n+1$, the mapping ${\psi }^c$ in (159) can be regarded as a particular case of ϕ in (210) if the canonical coordinates $(q^i,p_i,z)$ are decomposed as $\left(x^a\right)=(q^i,z)$ where $a=1,\dots ,n+1$, and $\left(x^{\rho }\right)=\left(p_i\right)$ where $\rho =1,\dots ,n$.

5.2. Equilibrium Thermodynamics

Obeying the geometry exhibited in the previous section, we take $m=3$ with coordinates $(S,V,N)\in {\mathbb{R}}^3$. Here, S stands for the entropy, V is the volume, and N is the mole number of classical ideal gas. The conjugate variables $(T,-P,\mu )\in {\left({\mathbb{R}}^3\right)}^{\ast }$ are the temperature, the pressure, and the chemical potential, respectively. By employing the internal energy U as the fiber coordinate, we complete the following realization of the extended cotangent bundle $(S,V,N,T,-P,\mu ,U)\in {\mathcal{T}}^{\ast }{\mathbb{R}}^3$. Consider the contact one-form

$${\eta }_{{\mathbb{R}}^3}=dU-TdS+PdV-\mu dN.$$

(212)

As a particular instance, we choose the internal energy

$$U(S,V,N)=U_0{V}^{-1/c}{N}^{(c+1)/c}exp\left(\frac{S}{cNR}\right),$$

(213)

as a function depending on the base coordinates $(S,V,N)\in {\mathbb{R}}^3$. Here, $U_0$ is a positive constant, c is the heat capacity and R is the universal gas constant. The first prolongation ${\mathcal{T}}^{\ast }U$ is a Legendrian submanifold $\mathcal{N}$ of the contact manifold $({\mathcal{T}}^{\ast }{\mathbb{R}}^3,{\eta }_{{\mathbb{R}}^3})$. By considering that the temperature $T=\partial U/\partial S$ and the pressure $P=-\partial U/\partial V$, we have the following set of equations

$$c{V}^{1/c}RT=U_0{N}^{1/c}exp\left(\frac{S}{cNR}\right),PV=NRT,\mu =(c+1)RT-TS/N$$

(214)

those realizing $\mathcal{N}$.

The Legendre Transformations. In the light of Theorem 6, and the transformation (211), we now present the Legendre transformation between the internal energy, the enthalpy, the Helmholtz function, and the Gibbs function. For a similar discussion, but in the framework of symplectic geometry, see [12]. We start with the Legendrian submanifold determined by the internal energy U in (213).

(1) We decompose the base variables as $(S,N)$ and V and apply Theorem 6 to the volume variable. This results with a quantomorphism computed to be

$${\varphi }^1:{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,(S,V,N,T,-P,\mu ,U)\mapsto (S,-P,N,T,-V,\mu ,U+PV).$$

(215)

Note that, on the image space, the fiber component is the enthalpy function $B=U+PV$. If we solve the pressure from the equation $P=-\partial U/\partial V$, the enthalpy function can be written as a function of the new base variables $(S,P,N)$, that is

$$B(S,P,N)=\overline{c}{U}_0^{(c+1)/c}{P}^{1/(1+c)}Nexp\left(\frac{S}{(c+1)NR}\right),$$

(216)

where $\overline{c}$ is a constant defined to be $c^{1/(1+c)}+c^{-c/(1+c)}$. So that, the enthalpy function is another generator of the same Legendrian submanifold. Indeed, the first prolongation of ${\mathcal{T}}^{\ast }B$ is given the system of equations in (214) so that ${\varphi }^1\circ {\mathcal{T}}^{\ast }B=\mathcal{N}$.

(2) We start once more with the internal energy, but this time we perform the transformation to the entropy variable S. For this case, we have the quantomorphism

$${\varphi }^2:{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,(S,V,N,T,-P,\mu ,U)\mapsto (T,V,N,-S,-P,\mu ,U-ST).$$

(217)

In this case, the fiber term $F=U-ST$ is the Helmholtz function. Using the identity $T=\partial U/\partial S$, we write the Helmholtz function as a function of the base components $(T,V,N)$ of the image space, that is

$$F(T,V,N)=cNRT\left(1+\frac{1}{c}logN-\frac{1}{c}logV+log\left(\frac{U_0}{cRT}\right)\right).$$

(218)

Thus, F is another generator of the same Legendrian submanifold determined by the equations (214), that is ${\varphi }^2\circ {\mathcal{T}}^{\ast }F=\mathcal{N}$.

(3) This time, we consider the Helmholtz function F in (218) and apply the transformation given in Theorem 6 to the volume variable, that is

$${\varphi }^3:{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,(T,V,N,-S,-P,\mu ,F)\mapsto (T,-P,N,-S,-V,\mu ,F+PV).$$

(219)

The fiber term $G=F+PV$ is the Gibbs function. By taking $P=-\partial F/\partial V$, we write the Gibbs function as a function of the base variables, that is

$$G(T,P,N)=NRT\left(1+c+logN-log\frac{NRT}{P}+clog\left(\frac{U_0}{cRT}\right)\right).$$

(220)

We have that ${\varphi }^3\circ {\mathcal{T}}^{\ast }G=\mathcal{N}$ is the Legendrian submanifold determined by Equation (214). Evidently, by iteratively applying ${\varphi }^2$ in (217) and ${\varphi }^3$ in (219), one can define a quantomorphism from the internal energy setting to the Gibbs function setting. In this case, the quantomorphism is determined as

$$\begin{array}{cc}\hfill & {\varphi }^3\circ {\varphi }^2:{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,\hfill \\ & (S,V,N,T,-P,U)\mapsto (T,-P,N,-S,-V,\mu ,U-TS+PV).\hfill \end{array}$$

(221)

This gives a direct passage from the internal energy to the Gibbs function.

(4) To complete the Legendre transformations, we consider the Gibbs function in (220) and define the quantomorphism to the chemical potential variable $\mu$, that is

$$\begin{array}{cc}\hfill & {\varphi }^4:{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,\hfill \\ & (T,-P,N,-S,-V,\mu ,G)\mapsto (T,-P,\mu ,-S,-V,-N,G-\mu N).\hfill \end{array}$$

(222)

The fiber variable $W=G-\mu N$ determines a new generator of the Legendrian submanifold $\mathcal{N}$ depending on the base variables of the image space. In this case, we have that

$$\begin{array}{cc}\hfill W(T,P,\mu )=& NRT\left(1+c+log\frac{ST}{(c+1)RT-\mu }-log\frac{NRT}{P}+clog\left(\frac{U_0}{cRT}\right)\right)\hfill \\ & -\frac{\mu ST}{(c+1)RT-\mu }.\hfill \end{array}$$

(223)

Composing the quantomorphism ${\varphi }^3\circ {\varphi }^2$ in (221) with ${\varphi }^4$ in (222), one can derive a Legendre transformation from the internal energy U setting to W as

$$\begin{array}{cc}\hfill & \varphi :{\mathcal{T}}^{\ast }{\mathbb{R}}^3⟶{\mathcal{T}}^{\ast }{\mathbb{R}}^3,\hfill \\ & (S,V,N,T,-P,\mu ,U)\mapsto (T,-P,\mu ,-S,-V,-N,U-TS+PV-\mu N).\hfill \end{array}$$

(224)

This is the full Legendre transformation of the ideal gas. If we insist that the dynamics of the gas stays in the Legendrian submanifold $\mathcal{N}$ determined through the equations (214), we can consider the fiber variable $U-TS+PV-\mu N$ for an arbitrary internal energy U as a Hamiltonian function of the dynamics for the ideal case. Let us discuss this motion in the following subsection.

5.3. Hamiltonian Flow and Its Legendrian Realization

We consider the quantomorphism in (224), but instead of a specific internal energy, we let U be an independent variable. Consider the Hamiltonian function

$$H(S,V,N,T,-P,\mu ,U)=TS-NRT+\mu N-U$$

(225)

on the extended cotangent bundle ${\mathcal{T}}^{\ast }{\mathbb{R}}^3$. Note that we determine the Hamiltonian function (225) by substituting the ideal gas equation $PV=NRT$ (that is the second equation in (214)) into the minus of the generator function $U-TS+PV-\mu N$ of the quantomorphism. The minus sign is to fit physical intuition, for example, to be sure that the entropy S is increasing along the motion. In fact, while computing the Hamiltonian dynamics, the minus sign will be compensated with the minus sign appearing in the formulation given in (147). According to (151), the image space of the first prolongation of the Hamiltonian function (225)

$$\begin{array}{cc}\hfill & im(-{\mathcal{T}}^{\ast }H)=\hfill \\ & \left\{(S,V,N,T,-P,\mu ,U,-T,0,RT-\mu ,-S+NR,0,-N,1,-H)\in {\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }{\mathbb{R}}^3\right\}\hfill \end{array}$$

(226)

is a Legendrian submanifold of the iterated extended cotangent bundle ${\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }{\mathbb{R}}^3$. On the other hand, the contact Hamiltonian vector field associated with the Hamiltonian function (225) is

$$X_H^c=(S-NR)\frac{\partial }{\partial S}+N\frac{\partial }{\partial N}+P\frac{\partial }{\partial P}+RT\frac{\partial }{\partial \mu }+U\frac{\partial }{\partial U}.$$

(227)

Notice that, the coefficients of $\partial /\partial T$ and $\partial /\partial V$ are zero. This gives that the Hamiltonian dynamics is isothermal and isochoric.

We consider the induced coordinates on the extended tangent bundle $\mathcal{T}{\mathcal{T}}^{\ast }{\mathbb{R}}^3$ as follows. For the base manifold we use $(S,V,N,T,-P,\mu ,U)$, for the fibers of the tangent bundle $T{\mathcal{T}}^{\ast }{\mathbb{R}}^3$ we refer $(\dot{S},\dot{V},\dot{N},\dot{T},-\dot{P},\dot{\mu },\dot{U})$ and for the extension $\mathbb{R}$ we use u as the standard coordinate on $\mathbb{R}$. Accordingly, for the present case, the lifted contact one-form (148) is computed to be

$${\eta }^{\mathcal{T}}=d\dot{U}+udU-(\dot{T}+uT)dS+(\dot{P}+uP)dV-(\dot{\mu }+u\mu )dN-Td\dot{S}+Pd\dot{V}-\mu d\dot{N}.$$

(228)

The Reeb field on ${\mathcal{T}}^{\ast }{\mathbb{R}}^3$ is $\partial /\partial U$, and the directional derivative of the Hamiltonian function (225) is $\mathcal{R}\left(H\right)=-1$. The image space of the coupling of the Hamiltonian vector field and $\mathcal{R}\left(H\right)$ determines a Legendrian submanifold, of the extended tangent bundle $\mathcal{T}{\mathcal{T}}^{\ast }{\mathbb{R}}^3$, given by

$$\begin{array}{cc}\hfill {\mathcal{N}}_{-H}& =im(X_H^c,\mathcal{R}\left(H\right))\hfill \\ & = \left\{(S,V,N,T,-P,\mu ,U;S-NR,0,N,0,-P,RT,U;-1)\in \mathcal{T}{\mathcal{T}}^{\ast }{\mathbb{R}}^3\right\}.\hfill \end{array}$$

(229)

Evidently, the Legendrian submanifolds in (226) and (229) are related with the contactomorphism ${\beta }^c:\mathcal{T}{\mathcal{T}}^{\ast }{\mathbb{R}}^3\mapsto {\mathcal{T}}^{\ast }{\mathcal{T}}^{\ast }{\mathbb{R}}^3$ by satisfying the relation ${\beta }^c\circ {\mathcal{N}}_{-H}=im(-{\mathcal{T}}^{\ast }H)$ as we have proven in (147).

Now, we wish to generate the Legendrian submanifold in (229) referring to the left wing of the contact Tulczyjew’s triple (176), that is to generate it via a Lagrangian function (probably as a Morse family) defined on the extended tangent bundle $\mathcal{T}{\mathbb{R}}^3$. In the light of Section 4.5, we now apply the inverse Legendre transformation. The first step is to apply ${\alpha }^c$ in (162) to the Legendrian submanifold ${\mathcal{N}}_{-H}$ in (229). The image space is

$${\alpha }^c\left({\mathcal{N}}_{-H}\right)= \left\{(S,V,N,S-NR,0,N,U;-T,0,RT-\mu ,T,-P,\mu ,1;U)\in {\mathcal{T}}^{\ast }\mathcal{T}{\mathbb{R}}^3\right\}$$

(230)

which is a Legendrian submanifold of ${\mathcal{T}}^{\ast }\mathcal{T}{\mathbb{R}}^3$.

We assume the coordinates $(S,V,N,\dot{S},\dot{V},\dot{N},U)$ on the extended tangent bundle $\mathcal{T}{\mathbb{R}}^3$ and define the Whitney sum of the extended tangent and the extended cotangent bundles $\mathcal{T}{\mathbb{R}}^3{\times }_{{\mathbb{R}}^3\times \mathbb{R}}{\mathcal{T}}^{\ast }{\mathbb{R}}^3$ with coordinates $(S,V,N,\dot{S},\dot{V},\dot{N},T,-P,\mu ,U)$. Note that, in the Whitney sum, we fix the base coordinates $(S,V,N)$ in ${\mathbb{R}}^3$ and the extension U in $\mathbb{R}$. The subscript ${\mathbb{R}}^3\times \mathbb{R}$ in the notation of Whitney sum manifests these choices. A calculation gives that the Legendrian submanifold ${\alpha }^c\left({\mathcal{N}}_{-H}\right)$ in (230) is generated by a Morse family on $\mathcal{T}{\mathbb{R}}^3{\times }_{{\mathbb{R}}^3\times \mathbb{R}}{\mathcal{T}}^{\ast }{\mathbb{R}}^3$. In other words, the Lagrangian function

$$L(S,V,N,\dot{S},\dot{V},\dot{N},U;T,-P,\mu )=T(\dot{S}-S+NR)+\mu (\dot{N}-N)+P\dot{V}+U.$$

(231)

can be understood as defined on the extended tangent bundle $\mathcal{T}{\mathbb{R}}^3$, but depending on the auxiliary variables $(T,-P,\mu )$. By merging the Lagrangian function with the left wing of the contact Tulzyjew’s triple (176), we get the following diagram:

(232)

According to the local realization of the generating family in (104), it is a direct calculation to show that the Legendrian submanifold ${\alpha }^c\left({\mathcal{N}}_{-H}\right)$ is generated by the Morse family L in (231).

We remark that in this formalism of thermodynamics, Hamiltonians are usually singular, as in the case above, so there is no Lagrangian formulation in the classical sense. Hence, we think that Tulcyzjew triples might be a useful tool in this situation.

5.4. Evolutionary Flow and Its Lagrangian Realization

We once more consider the Hamiltonian function H exhibited in (225) and defined on the extended cotangent bundle ${\mathcal{T}}^{\ast }{\mathbb{R}}^3$. The evolutionary vector field ${\epsilon }_H$ is defined in (119). A direct calculation determines the evolutionary vector field for the Hamiltonian function (225) as

$${\epsilon }_H=(S-NR)\frac{\partial }{\partial S}+N\frac{\partial }{\partial N}+P\frac{\partial }{\partial P}+RT\frac{\partial }{\partial \mu }+(TS-NRT+\mu N)\frac{\partial }{\partial U}$$

(233)

on ${\mathcal{T}}^{\ast }{\mathbb{R}}^3$. By referring to Corollary 1, we establish that the image space

$$\begin{array}{cc}\hfill im({\epsilon }_H,\mathcal{R}\left(H\right))=& \left\{\right(S,V,N,T,-P,\mu ,U;\hfill \\ & S-NR,0,N,0,-P,RT,TS-NRT+\mu N;-1)\in \mathcal{T}{\mathcal{T}}^{\ast }{\mathbb{R}}^3\}\hfill \end{array}$$

(234)

turns out to be a Lagrangian submanifold of $H{\mathcal{T}}^{\ast }{\mathbb{R}}^3\times \mathbb{R}$ defined in (179). Via ${\alpha }^0$ in (193), we map the Lagrangian submanifold (234) to a Lagrangian submanifold

$${\alpha }^0\left(im({\epsilon }_H,\mathcal{R}\left(H\right))\right) = \left\{(S,V,N,S-NR,0,N,U;-T,0,RT-\mu ,T,-P,\mu ,1)\in T^{\ast }\mathcal{T}{\mathbb{R}}^3\right\}$$

(235)

in the cotangent bundle $T^{\ast }\mathcal{T}{\mathbb{R}}^3$. According to (30), it is immediate to see that the Lagrangian function L given in (231) defined on the Whitney sum of the extended tangent and the extended cotangent bundles $\mathcal{T}{\mathbb{R}}^3{\times }_{{\mathbb{R}}^3\times \mathbb{R}}{\mathcal{T}}^{\ast }{\mathbb{R}}^3$ generates the Lagrangian submanifold in (235). In order to visualize the Lagrangian submanifold, we draw the following diagram by merging the Morse family (231) and the left wing of the evolution contact Tulczyjew’s triple (199)

(236)

where $\tilde{\mathcal{T}{\pi }_{{\mathbb{R}}^3}^0}$ is the mapping given in (196).

6. Discussion

In this paper, we have used the tangent contact structure, on the extended tangent bundle $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$, which was introduced in [70]. Referring to this, and by introducing the notion of special contact structure, we have constructed a Tulcyzjew’s triple for contact manifolds, see Diagram (176). This permits us to describe both the contact Lagrangian and the contact Hamiltonian dynamics as Legendrian submanifolds of $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$. In this formulation, the Legendre transformation is defined as a passage between two generators of the same Legendrian submanifold. Note that, this approach is free from the Hessian condition. That means, it is applicable for degenerate theories as well. We, further, present Tulcyzjew’s triple for evolutionary dynamics, see Diagram (199). Instead of contact structures, the evolution triple (199) consists of special symplectic structures. In this construction, the contact manifold $\mathcal{T}{\mathcal{T}}^{\ast }\mathcal{Q}$ is substituted by the extended horizontal bundle $H{\mathcal{T}}^{\ast }\mathcal{Q}\times \mathbb{R}$, which is symplectic. We have concluded the paper by applications of the theoretical results to geometrical foundations of some thermodynamical models.

Here are some further questions we wish to pursue:

• In Section 4.2, we have established that the image space of a contact Hamiltonian vector field is a Legendrian submanifold of the tangent contact manifold. Evidently, not all Legendrian submanifolds determine explicit dynamical equations. This observation motivates us to define the notion of an implicit Hamiltonian Contact Dynamics as a non-horizontal Legendrian submanifold of the tangent contact manifold. We refer to [75] for a similar discussion done for the case of symplectic dynamics and integrability of the non-horizontal Lagrangian submanifolds. We find it interesting to elaborate the integrability of implicit Hamiltonian contact dynamics.
• Following, the first question raised in this section, we plan to write a Hamilton–Jacobi theory for implicit Hamiltonian contact dynamics. Hamilton–Jacobi theory for (explicit) Hamiltonian contact dynamics is recently examined in [72,76]. Hamilton–Jacobi theory for implicit symplectic dynamics is discussed in [77,78].
• In the literature, Tulczyjew’s triple for higher order classical dynamical systems is already available [14,15]. Higher order contact dynamics is studied in [79]. As a future work, we plan to extend the geometry presented in the present paper to higher order contact framework.