pith. sign in

arxiv: 2604.18706 · v1 · submitted 2026-04-20 · 🧮 math.DS

Transfer Operators for Stochastic Hybrid Systems on Manifolds with Guard-Induced Resets

Tejaswi K.C. , William A. Clark , Taeyoung Lee This is my paper

Pith reviewed 2026-05-10 03:44 UTC · model grok-4.3

classification 🧮 math.DS
keywords transfer operatorsstochastic hybrid systemsKoopman operatorFrobenius-Perron operatormanifoldsguard-induced resetsuncertainty propagationfinite volume method
0
0 comments X

The pith

A transfer operator framework unifies uncertainty propagation for stochastic hybrid systems on manifolds with guard-induced resets.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a transfer operator framework that treats both the Koopman operator for observables and the Frobenius-Perron operator for densities in stochastic hybrid systems whose state space is a differentiable manifold. By using their duality, the authors obtain adjoint generators that evolve these quantities according to the backward and forward Kolmogorov equations while incorporating jumps at guards and the associated resets. The entire construction is kept global and intrinsic to the manifold geometry rather than relying on local coordinates. A finite-volume discretization is also introduced that conserves total probability mass and correctly accounts for probability fluxes across guards. If this construction holds, uncertainty in hybrid systems can be tracked without introducing coordinate artifacts or mass leaks at discrete transitions.

Core claim

The central claim is that the duality between observables and probability densities extends to stochastic hybrid systems on manifolds by defining adjoint generators for the backward and forward Kolmogorov equations that remain valid across guard-induced resets, yielding a unified evolution law together with a finite-volume scheme that preserves total probability while capturing reset-induced transfers.

What carries the argument

The duality between observables and probability densities realized through adjoint generators of the backward and forward Kolmogorov equations, extended intrinsically across guard-induced resets on manifolds.

If this is right

  • Observables and probability densities evolve consistently under a single pair of adjoint operators that incorporate both continuous stochastic flow and discrete resets.
  • Total probability mass is preserved exactly by the finite-volume scheme even when probability crosses guards.
  • Fluxes induced by resets are captured directly without requiring separate jump maps or coordinate changes.
  • The formulation remains valid on any differentiable manifold and does not depend on a particular choice of local coordinates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Coordinate-free implementations on manifolds such as spheres or tori become feasible without switching between charts at each step.
  • The same duality structure may allow direct transfer of existing numerical transfer-operator techniques from non-hybrid systems to the hybrid setting.

Load-bearing premise

The duality between observables and probability densities extends globally and intrinsically to manifolds with guard-induced resets without additional regularity conditions or coordinate singularities.

What would settle it

A concrete numerical example on a manifold such as the circle or sphere in which the computed density after a guard crossing and reset fails to integrate to one or yields different results when recomputed in an overlapping coordinate chart.

Figures

Figures reproduced from arXiv: 2604.18706 by Taeyoung Lee, Tejaswi K.C., William A. Clark.

Figure 1
Figure 1. Figure 1: Example 1: the state space is X = [a, b] (shaded in blue), and both endpoints are reflecting. The vector field X(x) = −γ(x − c) for γ > 0 and c > b points to the right toward x = c. and prescribe reflecting boundary conditions on Γ. Thus, in the terminology of Assumption 2, the entire boundary is reflecting: Γ ref = ∂X, Γ abs = ∅. No reset map is involved in this example. 0 0.5 1 1.5 2 0 0.5 1 1.5 2 2.5 x … view at source ↗
Figure 2
Figure 2. Figure 2: Example 1: numerical simulation of the Frobenius–Perron operator. [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Example 1: numerical simulation of the Koopman operator. [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Example 2: the trajectory is reflected at [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Example 2: numerical simulation of the Frobenius–Perron operator [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Example 2: numerical simulation of the Koopman operator for mixed [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The reset condition (22) implies continuity of the Koopman observable across the jump: u(t, b) = u(t, a) [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 7
Figure 7. Figure 7: Example 3: trajectories are reflected at [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Example 3: numerical simulation of the Frobenius–Perron operator [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Example 3: numerical simulation of the Koopman operator for [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Decomposition of the torus (solid), with periodic identification ++ [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Evolution of the density (τ = T /50 = 0.06) and the diffusion coefficient is set to H = 0.5. The initial densities in the two modes are chosen to be the Gaussian in ψ and uniform in θ, namely, v1(0, θ, ψ) = N (µ1, σ2 1 )·1θ, v2(0, θ, ψ) = N (µ2, σ2 2 )·1θ, where 1θ denotes the uniform density on S 1 , and µ1 = 3π 4 , σ1 = 0.2, µ2 = 5π 4 , σ2 = 0.2. The grid resolution is set to Nψ = 100 and Nθ = 100, and … view at source ↗
Figure 12
Figure 12. Figure 12: Evolution of the total probability mass the guard (red circle) and subsequently reset to ψ = 0 (blue circle) at the opposite phase of θ, as shown in [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Verification of the computed density (left column) against Monte [PITH_FULL_IMAGE:figures/full_fig_p013_13.png] view at source ↗
read the original abstract

This paper develops a transfer operator framework for stochastic hybrid systems with guard-induced resets, encompassing both the Koopman and Frobenius--Perron operators. Exploiting their duality, we derive a unified formulation in which observables and probability densities evolve under adjoint generators corresponding to the backward and forward Kolmogorov equations. The formulation is developed in a global and intrinsic manner on differentiable manifolds, ensuring consistency with the underlying geometric structure of the state space. In addition, we propose a finite volume computational scheme on manifolds that preserves total probability mass while accurately capturing fluxes across guards and reset-induced transfers. The proposed framework provides a unified and geometrically consistent approach to uncertainty propagation in stochastic hybrid systems, bridging continuous stochastic dynamics and hybrid transitions within a transfer operator perspective.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper develops a transfer operator framework for stochastic hybrid systems on manifolds with guard-induced resets. It encompasses both the Koopman and Frobenius-Perron operators, exploits their duality to derive adjoint generators for the backward and forward Kolmogorov equations, and presents the formulation intrinsically on differentiable manifolds. A finite-volume computational scheme is proposed that preserves total probability mass while capturing fluxes across guards and reset-induced transfers.

Significance. If the derivations hold, the framework unifies uncertainty propagation for stochastic hybrid systems by bridging continuous stochastic dynamics and hybrid transitions in a geometrically consistent, manifold-intrinsic manner. This could be significant for applications in control theory, robotics, and biological modeling, where state spaces are manifolds and resets occur at guards. The emphasis on duality and mass preservation is a potential strength, though verification details would strengthen the contribution.

major comments (1)
  1. The abstract asserts a global and intrinsic formulation on manifolds without coordinate singularities, but the weakest assumption (duality extending without additional regularity) requires explicit treatment of function spaces and conditions on the vector fields, diffusion coefficients, and reset maps to ensure the adjoint generators are well-defined across guard surfaces.
minor comments (1)
  1. The abstract could be strengthened by briefly indicating the function spaces (e.g., C^1 or L^1) in which the operators act and whether any numerical validation or example is included in the manuscript.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough reading, positive evaluation, and constructive suggestion. The comment identifies a point where greater explicitness will strengthen the rigor of the presentation. We address it below and have incorporated a targeted clarification in the revised manuscript.

read point-by-point responses

  1. Referee: The abstract asserts a global and intrinsic formulation on manifolds without coordinate singularities, but the weakest assumption (duality extending without additional regularity) requires explicit treatment of function spaces and conditions on the vector fields, diffusion coefficients, and reset maps to ensure the adjoint generators are well-defined across guard surfaces.

    Authors: We agree that an explicit statement of the functional-analytic setting is necessary for a fully rigorous justification of the adjoint generators. The manuscript develops the transfer operators intrinsically via the Riemannian structure and the duality between the Koopman and Frobenius–Perron operators, but the precise function spaces and regularity hypotheses were stated only implicitly through the geometric assumptions. In the revised version we have inserted a new paragraph at the end of Section 2 (Preliminaries) that specifies: (i) the underlying space is the Sobolev space W^{1,2}(M) for observables and the dual space of probability measures with finite first moments for densities; (ii) the drift vector fields are C^2 with globally bounded derivatives, the diffusion coefficients are C^2 and uniformly elliptic; and (iii) the reset maps are C^1 diffeomorphisms that map guard surfaces to target surfaces while preserving the manifold structure. These conditions guarantee that the generators remain well-defined in the weak sense across the guards and that the duality relation holds without coordinate charts. We believe this addition removes any ambiguity while preserving the intrinsic, coordinate-free character of the framework. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper develops a transfer operator framework by extending the standard duality between Koopman and Frobenius-Perron operators to stochastic hybrid systems on manifolds with guard-induced resets. It derives adjoint generators corresponding to the backward and forward Kolmogorov equations in an intrinsic geometric manner and proposes a mass-preserving finite-volume scheme. No derivation step reduces by construction to a fitted parameter, self-defined quantity, or load-bearing self-citation; the central claims rest on established operator duality and geometric consistency without renaming known results or smuggling ansatzes via prior author work. The formulation is self-contained as a theoretical proposal against external benchmarks in dynamical systems and hybrid systems literature.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The abstract invokes standard differential geometry on manifolds and the duality of transfer operators; no free parameters, new entities, or ad-hoc axioms are explicitly introduced in the summary.

axioms (2)
  • domain assumption State space is a differentiable manifold with guard sets and reset maps that preserve the geometric structure
    The framework is developed globally and intrinsically on manifolds.
  • domain assumption Duality between Koopman and Frobenius-Perron operators extends to the hybrid setting with resets
    Used to derive the unified adjoint generators corresponding to backward and forward Kolmogorov equations.

pith-pipeline@v0.9.0 · 5423 in / 1323 out tokens · 23084 ms · 2026-05-10T03:44:11.036348+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

  1. [1]

    Goebel, R

    R. Goebel, R. G. Sanfelice, and A. R. Teel,Hybrid Dynamical Systems: Modeling, Stability, and Robustness. Princeton University Press, 2012

    work page 2012

  2. [2]

    Stochastic hybrid systems: Application to communi- cation networks,

    J. P. Hespanha, “Stochastic hybrid systems: Application to communi- cation networks,” inHybrid Systems: Computation and Control (HSCC 2004), ser. Lecture Notes in Computer Science, vol. 2993, 2004, pp. 387–401

    work page 2004

  3. [3]

    Towards a theory of stochastic hybrid systems,

    J. Hu, J. Lygeros, and S. Sastry, “Towards a theory of stochastic hybrid systems,” inInternational Workshop on Hybrid Systems: Computation and Control. Springer, 2000, pp. 160–173

    work page 2000

  4. [4]

    Stability analysis for stochastic hybrid systems: A survey,

    A. R. Teel, A. Subbaraman, and A. Sferlazza, “Stability analysis for stochastic hybrid systems: A survey,”Automatica, vol. 50, pp. 2435– 2456, 2014

    work page 2014

  5. [5]

    Review of stochastic hybrid systems with applications in biological systems modeling and analysis,

    X. Li, O. Omotere, L. Qian, and E. R. Dougherty, “Review of stochastic hybrid systems with applications in biological systems modeling and analysis,”EURASIP Journal on Bioinformatics and Systems Biology, vol. 2017, no. 1, p. 8, 2017

    work page 2017

  6. [6]

    Stochastic hybrid systems in cellular neuroscience,

    P. C. Bressloff and J. N. Maclaurin, “Stochastic hybrid systems in cellular neuroscience,”The Journal of Mathematical Neuroscience, vol. 8, no. 1, p. 12, 2018

    work page 2018

  7. [7]

    Spectrum of the Koopman operator, spectral expansions in functional spaces, and state-space geometry

    I. Mezic, “Spectrum of the Koopman operator, spectral expansions in functional spaces, and state-space geometry.”J. Nonlinear Sci., vol. 30, no. 5, pp. 2091–2145, 2020

    work page 2091

  8. [8]

    Koopman and Perron– Frobenius operators on reproducing kernel Banach spaces,

    M. Ikeda, I. Ishikawa, and C. Schlosser, “Koopman and Perron– Frobenius operators on reproducing kernel Banach spaces,”Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 32, no. 12, 2022

    work page 2022

  9. [9]

    Risken,The Fokker-Planck Equation: Methods of Solution and Applications, 2nd ed

    H. Risken,The Fokker-Planck Equation: Methods of Solution and Applications, 2nd ed. Springer, 1996

    work page 1996

  10. [10]

    G. A. Pavliotis,Stochastic Processes and Applications: Diffusion Pro- cesses, the Fokker-Planck and Langevin Equations. Springer, 2014

    work page 2014

  11. [11]

    A study of the long-term behavior of hybrid systems with symmetries via reduction and the Frobenius–Perron operator,

    M. Oprea, A. Shaw, R. Huq, K. Iwasaki, D. Kassabova, and W. Clark, “A study of the long-term behavior of hybrid systems with symmetries via reduction and the Frobenius–Perron operator,”SIAM Journal on Applied Dynamical Systems, vol. 23, no. 4, pp. 2899–2938, 2024

    work page 2024

  12. [12]

    Learning Koopman representations for hybrid systems,

    C. Bakker, A. Bhattacharya, S. Chatterjee, C. J. Perkins, and M. R. Oster, “Learning Koopman representations for hybrid systems,” 2020. [Online]. Available: https://arxiv.org/abs/2006.12427

    work page arXiv 2020

  13. [13]

    A model for stochastic hybrid systems with application to communication networks,

    J. P. Hespanha, “A model for stochastic hybrid systems with application to communication networks,”Nonlinear Analysis: Theory, Methods & Applications, vol. 62, no. 8, pp. 1353–1383, 2005

    work page 2005

  14. [14]

    Spectral uncertainty propagation for generalized stochastic hybrid systems with applications to a bouncing ball,

    W. Wang and T. Lee, “Spectral uncertainty propagation for generalized stochastic hybrid systems with applications to a bouncing ball,” in Proceedings of the American Control Conference, July 2020, pp. 1803– 1808

    work page 2020

  15. [15]

    Uncertainty propagation for general stochastic hybrid systems on compact Lie groups,

    ——, “Uncertainty propagation for general stochastic hybrid systems on compact Lie groups,”SIAM Journal on Applied Dynamical Systems, vol. 21, no. 3, pp. 2215–2240, June 2022

    work page 2022

  16. [16]

    Uncertainty propagation of stochastic hybrid systems: a case study for types of jump,

    T. C, W. Clark, and T. Lee, “Uncertainty propagation of stochastic hybrid systems: a case study for types of jump,” inProceedings of the IFAC Symposium on Nonlinear Control Systems, July 2025

    work page 2025

  17. [17]

    Stochastic models for chemically reacting systems using polynomial stochastic hybrid systems,

    J. P. Hespanha and A. Singh, “Stochastic models for chemically reacting systems using polynomial stochastic hybrid systems,”International Journal of Robust and Nonlinear Control: IFAC-Affiliated Journal, vol. 15, no. 15, pp. 669–689, 2005

    work page 2005

  18. [18]

    Rubust approximation of the stochastic Koopman operator,

    M. Wanner and I. Mezic, “Rubust approximation of the stochastic Koopman operator,”SIAM Journal on Applied Dynamical Systems, vol. 21, no. 3, pp. 1930–1951, 2022

    work page 1930

  19. [19]

    N. A. Weiss, P. T. Holmes, and M. Hardy,A course in probability. Pearson Addison Wesley Boston, MA, USA:, 2006

    work page 2006

  20. [20]

    Geometric interpretation of Brownian motion on Riemannian manifolds,

    T. Lee and G. Chirikjian, “Geometric interpretation of Brownian motion on Riemannian manifolds,”arXiv, 2025

    work page 2025

  21. [21]

    Invariant forms in hybrid and impact systems and a taming of Zeno,

    W. Clark and A. Bloch, “Invariant forms in hybrid and impact systems and a taming of Zeno,”Archive for Rational Mechanics and Analysis, vol. 247, no. 2, p. 13, 2023

    work page 2023

  22. [22]

    Øksendal,Stochastic differential equations

    B. Øksendal,Stochastic differential equations. Springer, 2003

    work page 2003

  23. [23]

    Efficient implementation of weighted eno schemes,

    G.-S. Jiang and C.-W. Shu, “Efficient implementation of weighted eno schemes,”Journal of computational physics, vol. 126, no. 1, pp. 202– 228, 1996

    work page 1996

  24. [24]

    Weighted essentially non-oscillatory schemes,

    X.-D. Liu, S. Osher, and T. Chan, “Weighted essentially non-oscillatory schemes,”Journal of computational physics, vol. 115, no. 1, pp. 200– 212, 1994

    work page 1994

  25. [25]

    C. T. Kelley,Solving nonlinear equations with Newton’s method. SIAM, 2003. APPENDIX A. Divergence Theorem The divergence of a vector field can be obtained by the Lie derivative, i.e.,div µ(X)µ=L X µ. Also, from Cartan’s formula, we haveL X µ=di X µ+i X dµ, wheredµ= 0as it is a top form. By combining these,div µ(X)µ=di X µ. According to Stoke’s theorem, we...

    work page 2003