Hamiltonian Structure Research Articles

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.

Hamiltonian formulation of X-point collapse in an extended magnetohydrodynamics framework

The study of X-point collapse in magnetic reconnection has witnessed extensive research in the context of space and laboratory plasmas. In this paper, a recently derived mathematical formulation of X-point collapse applicable in the regime of extended magnetohydrodynamics is shown to possess a noncanonical Hamiltonian structure composed of five dynamical variables inherited from its parent model. The Hamiltonian and the noncanonical Poisson brackets are both derived, and the latter is shown to obey the requisite properties of antisymmetry and the Jacobi identity (an explicit proof of the latter is provided). In addition, the governing equations for the Casimir invariants are presented, and one such solution is furnished. The above features can be harnessed and expanded in future work, such as developing structure-preserving integrators for this dynamical system.

Hamiltonian Birkhoff Normal Form for Gravity-Capillary Water Waves with Constant Vorticity: Almost Global Existence

We prove an almost global existence result for space periodic solutions of the 1D gravity-capillary water waves equations with constant vorticity. The result holds for any value of gravity, vorticity and depth, a full measure set of surface tensions, and any small and smooth enough initial datum. The proof demands a novel approach—that we call paradifferential Hamiltonian Birkhoff normal form for quasi-linear PDEs—in presence of resonant wave interactions: the normal form is not integrable but it preserves the Sobolev norms thanks to its Hamiltonian nature. A major difficulty is that paradifferential calculus used to prove local well posedness (as the celebrated Alinhac good unknown) breaks the Hamiltonian structure. A major achievement of this paper is to correct (possibly) unbounded paradifferential transformations to symplectic maps, up to an arbitrary degree of homogeneity. Thanks to a deep cancellation, our symplectic correctors are smoothing perturbations of the identity. Thus we are able to preserve both the paradifferential structure and the Hamiltonian nature of the equations. Such Darboux procedure is written in an abstract functional setting applicable also in other contexts.

Levinson–Smith Dissipative Equations and Geometry of GENERIC Formalism and Contact Hamiltonian Mechanics

We apply Jacobi’s Last Multiplier theory to construct the non-standard Lagrangian and Hamiltonian structures for the Levinson–Smith equations satisfying the Chiellini integrability condition. Then after a brief exposition of the contact geometry, we explore its connection with the non-standard Hamiltonian structures. We present the formulation of the Levinson–Smith equation in terms of General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) method and also study the gradient-type flow. We give a geometric formulation of GENERIC and apply this to general Levinson–Smith equations.

Isomonodromic and isospectral deformations of meromorphic connections: the sl2(C) case

Abstract We consider non-twisted meromorphic connections in sl 2 ( C ) and the associated symplectic Hamiltonian structure. In particular, we provide explicit expressions of the Lax pair in the geometric gauge supplementing results of Marchal et al (2022 arXiv:2212.04833) where explicit formulas have been obtained in the oper gauge. Expressing the geometric Lax matrices requires the introduction of specific Darboux coordinates for which we provide the explicit Hamiltonian evolutions. These expressions allow to build bridges between the isomonodromic deformations and the isospectral ones. More specifically, we propose an explicit change of Darboux coordinates to obtain isospectral coordinates for which Hamiltonians match the spectral invariants. This result solves the issue left opened in Bertola et al (2023 J. Math. Phys. 64 083502) in the case of sl 2 ( C ) .

Identification of Hamiltonian systems using neural networks and first integrals approaches

This research introduces a class of non-parametric identifiers based on differential neural networks represented by Hamiltonian dynamics. The structure of the identifier corresponds to the form of a canonical Hamiltonian system that uses the evolution of generalized coordinates and momentums. The learning laws of the identifier come from applying the first integrals approach, which justifies the design of an exact identifier considering the time invariance of the Hamiltonian, with a finite number of activation functions in the identifier structure. The first integrals approach derives several learning laws for the proposed class of identifiers. The learning laws design uses the estimated derivative of the generalized momentum assessed via a super-twisting differentiator with multiple inputs and outputs. All proposed laws require the solution of differential continuous-time Riccati equations and nonlinear differential equations for the learning laws, which depend on the identification error and state constraints. The developed identifier was evaluated compared to an identifier that did not consider the Hamiltonian constraint using first integrals. This comparison included a numerical evaluation of the identifier considering its application to a classical Hamiltonian system associated with the Kepler dynamics representing satellite orbital evolution. This evaluation confirmed that the identification results were improved with the proposed learning laws regarding the class of Hamiltonian structures, and the quality indicators based on the mean square error were several times lower.

A reduced ideal MHD system for nonlinear magnetic field turbulence in plasmas with approximate flux surfaces

This paper studies the nonlinear evolution of magnetic field turbulence in proximity of steady ideal Magnetohydrodynamics (MHD) configurations characterized by a small electric current, a small plasma flow, and approximate flux surfaces, a physical setting that is relevant for plasma confinement in stellarators. The aim is to gather insight on magnetic field dynamics, to elucidate accessibility and stability of three-dimensional MHD equilibria, as well as to formulate practical methods to compute them. Starting from the ideal MHD equations, a reduced dynamical system of two coupled nonlinear partial differential equations for the flux function and the angle variable associated with the Clebsch representation of the magnetic field is obtained. It is shown that under suitable boundary and gauge conditions such reduced system preserves magnetic energy, magnetic helicity, and total magnetic flux. The noncanonical Hamiltonian structure of the reduced system is identified, and used to show the nonlinear stability of steady solutions against perturbations involving only one Clebsch potential. The Hamiltonian structure is also applied to construct a dissipative dynamical system through the method of double brackets. This dissipative system enables the computation of MHD equilibria by minimizing energy until a critical point of the Hamiltonian is reached. Finally, an iterative scheme based on the alternate solution of the two steady equations in the reduced system is proposed as a further method to compute MHD equilibria. A theorem is proven which states that the iterative scheme converges to a nontrivial MHD equilbrium as long as solutions exist at each step of the iteration.

Stability of a point charge for the repulsive Vlasov–Poisson system

We consider solutions of the repulsive Vlasov–Poisson system which are a combination of a point charge and a small gas, i.e., measures of the form \delta_{(\mathcal{X}(t),\mathcal{V}(t))}+\mu^{2}d\mathbf{x}d\mathbf{v} for some (\mathcal{X}, \mathcal{V})\colon \mathbb{R}\to\mathbb{R}^{6} and a small gas distribution \mu\colon \mathbb{R}\to L^{2}_{\mathbf{x},\mathbf{v}} , and study asymptotic dynamics in the associated initial value problem. If initially suitable moments on \mu_{0}=\mu(t=0) are small, we obtain a global solution of the above form, and the electric field generated by the gas distribution \mu decays at an almost optimal rate. Assuming in addition boundedness of suitable derivatives of \mu_{0} , the electric field decays at an optimal rate, and we derive modified scattering dynamics for the motion of the point charge and the gas distribution. Our proof makes crucial use of the Hamiltonian structure. The linearized system is transport by the Kepler ODE, which we integrate exactly through an asymptotic action-angle transformation. Thanks to a precise understanding of the associated kinematics, moment and derivative control is achieved via a bootstrap analysis that relies on the decay of the electric field associated to \mu . The asymptotic behavior can then be deduced from the properties of Poisson brackets in asymptotic action coordinates.

On the stability of fully nonlinear hydraulic-fall solutions to the forced water wave problem

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

Interconnection of irreversible port Hamiltonian systems

This paper shows how the interconnection of two controlled Irreversible port Hamiltonian Systems has to be state and co-state modulated in order to ensure the closed-loop Irreversible port Hamiltonian structure, satisfying the first and second laws of Thermodynamics. It proposes a precise parametrization of this modulation from the open-loop systems structures in order to guarantee the consistency of the closed loop energy and entropy balance equations. The results are illustrated by means of the examples of a heat-exchanger, a gas-piston system and a chemical reaction.

An Integrated Integrable Hierarchy Arising from a Broadened Ablowitz–Kaup–Newell–Segur Scenario

This study introduces a 4×4 matrix eigenvalue problem and develops an integrable hierarchy with a bi-Hamiltonian structure. Integrability is ensured by the zero-curvature condition, while the Hamiltonian structure is supported by the trace identity. Explicit derivations yield second-order and third-order integrable equations, illustrating the integrable hierarchy.

On Properties of Adjoint Systems for Evolutionary PDEs

We investigate the geometric structure of adjoint systems associated with evolutionary partial differential equations at the fully continuous, semi-discrete, and fully discrete levels and the relations between these levels. We show that the adjoint system associated with an evolutionary partial differential equation has an infinite-dimensional Hamiltonian structure, which is useful for connecting the fully continuous, semi-discrete, and fully discrete levels. We subsequently address the question of discretize-then-optimize versus optimize-then-discrete for both semi-discretization and time integration, by characterizing the commutativity of discretize-then-optimize methods versus optimize-then-discretize methods uniquely in terms of an adjoint-variational quadratic conservation law. For Galerkin semi-discretizations and one-step time integration methods in particular, we explicitly construct these commuting methods by using structure-preserving discretization techniques.

Large-signal stability analysis of passivity-based droop control of islanded microgrids

In this paper, the large-signal stability properties of a system formed of a single-phase AC Microgrid equipped with a two-level control scheme are analyzed. The aim is to show that it is possible to formally characterize the operative capacities of the controlled system when the nonlinear model that describes its dynamical behavior is considered. Thus, the presented analysis contributes to the formulation of the widely demanded methodology to analyze the large-signal stability properties of Microgrids without invoking small-signal arguments. The analyzed model includes the dynamic of the power converters associated with each distributed energy resource, an inner control law, and a power-sharing mechanism. To better illustrate the scope of the contribution, it is assumed that the Microgrid operates in islanded mode. The inner control schemes are designed using passivity-based ideas, while the power-sharing stage includes the most popular and exhaustively implemented Droop control. The large-signal stability analysis of the Microgrid is carried out by exploiting the Port-controlled Hamiltonian structure of the network and taking advantage of the stability and passivity properties of this class of dynamical systems together with the Input-to-State stability properties of the Droop scheme.

Pseudo grid-based physics-informed convolutional-recurrent network solving the integrable nonlinear lattice equations

Traditional discrete learning methods involve discretizing continuous equations using difference schemes, necessitating considerations of stability and convergence. Integrable nonlinear lattice equations possess a profound mathematical structure that enables them to revert to continuous integrable equations in the continuous limit, particularly retaining integrable properties such as conservation laws, Hamiltonian structure, and multiple soliton solutions. The pseudo grid-based physics-informed convolutional-recurrent network (PG-PhyCRNet) is proposed to investigate the localized wave solutions of integrable lattice equations, which significantly enhances the model’s extrapolation capability to lattice points beyond the temporal domain. We conduct a comparative analysis of PG-PhyCRNet with and without pseudo grid by investigating the multi-soliton solutions and rational solitons of the Toda lattice and self-dual network equation. The results indicate that the PG-PhyCRNet excels in capturing long-term evolution and enhances the model’s extrapolation capability for solitons, particularly those with steep waveforms and high wave speeds. Finally, the robustness of the PG-PhyCRNet method and its effect on the prediction of solutions in different scenarios are confirmed through repeated experiments involving pseudo grid partitioning.

Constraint structure of the generalized Proca model in the Lagrangian formalism

We introduce a novel Lagrangian approach to analyze the dynamical structure of the generalized Proca theory (GP). This method incorporates the A-Z constraint structure of the model within the Lagrangian formalism, providing a precise determination of the degrees of freedom. Furthermore, we present a comprehensive analysis of the model's Hamiltonian constraint structure.

Hamiltonian shocks

AbstractWave propagation in the form of fronts or kinks, a common occurrence in a wide range of physical phenomena, is studied in the context of models defined by their Hamiltonian structure. Motivated, for dispersive wave evolution equations such as a strongly nonlinear model of two‐layer internal waves in the Boussinesq limit, by the symmetric properties of a class of front‐propagating solutions, known as conjugate states or solibores, a generalized formulation based purely on the dispersionless reduction of a system is introduced, and a class of undercompressive shock solutions, here referred to as “Hamiltonian shocks,” is defined. This analysis determines whether a Hamiltonian shock, representing locally a kink for the parent dispersive equations, will interact with a sufficiently smooth background wave without inducing loss of regularity, which would take the form of a classical dispersive shock for the parent equations. This property is also related to an infinitude of conservation laws, drawing a parallel to the case of completely integrable systems.

Complex Fluid Models of Mixed Quantum–Classical Dynamics

Several methods in nonadiabatic molecular dynamics are based on Madelung’s hydrodynamic description of nuclear motion, while the electronic component is treated as a finite-dimensional quantum system. In this context, the quantum potential leads to severe computational challenges and one often seeks to neglect its contribution, thereby approximating nuclear motion as classical. The resulting model couples classical hydrodynamics for the nuclei to the quantum motion of the electronic component, leading to the structure of a complex fluid system. This type of mixed quantum–classical fluid models has also appeared in solvation dynamics to describe the coupling between liquid solvents and the quantum solute molecule. While these approaches represent a promising direction, their mathematical structure requires a certain care. In some cases, challenging higher-order gradients make these equations hardly tractable. In other cases, these models are based on phase-space formulations that suffer from well-known consistency issues. Here, we present a new complex fluid system that resolves these difficulties. Unlike common approaches, the current system is obtained by applying the fluid closure at the level of the action principle of the original phase-space model. As a result, the system inherits a Hamiltonian structure and retains energy/momentum balance. After discussing some of its structural properties and dynamical invariants, we illustrate the model in the case of pure-dephasing dynamics. We conclude by presenting some invariant planar models.

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.

Integrable deformations of Rikitake systems, Lie bialgebras and bi-Hamiltonian structures

Integrable deformations of a class of Rikitake dynamical systems are constructed by deforming their underlying Lie–Poisson Hamiltonian structures, which are considered linearizations of Poisson–Lie structures on certain (dual) Lie groups. By taking into account that there exists a one-to one correspondence between Poisson–Lie groups and Lie bialgebra structures, a number of deformed Poisson coalgebras can be obtained, which allow the construction of integrable deformations of coupled Rikitake systems. Moreover, the integrals of the motion for these coupled systems can be explicitly obtained by means of the deformed coproduct map. The same procedure can be also applied when the initial system is bi-Hamiltonian with respect to two different Lie–Poisson algebras. In this case, to preserve a bi-Hamiltonian structure under deformation, a common Lie bialgebra structure for the two Lie–Poisson structures has to be found. Coupled dynamical systems arising from this bi-Hamiltonian deformation scheme are also presented, and the use of collective ‘cluster variables’, turns out to be enlightening in order to analyse their dynamical behaviour. As a general feature, the approach here presented provides a novel connection between Lie bialgebras and integrable dynamical systems.

Nonlinear Reduced Order Modelling for Korteweg-de Vries Equation

Efficient computational techniques that maintain the accuracy and invariant preservation property of the Korteweg-de Vries (KdV) equations have been studied by a wide range of researchers. In this paper, we introduce a reduced order model technique utilizing kernel principle component analysis (KPCA), a nonlinear version of the classical principle component analysis, in a non-intrusive way. The KPCA is applied to the data matrix, which is formed by the discrete solution vectors of KdV equation. In order to obtain the discrete solutions, the finite differences are used for spatial discretization, and linearly implicit Kahan's method for the temporal one. The back-mapping from the reduced dimensional space, is handled by a non-iterative formula based on the idea of multidimensional scaling (MDS) method. Through KPCA, we illustrate that the reduced order approximations conserve the invariants, i.e., Hamiltonian, momentum and mass structure of the KdV equation. The accuracy of reduced solutions, conservation of invariants, and computational speed enhancements facilitated by classical (linear) PCA and KPCA are exemplified through one-dimensional KdV equation.