{
  "generated_at": "2026-10-22T12:59:09Z",
  "papers": [
    {
      "title": "Port-Hamiltonian systems with energy and power ports",
      "summary": "We extend the port-Hamiltonian framework defined with respect to a Lagrangian submanifold and a Dirac structure by augmenting the Lagrangian submanifold with the space of external variables. The new pair of conjugated variables is called energy port. We show that in the most general case, the extension describes constrained Hamiltonian systems whose Hamiltonian function depends on inputs.",
      "url": "http://arxiv.org/abs/2405.01241v2",
      "authors": [
        "Kaja Krhac",
        "Bernhard Maschke",
        "Arjan van der Schaft"
      ],
      "updated": "2026-10-21"
    },
    {
      "title": "Mesh Field Theory: Port-Hamiltonian Formulation of Mesh-Based Physics",
      "summary": "We present Mesh Field Theory (MeshFT) and its neural realization, MeshFT-Net: a structure-preserving framework for mesh-based continuum physics that cleanly separates the physics' topological structure from its metric structure. Imposing minimal physical principles (locality, permutation equivariance, orientation covariance, and energy balance/dissipation inequality), we prove a reduction theorem for mesh-based physics. Under these conditions, the physical dynamics admit a local factorization into a port-Hamiltonian form: the conservative interconnection is fixed uniquely by mesh topology, whereas metric effects enter only through constitutive relations and dissipation. This reduction clarifies what must be fixed and what should be learned, directly informing MeshFT-Net's design. Across evaluations on analytic and realistic datasets, physics-consistency tests, and out-of-distribution validation, MeshFT-Net achieves near-zero energy drift and strong physical fidelity (correct dispersion and momentum conservation) along with robust extrapolation and high data efficiency. By eliminating non-physical degrees of freedom and learning only metric-dependent structure, MeshFT provides a principled inductive bias for stable, faithful, and data-efficient learning-based physical simulation.",
      "url": "http://arxiv.org/abs/2605.00394v4",
      "authors": [
        "Satoshi Noguchi",
        "Yoshinobu Kawahara"
      ],
      "updated": "2026-10-16"
    },
    {
      "title": "Structure-Aware Spatial Kernels for Gaussian Process Distributed Port-Hamiltonian Systems",
      "summary": "Gaussian process distributed port-Hamiltonian systems (GP-dPHS) enable data-driven modeling of distributed dynamical systems while preserving their physical structure. However, standard stationary kernels can be inefficient for spatially heterogeneous systems because they neither exploit known input invariances nor reflect structured spatial dependencies. We address this limitation with a structure-aware kernel design that removes irrelevant state coordinates and decouples state and spatial dependencies through a tensor-product covariance. This modular construction allows qualitative knowledge of spatial heterogeneity to be incorporated by selecting an appropriate spatial kernel without requiring an analytical description of the underlying profile. Numerical experiments on a nonlinear string with periodic, smoothly varying, and semi-piecewise stiffness profiles demonstrate that structural decoupling consistently improves reconstruction accuracy, while profile-matched spatial kernels yield substantially larger gains in force-field and trajectory reconstruction.",
      "url": "http://arxiv.org/abs/2610.22575v1",
      "authors": [
        "Tianyi Yang",
        "Thomas Beckers"
      ],
      "updated": "2026-10-15"
    },
    {
      "title": "Optimal energy conversion -- A port-Hamiltonian approach",
      "summary": "This paper studies optimal energy conversion in port-Hamiltonian systems. We analyze state- and control-constrained finite-horizon optimal control problems that maximize energy conversion efficiencies. Under suitable reachability conditions, we show that the asymptotic optimal efficiencies exist and coincide with the corresponding optimal cyclic efficiencies. For a class of port-Hamiltonian systems describing coupled storage-conversion components, we further derive intrinsic limits on achievable energy conversion and establish a generalized Carnot-type efficiency bound. Our results provide an optimal-control framework connecting asymptotic optimal operation, cyclic operation, and fundamental efficiency limits for port-Hamiltonian systems.",
      "url": "http://arxiv.org/abs/2610.21790v1",
      "authors": [
        "Guanru Pan",
        "Karl Worthmann",
        "Manuel Schaller",
        "Timm Faulwasser"
      ],
      "updated": "2026-10-15"
    },
    {
      "title": "TorchSpin: Differentiable Electron Paramagnetic Resonance Simulation in PyTorch",
      "summary": "Interpreting electron paramagnetic resonance (EPR) spectra commonly requires repeated forward simulation and optimization of spin Hamiltonian parameters. Established software provides mature simulation and fitting methods, but these calculations are not generally available as differentiable components within machine-learning workflows. We present TorchSpin, an open-source PyTorch framework that ports the organization and numerical conventions of EasySpin into a tensor-based computational environment. Representative calculations spanning continuous-wave EPR, ENDOR, pulse EPR, spin dynamics, trajectory simulation, and magnetometry closely reproduced corresponding EasySpin results. In matched benchmarks, TorchSpin CPU execution was faster than EasySpin for six of nine workloads although performance depended strongly on the calculation. We therefore focused on the capabilities enabled by automatic differentiation in PyTorch. In a Cu(II) powder-spectrum demonstration, feedback propagated through TorchSpin improved the quality of a supervised convolutional neural network (CNN) without using parameter labels during refinement. Models trained from randomly initialized CNNs using only noisy spectra and TorchSpin reconstruction loss achieved median similarities of $\\geq$ 0.99. CNNs trained on TorchSpin-generated libraries also transferred without adaptation to independently generated EasySpin spectra with performance comparable to EasySpin-trained controls. Once trained, the CNN produced initial parameter estimates in approximately 1 ms per spectrum while a batch of 512 was processed in approximately 1.3 ms with GPU acceleration. TorchSpin therefore provides both a tested Python implementation of established EPR calculations and a differentiable physical layer for rapid learned inference, spectral refinement, and parameter-label-free training.",
      "url": "http://arxiv.org/abs/2610.18924v1",
      "authors": [
        "Alec H. Follmer"
      ],
      "updated": "2026-10-12"
    },
    {
      "title": "Mixed finite element discretization of intrinsic geometrically exact beams for multiplier-free multibody dynamics",
      "summary": "The Reissner-Simo and Hodges models are two equivalent continuous descriptions of finite-strain beam dynamics. The Reissner-Simo formulation uses displacements and rotations, while the Hodges formulation is intrinsic and avoids both variables. Although equivalent in theory, the two approaches behave differently after discretization and offer distinct numerical advantages. In this work, we develop a structure-preserving discretization of the intrinsic planar beam formulation. Because the intrinsic equations involve linear differential operators, both kinematic and dynamic boundary conditions can be imposed naturally using mixed finite elements. The resulting formulation also enables multibody systems to be assembled without algebraic constraints, avoiding the stiff differential-algebraic equations typically introduced by kinematic constraints. We demonstrate the approach on different examples, also showing that closed kinematic loops can be modeled without algebraic constraints. The resulting interconnected systems retain a port-Hamiltonian structure, with all nonlinearities confined to the interconnection operator. This structure allows exact energy preservation when combined with implicit midpoint time integration. Furthermore the scheme appears to require less Newton iterations compared to existing energy preserving scheme.",
      "url": "http://arxiv.org/abs/2607.20245v2",
      "authors": [
        "Andrea Brugnoli",
        "Philipp L. Kinon",
        "Francesco Sanfedino",
        "Peter Betsch",
        "Olivier A. Bauchau"
      ],
      "updated": "2026-10-12"
    },
    {
      "title": "Well-Posedness and Strong Convergence of Infinite-Dimensional Port-Hamiltonian Primal-Dual Gradient Dynamics",
      "summary": "We establish a functional-analytic framework that justifies the LaSalle convergence argument outlined in [van der Schaft, 2026] on port-Hamiltonian primal-dual gradient dynamics for finite-horizon optimal control of linear systems. In addition to strict convexity of the running cost as in \\cite{vanderSchaft2026PDGJuly}, our analysis requires global Lipschitz continuity of the cost gradient. This ensures that the nonlinear gradient operator is well-defined and Lipschitz on the entire underlying space. We establish global well-posedness of the optimizer dynamics: every initial value determines a unique global mild solution, and the associated solution operators form a contraction semigroup. Similar to [Gernandtand Schaller, 2025], this is proved using tools from maximal monotone operator theory. Assuming existence of an optimal solution, we combine a dissipation inequality with an infinite-dimensional LaSalle invariance principle to prove strong convergence of every trajectory with a relatively compact positive orbit to the unique optimal state, control, and costate. Moreover, we present a sufficient condition for compactness through uniform control regularity and show that convergence holds for all initial states. Finally, we present an example with a strictly but not strongly convex running cost for which every mild solution converges strongly to the unique equilibrium, while no global exponential estimate uniformly over initial values can be obtained.",
      "url": "http://arxiv.org/abs/2610.16263v1",
      "authors": [
        "O. V. Iftime",
        "A. J. van der Schaft"
      ],
      "updated": "2026-10-12"
    },
    {
      "title": "Demonstration Calibrated Port-Hamiltonian Retuning for Manipulation Policies",
      "summary": "Diffusion and VLA policies for manipulation are often deployed through downstream impedance controllers. The stiffness and damping gains of these controllers affect task success, yet are commonly inherited from data collection rather than selected for the deployed policy. Although empirical gain sweeps can improve performance, they require repeated evaluation rollouts. We introduce PHRetune, an offline method that derives controller gains for a frozen policy without evaluation rollouts or gain search. Our approach learns a port-Hamiltonian model from demonstrations to estimate the effort and energy associated with the policy's predicted actions. The policy is applied to recorded demonstration observations, and its predictions are assessed against demonstration-derived effort and energy budgets. From this comparison, we derive a single gain scale in closed form, adjusting the downstream controller while preserving the policy and its action representation. The gains are fixed before evaluation, without requiring a prior manipulator model, task rewards, policy retraining, or additional runtime computation. Across LIBERO suites, PHRetune improves Diffusion Policy success by up to 9.4 percentage points, with the derived gains achieving the highest observed success rates in empirical gain sweeps. On all four real-world manipulation tasks, PHRetuned Diffusion Policy outperforms the nominal policy, alternative gain-tuning methods, and a policy-retraining baseline. The same procedure improves success with SmolVLA and OpenVLA-OFT on every task, while reducing acceleration and jerk for both VLA backbones.",
      "url": "http://arxiv.org/abs/2610.05755v3",
      "authors": [
        "Yulong Yang",
        "Fan Wu",
        "Christine Allen-Blanchette",
        "Amit Chakraborty"
      ],
      "updated": "2026-10-09"
    },
    {
      "title": "Discovering Phase Space Structure in Learned Hamiltonian Systems",
      "summary": "Mechanics-inspired neural dynamical models improve interpretability and physical consistency, but they usually assume a known configuration space, fixed physical parameters, and deterministic prediction from a prescribed initial condition. We introduce Structured Phase Space GAN (SPS-GAN), a conditional generative model that learns structured dynamics from both Cartesian trajectories and video. SPS-GAN combines a port-Hamiltonian backbone with adversarial conditional generation, allowing conservative, dissipative, and forced systems to be represented in a single architecture. To infer latent mechanical structure from raw observations, we use a cyclic-coordinate loss that recovers both the system degrees of freedom and canonical coordinates without access to the true phase space. Experiments on simulated and real-world systems show that SPS-GAN improves trajectory accuracy over supervised dynamics baselines and video consistency over generative baselines, while consistently recovering the correct latent mechanical coordinates across observation modalities and dynamical regimes.",
      "url": "http://arxiv.org/abs/2509.23003v3",
      "authors": [
        "Yulong Yang",
        "Jiayin Liu",
        "Vineet Bansal",
        "Christine Allen-Blanchette"
      ],
      "updated": "2026-10-08"
    },
    {
      "title": "Nonlinear differential-algebraic systems with dissipative port-Hamiltonian structure",
      "summary": "Different representations of nonlinear dissipative Hamiltonian and port-Hamiltonian (pH) differential-algebraic equations (DAE) systems are derived and compared. Using global geometric and algebraic points of view, including singular Lagrange and Dirac structures, as well as Morse families and their skew counterparts, constraints are incorporated in the models. The results are illustrated specifically with constrained mechanical systems using kinematic and geometric constraints. The relation between the different representations is studied and in particular the conversion of Lagrange into Dirac algebraic constraints and vice versa is discussed. To incorporate dissipation also maximally monotone structures are studied.",
      "url": "http://arxiv.org/abs/2610.11455v1",
      "authors": [
        "Volker Mehrmann",
        "Arjan van der Schaft"
      ],
      "updated": "2026-10-08"
    },
    {
      "title": "On Port-Hamiltonian Formulation of Hysteretic Energy Storage Elements: The Backlash Case",
      "summary": "This paper presents a port-Hamiltonian formulation of backlash-driven hysteretic energy storage elements. First, we revisit the passivity property of backlash-driven storage elements by presenting a family of storage functions. We explicitly derive the corresponding available storage and required supply functions in the sense of Willems, and show the interlacing property of the aforementioned family of storage functions sandwiched between the available storage and required supply functions. Second, using the obtained family of storage functions, we present a port-Hamiltonian formulation of hysteretic inductors. In particular, we show how an appropriate Hamiltonian function is defined using the family of storage functions and how the hysteretic elements can be expressed as a port-Hamiltonian system with feedthrough term, where the feedthrough term represents energy dissipation. Correspondingly, we illustrate its applicability in describing an RC circuit (in parallel and in series) containing a backlash inductor.",
      "url": "http://arxiv.org/abs/2603.25211v3",
      "authors": [
        "Jurrien Keulen",
        "Bayu Jayawardhana",
        "Arjan van der Schaft"
      ],
      "updated": "2026-10-06"
    },
    {
      "title": "Port-Hamiltonian Neural Networks for Systems with Multiple Asymptotically Stable Equilibria",
      "summary": "Stable port-Hamiltonian neural networks certify asymptotic stability by construction. Yet, their Hamiltonian is a global Lyapunov function with a single global minimum, so they can represent only dynamic systems with one attractor. We demonstrate that this excludes even simple systems with energy landscapes forming a double well, and we overcome the restriction by parametrising the Hamiltonian as a product of Bregman divergences generated by one input-convex network. We prove that the resulting model is locally Lyapunov stable, that the coexistence of stable equilibria forces additional non-asymptotically-stable equilibria to exist, that all equilibria lie in a bounded region, and under a hyperbolicity assumption that almost-everywhere stability holds. On three systems our approach is able to recover the energy surface characteristics and improve the convergence speed by 1.8$\\times$-8.5$\\times$.",
      "url": "http://arxiv.org/abs/2610.01356v2",
      "authors": [
        "Simon Heilig",
        "Jens Püttschneider",
        "Mohammad Itani",
        "Asja Fischer",
        "Timm Faulwasser"
      ],
      "updated": "2026-10-04"
    },
    {
      "title": "Block numerical ranges for linear relations",
      "summary": "Linear relations are multivalued generalizations of linear operators whose spectral properties are closely linked to the stability of differential-algebraic equations. Consequently, reliable spectral enclosures are of considerable interest. In this note, we extend the established concept of the block numerical range from operators to linear relations by means of a range representation. Unlike the classical numerical range, the resulting block numerical range may be nonconvex and can therefore detect gaps in the spectrum. We establish fundamental properties of this new notion and derive corresponding spectral enclosure results. Finally, we illustrate the theory with a numerical example comparing the block numerical ranges obtained from different constructions. As an application, we obtain an explicit left-half-plane spectral bound for a class of port-Hamiltonian DAEs.",
      "url": "http://arxiv.org/abs/2610.04431v1",
      "authors": [
        "Rabeb Aydi",
        "Hannes Gernandt",
        "Omaima Kchaou",
        "Carsten Trunk",
        "Mitsuru Wilson"
      ],
      "updated": "2026-10-03"
    },
    {
      "title": "Stimulated Parametric Down-Conversion: From Foundational Coherence to Structured Light Applications",
      "summary": "Stimulated parametric down-conversion (StimPDC) is the seeded version of the process that underlies much of photonic quantum information. Replacing the vacuum at one input port with a coherent seed leaves the interaction Hamiltonian, the phase-matching function and the two-photon amplitude unchanged, while raising the output to levels accessible with ordinary classical detectors. This review treats the transverse spatial physics of that process. We show that a single relation, in which the idler amplitude is the pointwise product of the pump field with the complex conjugate of the seed, organizes results usually presented separately: the image and coherence transfer experiments of the 1990s; phase conjugation and its use for wavefront correction; the conservation and conjugation of orbital angular momentum, together with the generation of radial modes and geometric intensity patterns; and the processing of vector vortex beams carrying spatially varying polarization. We then consider StimPDC as an instrument rather than as an object of study, covering stimulated emission tomography, the experimental realization of Klyshko's advanced-wave picture, and the transfer of partial coherence, including the conservation of the twist phase of twisted Gaussian Schell-model beams. We close with the quantum regime, reached when the seed approaches the single-photon level or the idler is used as a herald, and with an assessment of open problems and prospective applications.",
      "url": "http://arxiv.org/abs/2610.03241v1",
      "authors": [
        "M. F. Z. Arruda",
        "M. G. Damaceno",
        "G. H. dos Santos",
        "R. Medeiros de Araújo",
        "N. Rubiano da Silva",
        "S. P. Walborn",
        "P. H. Souto Ribeiro"
      ],
      "updated": "2026-10-02"
    },
    {
      "title": "Nonlinear controlled port-Hamiltonian systems: Existence of (optimal) solutions",
      "summary": "We investigate the existence of solutions of reversible and irreversible port-Hamilto\\-nian systems. To this end, we utilize the associated exergy, a function that is composed of the system's Hamiltonian and entropy, to prove global existence in time for bounded control functions. Then, we rigorously verify our existence conditions for two examples, the gas-piston system and a network of heat exchangers. Last, we explore model predictive control tailored to irreversible port-Hamiltonian systems by means of a numerical case study with a heat exchanger network.",
      "url": "http://arxiv.org/abs/2410.18888v2",
      "authors": [
        "Willem Esterhuizen",
        "Bernhard Maschke",
        "Till Preuster",
        "Manuel Schaller",
        "Karl Worthmann"
      ],
      "updated": "2026-10-02"
    },
    {
      "title": "Identify then Realize: Contrastive Learning of Latent Port-Hamiltonian Dynamics from Partial Observations",
      "summary": "Identifying latent state representations and dynamics is essential when direct modeling in observation space is infeasible, particularly under partial and high-dimensional observations. In such settings, representation learning and physics-aware modeling are inherently coupled. We propose CIPHER, a two-stage identify-then-realize framework for learning latent port-Hamiltonian models of conservative and dissipative systems. First, a contrastive teacher jointly learns an encoder and a neural ODE to obtain predictive state representations from observation histories. Second, a student learns an invertible nonlinear coordinate transformation and port-Hamiltonian dynamics by matching encoded observed futures. We establish sufficient conditions for recovering the physical state up to a diffeomorphism and realizing its dynamics in port-Hamiltonian form. Across ten clean and noisy datasets, CIPHER is competitive with or outperforms state-of-the-art baselines by improving long-horizon forecast on complex systems. Further studies show robustness to over-specified latent dimensions and additional predictive gains from pH parameterization.",
      "url": "http://arxiv.org/abs/2605.16682v2",
      "authors": [
        "Peilun Li",
        "Kaiyuan Tan",
        "Daniel Moyer",
        "Thomas Beckers"
      ],
      "updated": "2026-10-01"
    },
    {
      "title": "Petrov-Galerkin operator inference with application to stability-encouraging identification",
      "summary": "Data-driven model order reduction methods such as operator inference enable the efficient construction of reduced-order models directly from high-dimensional time-domain data. Standard operator inference typically seeks a Galerkin-type reduced model in a prescribed low-dimensional subspace by identifying its reduced operators from projected snapshot data. The resulting inference problem is formulated as a least-squares problem admitting an efficient closed-form solution. However, it is well known from intrusive model order reduction for linear time-invariant systems that Petrov-Galerkin projections can additionally preserve important system properties such as stability and passivity. To overcome the limitations of standard operator inference, we extend the framework for linear time-invariant systems to incorporate Petrov-Galerkin projections and provide explicit error expressions and bounds between the intrusive and nonintrusive reduced operators, thus generalizing results from the literature. We demonstrate the proposed approach in the context of dissipative and port-Hamiltonian systems. Furthermore, we introduce a novel convex optimization formulation that explicitly enforces the port-Hamiltonian structure on the inferred operators. The effectiveness of the proposed methods is demonstrated on several well-established benchmark problems, including the CD player, an atmospheric model, a mass-spring-damper system, and a poroelasticity system.",
      "url": "http://arxiv.org/abs/2610.01295v1",
      "authors": [
        "Johannes Rettberg",
        "Jonas Nicodemus",
        "Harsh Sharma",
        "Boris Kramer",
        "Jörg Fehr",
        "Benjamin Unger"
      ],
      "updated": "2026-10-01"
    },
    {
      "title": "Input-to-state stability of second-order port-Hamiltonian systems under nonlinear dynamic boundary feedback",
      "summary": "The input-to-state stability of a class of infinite-dimensional second-order port-Hamiltonian systems on a one-dimensional spatial domain is analyzed under nonlinear dynamic boundary feedback and boundary disturbances. Using energy methods in the port-Hamiltonian framework, we obtain sufficient conditions that guarantee uniform input-to-state stability of the closed-loop system.",
      "url": "http://arxiv.org/abs/2609.39948v1",
      "authors": [
        "Bouchra Elghazi",
        "Birgit Jacob",
        "Christophe Prieur"
      ],
      "updated": "2026-09-30"
    },
    {
      "title": "Robust forms for passive descriptor systems",
      "summary": "We show how to compute the passivity radius of a given generalized state-space model of a rational transfer function that represents a passive system, both in the continuous-time and discrete-time cases. We then use this to construct robust descriptor realizations. For the continuous-time case, these realizations have a port-Hamiltonian descriptor form. For the discrete-time case a similar result holds for the so-called normalized realizations. The computation of such robust realizations is linked in both these cases to a particular solution of a linear matrix inequality that characterizes the passivity of the transfer function.",
      "url": "http://arxiv.org/abs/2609.39465v1",
      "authors": [
        "Peter Benner",
        "Anshul Prajapati",
        "Paul Van Dooren"
      ],
      "updated": "2026-09-30"
    },
    {
      "title": "Safe-by-Design Learning via Energy-based Neural Networks",
      "summary": "Learning neural-network models of dynamical systems with safety guarantees is a fundamental requirement for their deployment in safety-critical settings. Safety is commonly established by proving the invariance of a desired subset in state-space, ensuring that every trajectory initialized in this subset remains confined to it for all time under admissible inputs. Existing frameworks, however, either rely on computationally expensive post-hoc verification or employ safety-enforcing mechanisms without formal correctness guarantees. In this paper, we introduce a novel neural architecture grounded in energy-based modern Hopfield networks to guarantee safety-by-design while retaining sufficient expressiveness to model complex nonlinear dynamics. Specifically, we integrate modern Hopfield networks with a port-Hamiltonian neural ODE, enabling by design the construction of barrier functions yielding explicit admissible-input sets and quantitative robustness radii. Across several benchmarks, including an 12-dimensional nanodrone model, our framework achieves state-of-the-art performance while producing certified invariant sets that are more robust to external solicitations than comparable existing approaches.",
      "url": "http://arxiv.org/abs/2609.36942v2",
      "authors": [
        "Simone Betteti",
        "Morteza Lahijanian",
        "Luca Laurenti"
      ],
      "updated": "2026-09-30"
    },
    {
      "title": "Port-Hamiltonian Latent Deliberation: Mitigating the Deliberation Drift Cliff in Test-Time Compute Scaling",
      "summary": "Test-time compute scaling has emerged as a cornerstone of advanced machine reasoning, yet performing iterative deliberation directly within continuous latent representation spaces reveals a catastrophic pathology: the Deliberation Drift Cliff. While unconstrained recurrent latent models achieve initial reasoning gains at short horizons (K <= 4), their reasoning collapses when extrapolated to deeper thinking steps (K >= 16), dropping by 22% to 62% across standard logical benchmarks. We resolve the trilemma among expressivity, Lyapunov stability, and computational efficiency in test-time latent reasoning through a 22-round empirical and theoretical investigation. We demonstrate that strictly conservative scalar potential gradient flows suppress long-range drift (cliff 3.40%) but bottleneck peak reasoning accuracy at 32.73%, whereas unconstrained rotational flows achieve high symbolic expressivity (82.33%) but suffer a severe 36.87% drift cliff. To resolve this geometric duality, we establish Port-Hamiltonian Latent Deliberation (PH-LD) and propose the Direct-Gradient Pure-Tensor Helmholtz-Hodge Decomposition (DG-HHD). DG-HHD parameterizes the attracting flow as a tangent projection tensor network while orthogonally decoupling non-zero circulation (Hodge machine error 1.65e-17, contraction error 5.55e-17), eliminating runtime autograd dependencies to achieve 1.84x vector field and 2.09x RK45 rollout speedups. In a 15-arm symmetrical Pareto benchmark, DG-HHD achieves 58.67% peak accuracy (+25.94% absolute gain over conservative HHD) and retains 35.27% at K=32. Transferred to small language model (SLM) multi-hop causal reasoning, DG-HHD delivers monotonic compute scaling (49.33% to 51.56%) and suppresses out-of-distribution drift (cliff -0.66%). All 30 Level 0 deterministic invariants are certified.",
      "url": "http://arxiv.org/abs/2609.37351v1",
      "authors": [
        "Zeyu Jia"
      ],
      "updated": "2026-09-29"
    },
    {
      "title": "Mathematical modeling, simulation and optimization via (port-Hamiltonian) differential-algebraic equations",
      "summary": "Differential-algebraic equations (DAEs) present a ubiquitous model class in the automated modeling, simulation, control and optimization of complex finite or infinite dimensional dynamical systems, in particular in the construction of digital twins. When combined with adequate space and time discretization methods they present the work-horse of developments in many engineering domains. However, this approach comes with certain challenges for the numerical methods and in particular in the coupling of systems from different physical domains with different scales or regularity requirements. We survey the modeling with DAEs and address how to overcome some of the challenges via the framework of port-Hamiltonian differential-algebraic systems (pHDAEs) and their simulation and optimization with space-time-model adaptive methods.",
      "url": "http://arxiv.org/abs/2609.33954v1",
      "authors": [
        "Volker Mehrmann"
      ],
      "updated": "2026-09-27"
    },
    {
      "title": "Hamiltonian JEPA: Action-Conditioned World Models with an Inherited Control State",
      "summary": "Planning from pixels needs more than a latent space that is stable and predictable. The state the planner scores must also be organized by how actions move the system. Joint-embedding predictive architectures (JEPAs) avoid pixel reconstruction by predicting future representations, but existing action-conditioned JEPAs ask one embedding to serve both perception and control. We introduce H-JEPA, which separates the two. A wide perceptual code is regularized toward a well-scaled isotropic geometry with a Bures-Wasserstein prior, and a fixed orthonormal slice of that code is the control state, which inherits the code's covariance without any objective of its own. The state evolves under phase-conditioned dissipative port-Hamiltonian dynamics whose input port has orthonormal columns. Port-inverse consistency (PIC) reads the executed action back through the transpose of that port. We show that this readout is exactly the rollout error projected onto the port directions, so PIC is a parameter-free reweighting of prediction error and not an auxiliary action decoder. Untying the readout from the port breaks this identity and loses half of the gain. H-JEPA matches or exceeds reconstruction-free baselines, including the action-decoding Delta-JEPA, on four pixel-based control benchmarks after at most $10$ training epochs, and its largest gain is on OGB-Cube ($91.9$ against $79.3$ percent). Ablations on PushT and OGB-Cube separate the contributions of the structured predictor, PIC, the prediction horizon, the state rank, and the anti-collapse prior.",
      "url": "http://arxiv.org/abs/2609.33497v1",
      "authors": [
        "Tamim Zoabi",
        "Ameen Ali",
        "Lior Wolf"
      ],
      "updated": "2026-09-27"
    },
    {
      "title": "SIPHy: Sparse identification of port-Hamiltonian systems from noisy data",
      "summary": "We propose sparse identification of port-Hamiltonian systems (SIPHy), enabling structure-preserving symbolic regression from noisy trajectory observations. The method applies to port-Hamiltonian systems, which provide a general framework for describing dynamical systems in terms of energy exchange, dissipation and control. Our algorithm can jointly identify the Hamiltonian as well as the dissipation and input matrices. Furthermore, we introduce Hamiltonian flow splines to better approximate derivatives of trajectory data corrupted by noise or with missing time points, a major challenge of system identification for differential equations. This method assembles flows of piecewise polynomial Hamiltonians to produce a smooth, differentiable trajectory necessary for sparse regression. Combining flow splines with SIPHy yields interpretable, physically grounded models from noisy and incomplete trajectory data. Because the Hamiltonian, dissipation, and input matrices are identified as separate components, the resulting models can be simulated under control inputs and dissipation regimes never observed during training, thus improving the robustness and generalization capabilities of model discovery methods.",
      "url": "http://arxiv.org/abs/2609.20399v2",
      "authors": [
        "Håkon Noren Myhr",
        "Sølve Eidnes",
        "J. Nathan Kutz"
      ],
      "updated": "2026-09-23"
    },
    {
      "title": "Energy-optimal predictive control of discrete-time port-Hamiltonian systems: Closed-loop practical stability",
      "summary": "We study a dissipativity-based model predictive control (MPC) scheme for energy-optimal constrained output stabilization of a discrete-time nonlinear SISO port-Hamiltonian system described by difference and differential representation. For the optimal control problem to be solved in each MPC step, we establish a measure turnpike behavior. The turnpike set in our work is obtained from the underlying port-Hamiltonian structure without assuming existence of steady-states or periodic orbits. In particular, unlike existing results which typically establish the turnpike property w.r.t. a controlled forward invariant set such as optimal steady-states or optimal periodic orbits, we do not require the turnpike set to be forward control invariant. For MPC, we establish recursive feasibility and show practical stability to the turnpike set w.r.t. the MPC closed loop leveraging both dissipativity and the port-Hamiltonian structure. Finally, we demonstrate our results using two numerical examples.",
      "url": "http://arxiv.org/abs/2609.25920v1",
      "authors": [
        "Vaibhav Kumar Singh",
        "Manuel Schaller",
        "Timm Faulwasser",
        "Karl Worthmann"
      ],
      "updated": "2026-09-22"
    }
  ]
}