Table of Contents
Fetching ...

Simulation of infinite-dimensional diffusion bridges

Thorben Pieper-Sethmacher, Frank van der Meulen, Aad van der Vaart

TL;DR

This work develops a principled method to sample infinite-dimensional diffusion bridges for semilinear SPDEs conditioned on a finite-dimensional observation $L X_T = y$. By leveraging an exponential change of measure with a Doob $h$-transform, the authors construct a tractable guided process $X^{\circ}$ whose law is absolutely continuous w.r.t. the true bridge $X^{\star}$, enabling efficient Metropolis-Hastings and related inference techniques. They prove limit-behavior and absolute-continuity results, provide explicit constructions using the OU process, and demonstrate the approach on stochastic reaction-diffusion equations, including diagonalisable and neural-field models, with numerical illustrations and practical MH algorithms. The framework advances sampling and smoothing for SPDE bridges and has potential applications in partially observed SPDE parameter estimation and state estimation. Overall, the paper extends finite-dimensional diffusion-bridge methods to the SPDE setting and offers a versatile toolkit for inference in infinite-dimensional stochastic systems.

Abstract

Let X be the mild solution to a semilinear stochastic partial differential equation. In this article, we develop methodology to sample from the infinite-dimensional diffusion bridge that arises from conditioning X on a linear transformation LXT of the final state XT at some time T > 0. This solves a problem that has so far not been attended to in the literature. Our main contribution is the derivation of a path measure that is absolutely continuous with respect to the path measure of the infinite-dimensional diffusion bridge. This lifts previously known results for stochastic ordinary differential equations to the setting of infinite-dimensional diffusions and stochastic partial differential equations. We demonstrate our findings through numerical experiments on stochastic reaction-diffusion equations.

Simulation of infinite-dimensional diffusion bridges

TL;DR

This work develops a principled method to sample infinite-dimensional diffusion bridges for semilinear SPDEs conditioned on a finite-dimensional observation . By leveraging an exponential change of measure with a Doob -transform, the authors construct a tractable guided process whose law is absolutely continuous w.r.t. the true bridge , enabling efficient Metropolis-Hastings and related inference techniques. They prove limit-behavior and absolute-continuity results, provide explicit constructions using the OU process, and demonstrate the approach on stochastic reaction-diffusion equations, including diagonalisable and neural-field models, with numerical illustrations and practical MH algorithms. The framework advances sampling and smoothing for SPDE bridges and has potential applications in partially observed SPDE parameter estimation and state estimation. Overall, the paper extends finite-dimensional diffusion-bridge methods to the SPDE setting and offers a versatile toolkit for inference in infinite-dimensional stochastic systems.

Abstract

Let X be the mild solution to a semilinear stochastic partial differential equation. In this article, we develop methodology to sample from the infinite-dimensional diffusion bridge that arises from conditioning X on a linear transformation LXT of the final state XT at some time T > 0. This solves a problem that has so far not been attended to in the literature. Our main contribution is the derivation of a path measure that is absolutely continuous with respect to the path measure of the infinite-dimensional diffusion bridge. This lifts previously known results for stochastic ordinary differential equations to the setting of infinite-dimensional diffusions and stochastic partial differential equations. We demonstrate our findings through numerical experiments on stochastic reaction-diffusion equations.

Paper Structure

This paper contains 25 sections, 14 theorems, 119 equations, 11 figures, 2 algorithms.

Key Result

Theorem 2.2

Let $h: [0,T) \times H \longrightarrow \mathbb{R}_{>0}$ satisfy the following for any $S < T$: Then $h$ defines a unique change of measure $\mathbb{P}^h$ on $(\Omega, \mathcal{F}_T)$ via Additionally, let $h$ satisfy the following: Then $X$ under $\mathbb{P}^h$ is a mild solution to the SPDE where $W^h$ is a $\mathbb{P}^h$-cylindrical Wiener process.

Figures (11)

  • Figure 1: Forward simulation of $X$. Left: Heatmap of a sample path $X(t,\xi)$ of Equation \ref{['eq: MichaelisMenten']} with parameters as specified in \ref{['eq: MM_params']}. Right: The states $X_0$ and $X_T$ and the observation $y = P_k X_T$, $k=10$.
  • Figure 2: Heatmap of a sample path $X^{\circ}(t,\xi)$ of the guided process corresponding to \ref{['eq: MichaelisMenten']} with conditioning state $y = P_k X_T, k = 10$.
  • Figure 3: Heatmap of the mean sample path of $100$ samples of the estimated diffusion bridge $X^{\star}(t,\xi)$ of Equation \ref{['eq: MichaelisMenten']} returned by Algorithm \ref{['alg: MH_sampler1']}.
  • Figure 4: Paths of two spectral modes of the diffusion bridge samples of Equation \ref{['eq: MichaelisMenten']} returned by Algorithm \ref{['alg: MH_sampler1']}. Every $500$-th sample is shown. Green paths represent 'earlier' samples in the Markov chain, whereas blue paths represent 'later' samples. The orange path is the spectral mode of the data generating process $X$.
  • Figure 5: Sample paths in the spectral modes $1$ and $2$ obtained by forward sampling the SPDE \ref{['eq: MichaelisMenten']} with parameters as specified in \ref{['eq: MM_params']}. Only paths that satisfy \ref{['eq: MMexperimentsconditioning']} are kept.
  • ...and 6 more figures

Theorems & Definitions (34)

  • Theorem 2.2: Theorem 3.5, Piepersethmacher24Classexponentialchangesmeasure
  • Remark 2.3
  • Proposition 2.4: Existence of the diffusion bridge
  • Remark 2.5
  • Proposition 2.6: Existence of the guided process
  • Remark 2.7
  • Theorem 3.2
  • Theorem 3.4
  • Remark 3.5
  • Lemma 3.6
  • ...and 24 more