---
title: Probabilistic Tube Bounds in Stochastic Systems
url: https://www.emergentmind.com/topics/probabilistic-tube-bounds
type: topic
---

# Probabilistic Tube Bounds in Stochastic Systems

A probabilistic tube bound refers to a mathematically rigorous enclosure of all possible trajectories of a dynamical system, random field, or stochastic process, such that the trajectories remain within a "tube" or region around a reference path, set, or manifold with high probability. These constructions quantify uncertainty and robustness, yielding finite-sample or asymptotic probability guarantees for verification, control, learning, and geometric estimation in contexts ranging from stochastic differential equations to reachability analysis and prediction interval regression.

## 1. Core Mathematical Formulations

Probabilistic tube bounds come in several canonical forms, but always involve the high-probability containment of random trajectories within a time-/space-varying set, typically parameterized by a tube width $r(t)$ or set $\Omega_t$, and measured against a reference trajectory or manifold.

### General Template

For a dynamical system $X_t$ (with or without stochasticity), the aim is to certify
\[
\Pr\left\{ X_t \in \Omega_t \text{ for all } t\in[0,T] \right\} \geq 1-\delta.
\]
Key settings include:

- **Continuous-depth ODEs (Neural ODEs):**
  - For $x'(t) = f(x(t)), x(t_0) \in B_0 = B(x_0, r_0)$, the reachable set at time $t$ is $R(t; B_0) = \{\phi(t;x) : x\in B_0\}$.
  - A probabilistic tube is constructed so that $R(t_j; B_0) \subseteq B(\phi(t_j; x_0), \hat M_j)$ holds simultaneously for all $j$ with total probability at least $1-\delta$ [2107.08467].

- **Stochastic Differential Equations:**
  - For $dX_t = f(X_t) dt + g(X_t) dW_t$ and reference $x_t$ from $\dot x_t = f(x_t)$,
    \[
    \Pr\left( \sup_{t\in[0,T]} \| X_t - x_t \| \leq r \right) \geq 1-\delta
    \]
    where $r$ is linked to system noise, contraction, and time horizon via martingale concentration [2503.03328].

- **Target Tube Reachability:**
  - For stochastic/deterministic systems $x_{k+1} = f_k(x_k,u_k,w_k)$ and a "target tube" $\mathcal{T} = \{\mathcal{T}_k\}$,
    \[
    \Pr^\pi\left\{ x_k \in \mathcal{T}_k ~\forall k=0,\ldots,N \right\} \geq \alpha
    \]
    is maximized or underapproximated over initial states [1810.05217, 2405.02827, 2505.03482].

- **Gaussian/Random Fields:**
  - The tube method (Euler characteristic heuristic) estimates $P\left\{\sup_{u\in M} f(u) \geq c \right\}$ by integrating the probability of hitting threshold $c$ over tubes around $M$; the relative error decays exponentially in the subexponential tail case [2507.11154].

## 2. Probabilistic Tube Bounds in Reachability and Verification

### Sampling-Based and Optimization Approaches

- **GoTube Algorithm (ODEs and continuous-depth NNs):**
  - At each time $t_j$, draws $N$ i.i.d. initial conditions, simulates trajectories, and computes the maximal deviation $m_n$. A concentration-based padding $\varepsilon_j$ (from DKW inequality) is added:
    \[
    \hat{M}_j = m_n + \sqrt{(1/2N)\ln(1/\delta_j)}
    \]
    so that
    \[
    \Pr(M_j > \hat{M}_j) \leq \delta_j.
    \]
  - A union bound over time steps gives a full-horizon tube with overall confidence $1-\delta$ [2107.08467].

- **Convex Optimization (Discrete-Time):**
  - For control systems, a Bellman recursion or open-loop convex program yields polytopic underapproximations of the set of initial states from which staying in a target tube is achievable with level $\alpha$. The method leverages log-concavity for set interpolation and convexity of reachable regions [1810.05217].

- **Set-Erosion in Stochastic Safety:**
  - A "probabilistic tube" $r(t)$ quantifies deviation due to stochastic noise. Eroding the safety set by $r(t)$ reduces stochastic verification to a deterministic one, with the tube bound
    \[
    \Pr\left(\sup_{t\in[0,T]} \| X_t - x_t \| \leq r \right) \geq 1-\delta
    \]
    where typically $r = \sqrt{2\Lambda(T)\ln(1/\delta)}$ for suitable system-dependent $\Lambda(T)$ [2503.03328].

### Techniques for Multi-Agent and MPC Settings

- **Probabilistic Reachable Tubes for Multi-Agent Systems:**
  - The error components of each agent evolve as decoupled processes, for which confidence regions (ellipsoidal, via Chebyshev) are propagated through time. By recursively composing one-step probabilistic reachable sets, a high-probability tube for the entire network state is achieved, with an explicit lower bound on joint probability via Boole's inequality [2405.02827].

- **Learning-Based Tube MPC:**
  - Tube boundaries are learned via quantile regression or "Tube Loss" objectives to encapsulate trajectory uncertainty, yielding tubes that (with prescribed coverage) contain a fraction $t$ of all possible regression targets or system states [2412.06853, 2002.01587, 2505.03482].

## 3. Probabilistic Tube Bounds in Hypoelliptic and Geometric Settings

### Degenerate Diffusions and Non-Isotropic Tubes

- **Diffusions under Hörmander Condition:**
  - For SDEs where the diffusion matrix may be degenerate, but the Lie algebra generated by diffusion and bracket terms is full-rank, the tube is defined via a non-isotropic norm $\|y\|_{A_R(x)}$ reflecting different propagation speeds in direct and bracket directions. Two-sided exponential bounds for trajectory containment probabilities are established:
    \[
    \exp\left(-C\int_0^T H_t\left(\frac{1}{R_t} + |\phi_t|^2\right)dt\right)
    \leq
    \Pr\left(\sup_{0\leq t\leq T} \|X_t - x_t(\phi)\|_{A_{R_t}(t,x_t(\phi))} \leq 1\right)
    \]
    where $H_t$ captures local geometry and degeneracy [1607.04544, 1202.4771].

### Tube Volume and Random Geometry

- **Volume of Tubular Neighborhoods (Algebraic Geometry):**
  - For a manifold or variety $V$ and tube of radius $\varepsilon$, the volume (and thus probability that a random point is within $\varepsilon$ of $V$) can be bounded in terms of the variety's algebraic complexity (degree, codimension), via Weyl's tube formula, Crofton formula, and Bézout-type degree bounds. These estimates provide upper bounds for the probability of small-distance proximity in high-dimensional spaces [1210.3742].

- **Random Cantor Sets and Tube Intersections:**
  - In geometric measure theory, random fractal sets can satisfy uniform upper bounds on Hausdorff measure within tubes, i.e.,
    \[
    \sup_{T} \frac{\mathcal{H}^{d-1}(E \cap T)}{w(T)^t} < \infty
    \]
    for all $t < d-1$, directly quantifying how mass concentrates within arbitrarily thin tubes [1410.1183].

## 4. Practical Construction and Computation of Probabilistic Tubes

### Algorithmic Elements

- **Sampling, Padding, and Concentration:**
  - Sampling-based methods (e.g., GoTube) utilize random or quasi-random exploration of initial condition balls, updating tube radii with statistical concentration inequalities (DKW/Hoeffding) to maintain probabilistic validity [2107.08467].

- **Convex Optimization:**
  - Reach set underapproximations leverage epigraphical and line search convex programs, iteratively computing polytopic boundaries and interpolating coverage levels to obtain anytime probabilistic certificates [1810.05217].

- **Receding Horizon and Learning-Based Updates:**
  - In control, tube cross-sections are refined via recursive propagation or online data-driven approaches (linear programs or neural quantile regressors), providing adaptive, scenario-based guarantees on containment and recursive feasibility [2405.02827, 2505.03482, 2002.01587].

### Parameterization and Tightness

- **Control of Tightness and Level:**
  - Parameters such as the tightness factor $\mu$, number of samples $N$, or tube width regularization in loss objectives allow practitioners to trade off between statistical conservatism and tube width. Probability mass can be allocated across time or state dimensions using the union bound [2107.08467, 2412.06853].

- **Explicit Formulas:**
  - In certain linear/stochastic settings, closed-form expressions relate tube radius explicitly to system noise, contraction rates, time horizon, and target failure rates [2503.03328].

## 5. Theoretical Guarantees and Limitations

### Asymptotic and Finite-Sample Guarantees

- **Exponential Decay and Convergence:**
  - As the number of samples $N \to \infty$, statistical padding $\varepsilon_j$ vanishes, and the coverage probability approaches the target level for both GoTube and convex polytope reach set methods [2107.08467, 1810.05217].
  - For random fields and tube methods, the relative error $\Delta(c)$ of tube formula approximation decays exponentially for subexponential tails but not for regularly varying heavy-tailed marginals, where the method may remain biased [2507.11154].

- **Non-Isotropic, Degenerate Cases:**
  - In hypoelliptic diffusions, tube bounds account for different scaling in tangent and bracket directions; the result is strictly stronger than results obtained through isotropic (Euclidean) norm tubes and is locally equivalent to the intrinsic control (Carathéodory) metric [1607.04544, 1202.4771].

### Limitations

- **Sharpness and Failure Regimes:**
  - In geometric measure settings, bounds on the mass in tubes are sharp up to the threshold $t = d-1$; above this, tube-nullity intervenes and upper bounds cannot hold [1410.1183].
  - For random field tube methods, regularly varying tails invalidate tail-accuracy, requiring alternative Bonferroni-type bounds or correction terms [2507.11154].
  - Volume-of-tube bounds become ineffective for varieties with high complexity or in high dimension due to exponential constants, although they remain polynomial in degree for fixed $n$ [1210.3742].

## 6. Applications and Impact

- **Safety Verification & Model Checking:** Probabilistic tubes are used for formal safety verification in neural ODEs, nonlinear stochastic systems, and controller synthesis, allowing tractable reduction of chance constraints to deterministic tightening [2107.08467, 2503.03328, 2405.02827, 2505.03482].

- **Learning and Forecasting:** Tube Loss and quantile-based algorithms deliver direct construction of optimal prediction intervals and forecast tubes in neural networks, achieving correct empirical coverage and minimal prediction width efficiently [2412.06853, 2002.01587].

- **Geometric Probability and Statistics:** Probabilistic tube bounds control rare event probabilities, estimate small-volume neighborhoods in algebraic geometry, and enable precise calculation of tail exceedance probabilities for high-threshold random fields [1210.3742, 2507.11154].

- **Random Graphs and Percolation Theory:** Fluctuation lower bounds for first-passage percolation use tube-slicing and anti-concentration arguments to provide the first polynomial variance lower bounds in high-dimensional geometries [2204.06592].

- **Geometric Measure Theory:** The intersection behavior of random fractals with tubular neighborhoods answers open questions on mass concentration, applying probabilistic tube bounds in measure-theoretic and fractal contexts [1410.1183].

## 7. Summary Table: Main Approaches to Probabilistic Tube Bounds

| Context/Setting                      | Construction/Method               | Probabilistic Guarantee/Bound                                        |
|--------------------------------------|-----------------------------------|---------------------------------------------------------------------|
| Neural ODEs, Continuous Depth        | GoTube: Sampling + DKW Padding    | $\Pr(\text{all } t_j:~R(t_j;B_0)~\text{in tube}) \geq 1-\delta$     |
| Stochastic Control, Reachability     | Convex Polytope, Dynamic Prog.    | $\Pr(x_k \in \mathcal{T}_k~\forall k) \geq \alpha$                  |
| SDEs, Nonlinear Stochastic Systems   | Martingale/Concentration Analysis | $\Pr(\sup_{t} \|X_t - x_t\| \leq r ) \geq 1-\delta$                 |
| Learning-Based Forecasting           | Tube Loss/Quantile Regression     | Empirical coverage $\to t$, width minimized in a single objective   |
| Diffusions, Hypoelliptic Settings    | Non-Isotropic Norm Tubes          | Two-sided exponential probability bounds in bracket-adapted norm     |
| Algebraic Geometry, Random Fields    | Tube/Euler Characteristic Formula | Excursion probability, with relative error $\to 0$ (subexp. tails)  |

Probabilistic tube bounds thus provide a unifying principle for the explicit, computable, and probability-quantified enclosure of system trajectories, spanning domains from stochastic verification, robust control, learning-based forecasting, to geometric analysis of high-dimensional or fractal sets. The mathematical techniques combine sampling, convex optimization, concentration inequalities, martingale theory, geometric measure, and curvature integration, yielding robust guarantees in both classical and contemporary high-dimensional regimes.

Source: https://www.emergentmind.com/topics/probabilistic-tube-bounds