Hamiltonian Flows Research Articles

Proof of bulk-edge correspondence for band topology by Toeplitz algebra

Abstract We rigorously yet concisely prove the bulk-edge correspondence for general d-dimensional (dD) topological insulators in complex Altland-Zirnbauer classes, which states that the bulk topological number equals to the edge-mode index. Specifically, an essential formula is discovered that links the quantity expressed by Toeplitz algebra, i.e., hopping terms on the lattice with an edge, to the Fourier series on the bulk Brillouin zone. We then apply it to chiral models and utilize exterior differential calculations, instead of the sophisticated K -theory, to show that the winding number of bulk system equals to the Fredholm index of 1D edge Hamiltonian, or to the sum of edge winding numbers for higher odd dimensions. Moreover, this result is inherited to the even dimensional Chern insulators as each of them can be mapped to an odd dimensional chiral model. It is revealed that the Chern number of bulk system is identical to the spectral flow of 2D edge Hamiltonian, or to the negative sum of edge Chern numbers for higher even dimensions. Our methods and conclusions are friendly to physicists and could be easily extended to other physical scenarios.

Symmetry Preservation in Hamiltonian Systems: Simulation and Learning

This work presents a general geometric framework for simulating and learning the dynamics of Hamiltonian systems that are invariant under a Lie group of transformations. This means that a group of symmetries is known to act on the system respecting its dynamics and, as a consequence, Noether’s theorem, conserved quantities are observed. We propose to simulate and learn the mappings of interest through the construction of G-invariant Lagrangian submanifolds, which are pivotal objects in symplectic geometry. A notable property of our constructions is that the simulated/learned dynamics also preserves the same conserved quantities as the original system, resulting in a more faithful surrogate of the original dynamics than non-symmetry aware methods, and in a more accurate predictor of non-observed trajectories. Furthermore, our setting is able to simulate/learn not only Hamiltonian flows, but any Lie group-equivariant symplectic transformation. Our designs leverage pivotal techniques and concepts in symplectic geometry and geometric mechanics: reduction theory, Noether’s theorem, Lagrangian submanifolds, momentum mappings, and coisotropic reduction among others. We also present methods to learn Poisson transformations while preserving the underlying geometry and how to endow non-geometric integrators with geometric properties. Thus, this work presents a novel attempt to harness the power of symplectic and Poisson geometry toward simulating and learning problems.

Generalized Hamiltonian and Lagrangian aspects of a model for virus–tumor interaction in oncolytic virotherapy

We analyze the generalized Hamiltonian structure of a system of first‐order ordinary differential equations for the Jenner et al. system (Letters in Biomathematics 5 (2018), no. S1, S117–S136). The system of equations is used for modeling the interaction of an oncolytic virus with a tumor cell population. Our analysis is based on the existence of a Jacobi last multiplier and a time‐dependent first integral. Suitable conditions on the model parameters allow for the reduction of the problem to a planar system of equations, and the time‐dependent Hamiltonian flows are described. The geometry of the Hamiltonian flows is also investigated using the symplectic and cosymplectic methods.

Understanding of linear operators through Wigner analysis

In this work, we extend Wigner's original framework to analyze linear operators by examining the relationship between their Wigner and Schwartz kernels. Our approach includes the introduction of (quasi-)algebras of Fourier integral operators (FIOs), which encompass FIOs of type I and II. The symbols of these operators belong to (weighted) modulation spaces, particularly in Sjöstrand's class, known for its favorable properties in time-frequency analysis. One of the significant results of our study is demonstrating the inverse-closedness of these symbol classes.Our analysis includes fundamental examples such as pseudodifferential operators and Fourier integral operators related to Schrödinger-type equations. These examples typically feature classical Hamiltonian flows governed by linear symplectic transformations S∈Sp(d,R). The core idea of our approach is to utilize the Wigner kernel to transform a Fourier integral operator T on Rd into a pseudodifferential operator K on R2d. This transformation involves a symbol σ well-localized around the manifold defined by z=Sw.

From Geometry of Hamiltonian Dynamics to Topology of Phase Transitions: A Review.

In this review work, we outline a conceptual path that, starting from the numerical investigation of the transition between weak chaos and strong chaos in Hamiltonian systems with many degrees of freedom, comes to highlight how, at the basis of equilibrium phase transitions, there must be major changes in the topology of submanifolds of the phase space of Hamiltonian systems that describe systems that exhibit phase transitions. In fact, the numerical investigation of Hamiltonian flows of a large number of degrees of freedom that undergo a thermodynamic phase transition has revealed peculiar dynamical signatures detected through the energy dependence of the largest Lyapunov exponent, that is, of the degree of chaoticity of the dynamics at the phase transition point. The geometrization of Hamiltonian flows in terms of geodesic flows on suitably defined Riemannian manifolds, used to explain the origin of deterministic chaos, combined with the investigation of the dynamical counterpart of phase transitions unveils peculiar geometrical changes of the mechanical manifolds in correspondence to the peculiar dynamical changes at the phase transition point. Then, it turns out that these peculiar geometrical changes are the effect of deeper topological changes of the configuration space hypersurfaces ∑v=VN-1(v) as well as of the manifolds {Mv=VN-1((-∞,v])}v∈R bounded by the ∑v. In other words, denoting by vc the critical value of the average potential energy density at which the phase transition takes place, the members of the family {∑v}v<vc are not diffeomorphic to those of the family {∑v}v>vc; additionally, the members of the family {Mv}v>vc are not diffeomorphic to those of {Mv}v>vc. The topological theory of the deep origin of phase transitions allows a unifying framework to tackle phase transitions that may or may not be due to a symmetry-breaking phenomenon (that is, with or without an order parameter) and to finite/small N systems.

Learning Hamiltonian dynamics with reproducing kernel Hilbert spaces and random features

A method for learning Hamiltonian dynamics from a limited and noisy dataset is proposed. The method learns a Hamiltonian vector field on a reproducing kernel Hilbert space (RKHS) of inherently Hamiltonian vector fields, and in particular, odd Hamiltonian vector fields. This is done with a symplectic kernel, and it is shown how the kernel can be modified to an odd symplectic kernel to impose the odd symmetry. A random feature approximation is developed for the proposed odd kernel to reduce the problem size. The performance of the method is validated in simulations for three Hamiltonian systems. It is demonstrated that the use of an odd symplectic kernel improves prediction accuracy and data efficiency, and that the learned vector fields are Hamiltonian and exhibit the imposed odd symmetry characteristics.

Phase-space path integral approach to the kinetics of black hole phase transition in massive gravity

The dynamics of the state-switching process of black holes in dRGT massive gravity theory is presented using free energy landscape and stochastic Langevin equations. The free energy landscape is constructed using the Gibbons-Hawking path integral method. The black hole phases are characterized by taking its horizon radius as the order parameter. The free energy landscape provides three black hole phases: small, intermediate, and large. The small and large black holes are thermodynamically stable whereas the intermediate one is unstable. The Martin–Siggia–Rose–Janssen–de Dominicis (MSRJD) functional describes the stochastic dynamics of black hole phase transition. The Hamiltonian flow lines are obtained from the MSRJD functional and are used to analyze the stability and the phase transition properties. The dominant kinetic path between different phases is discussed for various configurations of the free energy landscape. We discuss the effect of black hole charge and the graviton mass on the critical behavior of black hole phase transition.

Cosymplectic Geometry, Reductions, and Energy-Momentum Methods with Applications

Classical energy-momentum methods study the existence and stability properties of solutions of t-dependent Hamilton equations on symplectic manifolds whose evolution is given by their Hamiltonian Lie symmetries. The points of such solutions are called relative equilibrium points. This work devises a new cosymplectic energy-momentum method providing a new and more general framework to study t-dependent Hamilton equations. In fact, cosymplectic geometry allows for using more types of distinguished Lie symmetries (given by Hamiltonian, gradient, or evolution vector fields), relative equilibrium points, and reduction methods, than symplectic techniques. To make our work more self-contained and to fill some gaps in the literature, a review of the cosymplectic formalism and the cosymplectic Marsden–Weinstein reduction is included. Known and new types of relative equilibrium points are characterised and studied. Our methods remove technical conditions used in previous energy-momentum methods, like the Ad∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extrm{Ad}^*$$\\end{document}-equivariance of momentum maps. Eigenfunctions of t-dependent Schrödinger equations are interpreted in terms of relative equilibrium points in cosymplectic manifolds. A new cosymplectic-to-symplectic reduction is developed and a new associated type of relative equilibrium points, the so-called gradient relative equilibrium points, are introduced and applied to study the Lagrange points and Hill spheres of a restricted circular three-body system by means of a not Hamiltonian Lie symmetry of the system.

A Resonant Lyapunov Centre Theorem with an Application to Doubly Periodic Travelling Hydroelastic Waves

We present a Lyapunov centre theorem for an antisymplectically reversible Hamiltonian system exhibiting a nondegenerate 1 : 1 or 1:-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1:-1$$\\end{document} semisimple resonance as a detuning parameter is varied. The system can be finite- or infinite-dimensional (and quasilinear) and have a non-constant symplectic structure. We allow the origin to be a ‘trivial’ eigenvalue arising from a translational symmetry or, in an infinite-dimensional setting, to lie in the continuous spectrum of the linearised Hamiltonian vector field provided a compatibility condition on its range is satisfied. As an application, we show how Kirchgässner’s spatial dynamics approach can be used to construct doubly periodic travelling waves on the surface of a three-dimensional body of water (of finite or infinite depth) beneath a thin ice sheet (‘hydroelastic waves’). The hydrodynamic problem is formulated as a reversible Hamiltonian system in which an arbitrary horizontal spatial direction is the time-like variable, and the infinite-dimensional phase space consists of wave profiles which are periodic (with fixed period) in a second, different horizontal direction. Applying our Lyapunov centre theorem at a point in parameter space associated with a 1 : 1 or 1:-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1:-1$$\\end{document} semisimple resonance yields a periodic solution of the spatial Hamiltonian system corresponding to a doubly periodic hydroelastic wave.

Enhanced Dissipation for Two-Dimensional Hamiltonian Flows

Let H∈C1∩W2,p\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H\\in C^1\\cap W^{2,p}$$\\end{document} be an autonomous, non-constant Hamiltonian on a compact 2-dimensional manifold, generating an incompressible velocity field b=∇⊥H\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$b=\ abla ^\\perp H$$\\end{document}. We give sharp upper bounds on the enhanced dissipation rate of b in terms of the properties of the period T(h) of the closed orbit {H=h}\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\{H=h\\}$$\\end{document}. Specifically, if 0<ν≪1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0<\ u \\ll 1$$\\end{document} is the diffusion coefficient, the enhanced dissipation rate can be at most O(ν1/3)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$O(\ u ^{1/3})$$\\end{document} in general, the bound improves when H has isolated, non-degenerate elliptic points. Our result provides the better bound O(ν1/2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$O(\ u ^{1/2})$$\\end{document} for the standard cellular flow given by Hc(x)=sinx1sinx2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_{\ extsf{c}}(x)=\\sin x_1 \\sin x_2$$\\end{document}, for which we can also prove a new upper bound on its mixing rate and a lower bound on its enhanced dissipation rate. The proofs are based on the use of action-angle coordinates and on the existence of a good invariant domain for the regular Lagrangian flow generated by b.

Twisting in Hamiltonian flows and perfect fluids

Twisting in Hamiltonian flows and perfect fluids

Variational aspect and kinetic theory of locally conformal dynamics

Abstract We present the locally conformal generalization of the Euler–Lagrange equations. We determine the dual space of the LCS Hamiltonian vector fields. Within this dual space, we formulate the Lie–Poisson equation that governs the kinetic motion of Hamiltonian systems in the context of local conformality. By expressing the Lie–Poisson dynamics in terms of density functions, we derive locally conformal Vlasov dynamics. In addition, we outline a geometric pathway that connects LCS Hamiltonian particle motion to locally conformal kinetic motion.

Generalized Hamiltonian systems on subvarieties: constant rank case

For the constraint variety in symplectic manifold, the solvable Hamiltonian vector fields on the constraint are investigated. According to P.A.M. Dirac [3], the space of solvable Hamiltonian systems is determined by the geometric restriction of the symplectic form to the constraint. Solvability condition of the generalized Hamiltonian systems is extended to singular varieties and applied under some assumption on singularities. The constraint being a smooth submanifold in a symplectic space was considered in [6]. In this paper, we investigate the solvability of generalized Hamiltonian systems and the constraint invariants on singular constraints in the constant rank case.

Invariant manifolds in a reversible Hamiltonian system: The tentacle-like geometry

We study a one-parameter family of time-reversible Hamiltonian vector fields in R4, which has received great attention in the literature. On the one hand, it is due to the role it plays in the context of certain applications in the field of Physics or Engineering and, on the other hand, we especially highlight its relevance within the framework of generic unfoldings of the four-dimensional nilpotent singularity of codimension four.The system exhibits a bifocal equilibrium point for a range of parameter values. The associated two-dimensional invariant manifolds, stable and unstable, fold into the phase space in such a way that they produce intricate patterns. This entangled geometry has previously been called tentacular geometry.We consider a three-dimensional level set containing the bifocal equilibrium point to gain insight into the folding behavior of these invariant manifolds. Our method consists of describing the traces left by invariant manifolds when crossing an invariant cross section by the reversibility map. With this new approach, we provide a better understanding of how tentacular geometry evolves with respect to the parameter.Our techniques enables us to link the tentacular geometry on the cross section with the study of cocooning cascades of homoclinic tangencies. Indeed, we present a general theory to extend to R4 the phenomena related to cocoon bifurcation that classically develop within families of three-dimensional reversible vector fields. On the basis of these results, we conjecture the existence of heteroclinic cycles consisting of two orbits connecting the bifocus with a saddle node periodic orbit.

Action principles for dissipative, non-holonomic Newtonian mechanics

A methodology for deriving dual variational principles for the classical Newtonian mechanics of mass points in the presence of applied forces, interaction forces and constraints, all with a general dependence on particle velocities and positions, is presented. Methods for incorporating constraints are critically assessed. General theory, as well as explicitly worked out variational principles for a dissipative system (due to Lorenz) and a system with anholonomic constraints (due to Pars) are demonstrated. Conditions under which a (family of) dual Hamiltonian flow(s), as well as a constant(s) of motion, may be associated with a conservative or dissipative, and possibly constrained, primal system naturally emerge in this work.

Coarse-graining Hamiltonian systems using WSINDy

Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth ε\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varepsilon$$\\end{document}-dependent Hamiltonian vector field Xε\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$X_\\varepsilon$$\\end{document} possesses an approximate symmetry if the limiting vector field X0=limε→0Xε\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$X_0=\\lim _{\\varepsilon \\rightarrow 0}X_\\varepsilon$$\\end{document} possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the ε>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varepsilon >0$$\\end{document} regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

A Nekhoroshev theorem for some perturbations of the Benjamin-Ono equation with initial data close to finite gap tori

We consider a perturbation of the Benjamin Ono equation with periodic boundary conditions on a segment. We consider the case where the perturbation is Hamiltonian and the corresponding Hamiltonian vector field is analytic as a map from the energy space to itself. Let ϵ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\epsilon $$\\end{document} be the size of the perturbation. We prove that for initial data close in energy norm to an N-gap state of the unperturbed equation all the actions of the Benjamin Ono equation remain O(ϵ12(N+1))\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {O}}(\\epsilon ^{\\frac{1}{2(N+1)}})$$\\end{document} close to their initial value for times exponentially long with ϵ-12(N+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\epsilon ^{-\\frac{1}{2(N+1)}}$$\\end{document}.

Research on shape-controllable generalized multi-cluster Hamiltonian conservative chaotic flow systems and their FPGA implementation

Research on shape-controllable generalized multi-cluster Hamiltonian conservative chaotic flow systems and their FPGA implementation

Nonlinear Stability of Equilibria in Hamiltonian Systems with Multiple Resonances without Interactions

In this work, we advance in the study of the Lyapunov stability and instability of equilibrium solutions of Hamiltonian flows. More precisely, we study the nonlinear stability in the Lyapunov sense of equilibrium solutions in autonomous Hamiltonian systems with n-degrees of freedom, assuming the existence of two resonance vectors k1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extbf{k}^1$$\\end{document} and k2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extbf{k}^2$$\\end{document} without interaction (|k1|≤|k2|\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$|\ extbf{k}^1|\\le |\ extbf{k}^2|$$\\end{document}). We provide conditions to obtain a type of formal stability, called Lie stability. In particular, we need to normalize the Hamiltonian function to any arbitrary order, and our results take into account the sign of the components of the resonance vectors. Subsequently, we guarantee some sufficient conditions to obtain exponential stability in the sense of Nekhoroshev for Lie stable systems. In addition, we give sufficient conditions for the instability in the Lyapunov sense of the full system. For this, it is necessary to normalize the Hamiltonian function to an adequate order, and assuming that the components of at least one resonance vector change of sign.

A new canonical affine bracket formulation of Hamiltonian classical field theories of first order

It has been a long standing question how to extend, in the finite-dimensional setting, the canonical Poisson bracket formulation from classical mechanics to classical field theories, in a completely general, intrinsic, and canonical way. In this paper, we provide an answer to this question by presenting a new completely canonical bracket formulation of Hamiltonian Classical Field Theories of first order on an arbitrary configuration bundle. It is obtained via the construction of the appropriate field-theoretic analogues of the Hamiltonian vector field and of the space of observables, via the introduction of a suitable canonical Lie algebra structure on the space of currents (the observables in field theories). This Lie algebra structure is shown to have a representation on the affine space of Hamiltonian sections, which yields an affine analogue to the Jacobi identity for our bracket. The construction is analogous to the canonical Poisson formulation of Hamiltonian systems although the nature of our formulation is linear-affine and not bilinear as the standard Poisson bracket. This is consistent with the fact that the space of currents and Hamiltonian sections are respectively, linear and affine. Our setting is illustrated with some examples including Continuum Mechanics and Yang–Mills theory.