Table of Contents
Fetching ...

Kernel Methods for Stochastic Dynamical Systems with Application to Koopman Eigenfunctions: Feynman-Kac Representations and RKHS Approximation

Boumediene Hamzi, Houman Owhadi, Umesh Vaidya

Abstract

We extend the unified kernel framework for transport equations and Koopman eigenfunctions, developed in previous work by the authors for deterministic systems, to stochastic differential equations (SDEs). In the deterministic setting, three analytically grounded constructions-Lions-type variational principles, Green's function convolution, and resolvent operators along characteristic flows--were shown to yield identical reproducing kernels. For stochastic systems, the Koopman generator includes a second-order diffusion term, transforming the first-order hyperbolic transport equation into a second-order elliptic-parabolic PDE. This fundamental change necessitates replacing the method of characteristics with probabilistic representations based on the Feynman--Kac formula. Our main contributions include: (i) extension of all three kernel constructions to stochastic systems via Feynman--Kac path-integral representations; (ii) proof of kernel equivalence under uniform ellipticity assumptions; (iii) a collocation-based computational framework incorporating second-order differential operators; (iv) error bounds separating RKHS approximation error from Monte Carlo sampling error; (v) analysis of how diffusion affects numerical conditioning; and (vi) connections to generator EDMD, diffusion maps, and kernel analog forecasting. Numerical experiments on Ornstein--Uhlenbeck processes, nonlinear SDEs with varying diffusion strength, and multi-dimensional systems validate the theoretical developments and demonstrate that moderate diffusion can improve numerical stability through elliptic regularization.

Kernel Methods for Stochastic Dynamical Systems with Application to Koopman Eigenfunctions: Feynman-Kac Representations and RKHS Approximation

Abstract

We extend the unified kernel framework for transport equations and Koopman eigenfunctions, developed in previous work by the authors for deterministic systems, to stochastic differential equations (SDEs). In the deterministic setting, three analytically grounded constructions-Lions-type variational principles, Green's function convolution, and resolvent operators along characteristic flows--were shown to yield identical reproducing kernels. For stochastic systems, the Koopman generator includes a second-order diffusion term, transforming the first-order hyperbolic transport equation into a second-order elliptic-parabolic PDE. This fundamental change necessitates replacing the method of characteristics with probabilistic representations based on the Feynman--Kac formula. Our main contributions include: (i) extension of all three kernel constructions to stochastic systems via Feynman--Kac path-integral representations; (ii) proof of kernel equivalence under uniform ellipticity assumptions; (iii) a collocation-based computational framework incorporating second-order differential operators; (iv) error bounds separating RKHS approximation error from Monte Carlo sampling error; (v) analysis of how diffusion affects numerical conditioning; and (vi) connections to generator EDMD, diffusion maps, and kernel analog forecasting. Numerical experiments on Ornstein--Uhlenbeck processes, nonlinear SDEs with varying diffusion strength, and multi-dimensional systems validate the theoretical developments and demonstrate that moderate diffusion can improve numerical stability through elliptic regularization.
Paper Structure (66 sections, 21 theorems, 166 equations, 1 figure, 2 algorithms)

This paper contains 66 sections, 21 theorems, 166 equations, 1 figure, 2 algorithms.

Key Result

Theorem 2.4

If $w$ is a left eigenvector of $A$ with eigenvalue $\lambda$, then $h(x)$ satisfies

Figures (1)

  • Figure 1: Numerical experiments for SDE Koopman eigenfunctions. (a) OU eigenfunction: exact recovery of $\phi(x)=x$. (b) OU semigroup property verification. (c) Quadratic nonlinearity eigenfunctions for different noise levels. (d) Semigroup verification for quadratic case ($\sigma=0.3$). (e) 2D linear eigenfunction contours with zero level set. (f) 2D linear semigroup verification. All semigroup errors are below 10%.

Theorems & Definitions (67)

  • Theorem 2.4: PDE for Nonlinear Correction
  • proof
  • Remark 2.5: Deterministic Limit
  • Definition 3.1: Admissible Kernels for Stochastic Problems
  • Example 3.2: Admissible Kernels
  • Definition 3.3: Variational Kernel
  • Theorem 3.4: Feynman--Kac Formula
  • proof
  • Corollary 3.5: Feynman--Kac Representation of $h$
  • proof
  • ...and 57 more