---
title: Tempered Stable Subordinator Overview
url: https://www.emergentmind.com/topics/tempered-stable-subordinator
type: topic
---

# Tempered Stable Subordinator Overview

A tempered stable subordinator is an increasing, pure-jump Lévy process that generalizes the classical stable subordinator by introducing exponential tempering of the Lévy measure. It retains the heavy-tailed, infinite-activity structure for small jumps but induces exponentially damped tails for large jumps, thereby ensuring the existence of all moments. This property makes tempered stable subordinators fundamental in advanced stochastic modeling in areas such as anomalous diffusion, finance, insurance, queueing theory, and applied probability.

## 1. Definition, Lévy–Khintchine Representation, and Laplace Exponent

Let $S(t)$, $t\ge0$, denote a tempered stable subordinator with stability index $\alpha\in(0,1)$ and tempering parameter $\lambda>0$. $S$ is a nondecreasing, pure-jump Lévy process whose Lévy measure is absolutely continuous and given by
\[
\nu(dx) = C x^{-1-\alpha} e^{-\lambda x}\, dx, \qquad x>0,
\]
with scale $C>0$ (often $C = \alpha/\Gamma(1-\alpha)$). The tails of the Lévy measure,
$\nu(x)\sim x^{-1-\alpha}$ for small $x$ and $\nu(x)\sim e^{-\lambda x}$ for large $x$, reflect the mixture of power-law and exponential decay [2303.07060, 1110.2868, 1912.12190, 1905.00192, 2004.02267].

The Laplace exponent is 
\[
\phi(s) = (s+\lambda)^\alpha - \lambda^\alpha,
\]
so that
\[
\mathbb{E}[e^{-s S(t)}] = \exp\left(- t \phi(s) \right).
\]
This Laplace exponent is a complete Bernstein function, ensuring infinite divisibility and strict monotonicity [2405.09966, 1203.0892]. In the $\lambda\to 0$ limit, one recovers the classical $\alpha$-stable subordinator. For $\alpha\to1$, $\phi(s)\to s,$ so $S_{1,\lambda}(t)=t$ is deterministic.

## 2. Distributional Properties and Series/Integral Representations

The one-dimensional marginal density $f_{\alpha,\lambda}(x,t)$ may be expressed as an exponential tilt of the stable density $f_\alpha(x,t)$:
\[
f_{\alpha,\lambda}(x,t) = \exp(-\lambda x + \lambda^\alpha t) f_\alpha(x,t),
\]
where $f_\alpha(x,t)$ is not available in closed form but can be represented by Fox $H$-functions or convergent series [1905.00192, 2004.02267, 1110.2868]. For the classical tempered stable subordinator:
\[
\int_0^\infty e^{-s x} f_{\alpha,\lambda}(x,t)\,dx = \exp(-t[(s+\lambda)^\alpha-\lambda^\alpha]).
\]
Laplace inversion or expansions using generalized Mittag–Leffler functions are common for analytic and numerical work.

### Asymptotic Behaviors

- **Small $x$:** $f_{\alpha,\lambda}(x,t)\sim x^{-1-\alpha}$, matching the behavior of the underlying stable process.
- **Large $x$:** $f_{\alpha,\lambda}(x,t)\sim e^{-\lambda x} x^{-1-\alpha}$, yielding exponentially damped tails [2004.02267, 1912.12190].
- **Moments:** All moments exist. 
  - Mean: $\mathbb{E}[S(t)] = \alpha\, t\, \lambda^{\alpha-1}$.
  - Variance: $\mathrm{Var}[S(t)] = \alpha(1-\alpha)t\,\lambda^{\alpha-2}$.

This exponential tempering is critical in distinguishing tempered from pure stable subordinators, whose moments may diverge.

## 3. Governing Equations and Fractional Dynamics

The PDF $f_{\alpha,\lambda}(x,t)$ solves a shifted, fractional Kolmogorov equation:
\[
\partial_t p(x,t) = -[(\lambda - \partial_x)^\alpha - \lambda^\alpha]\, p(x,t), \qquad p(x,0) = \delta(x),\; p(0,t) = 0,
\]
where $(\lambda - \partial_x)^\alpha$ denotes a shifted Riemann–Liouville fractional derivative. This equation encodes nonlocal, history-dependent dynamics, and interpolation between pure fractional and classical regimes [1905.00192, 2405.09966, 1912.12190].

For the inverse subordinator $E_{\alpha,\lambda}(t) = \inf\{u\geq 0: S(u) > t\}$, the density $h_{\alpha,\lambda}(u, t)$ satisfies a tempered time-fractional evolution equation. The Caputo–temperately tempered derivative:
\[
{}^{C}D^{\alpha,\lambda}_{t}u(t) = \frac{1}{\Gamma(1-\alpha)} e^{-\lambda t} \int_0^t (t-s)^{-\alpha} e^{\lambda s} u'(s)\, ds.
\]
Then
\[
{}^{C}D^{\alpha,\lambda}_{t} h_{\alpha,\lambda}(u,t) = -\partial_{u} h_{\alpha,\lambda}(u,t), \qquad h_{\alpha,\lambda}(u,0) = \delta(u).
\]
This operator interpolates between fractional and classical derivatives depending on the tempering parameter. 

## 4. Simulation and Numerical Methods

The absence of closed-form densities necessitates specialized simulation methods:

- **Poisson Mixture Scheme**: Approximates the subordinator law using scaled Poisson mixtures, with explicit error bounds in $L^p$ (Kolmogorov, Wasserstein metrics) [2409.09909].
- **Series Expansion Sampling**: Expansions in Fox $H$-function or Mittag–Leffler series allow for efficient approximation of both the subordinator and its inverse [1110.2868, 2405.09966].
- **Exact First-Passage Simulation**: For barrier-crossing functionals, a combination of Chambers–Mallows–Stuck sampling, Esscher tilting, and accept/reject (Esscher-rejection) yields fast exact samples for the first-passage time, undershoot, and overshoot, with explicit complexity bounds [2303.11964].
- **Finite Difference and Euler Schemes**: For path-level simulation and numerical evaluation of functionals, discretized variants of the integral representations and weak approximation of the Lévy measure are effective [2009.02583, 2409.09909].

These methods leverage the explicit Lévy structure and the complete monotonicity of the Laplace exponent, allowing both accuracy guarantees and stable asymptotic error rates.

## 5. Applications and Extensions

Tempered stable subordinators underpin a range of stochastic modeling frameworks:

- **Anomalous Diffusion**: Subordinate Brownian motion $Y(t) = B(S(t))$ exhibits subdiffusive-to-normal diffusive crossover, with mean squared displacement scaling as $t^\alpha$ at short times (subdiffusion) and $t$ at long times (normal) [1110.2868, 2004.02267].
- **Finance**: The Normal-Tempered-Stable (NTS) model and its Additive (ATS) extension use $S(t)$ as the stochastic clock for Brownian motion or more general semimartingale-based returns, enabling realistic modeling of volatility clustering and implied volatility surfaces with consistent moment structure [1909.07139, 2303.07060, 2009.02583].
- **Queueing and Insurance**: The ability to control tail behavior and finite activity of large jumps makes tempered stable subordinators useful as input processes for arrival models and claims modeling.
- **Structural Degradation**: The running average of a tempered stable subordinator defines a new infinite-activity process with closed-form cumulant and density representations, suitable for modeling cumulative degradation under random environments [2009.02583].

Extension to **mixtures of tempered stable subordinators** further generalizes the tail behavior. For a mixture of $n$ subordinators with different $(\alpha_i, \lambda_i)$ parameters, the Laplace exponent becomes a convex combination: $\sum_{i=1}^{n} c_i [(s + \lambda_i)^{\alpha_i} - \lambda_i^{\alpha_i}]$ [1905.00192].

## 6. Statistical Inference and Parameter Estimation

Since the density is not explicit, parametric estimation for tempered stable subordinators relies on:

- **Maximum Likelihood Estimation (MLE)**: Relies on series representation or fast Fourier inversion for evaluating the likelihood.
- **Generalized Method of Moments (GMM)**: Uses matching of theoretical cumulants, which are analytically available via:
  \[
  \kappa_m = \Gamma(m-\alpha)\, \frac{\delta}{\lambda^{m-\alpha}}, \quad m=1,2,\dots
  \]
- **Empirical Characteristic Function Estimation**: Fits observed characteristic functions to the model form.
- **Simulation-Based Methods**: For instance, calibration of NTS models for volatility surface fitting as in equity and energy markets [2303.07060].

MLEs achieve asymptotic normality and consistency under standard conditions; GMM provides computational advantages when adequate moment conditions and weighting are met. All moments are finite for $\lambda > 0$. For small samples, boundary solutions for $\alpha$ can result, highlighting the need for careful regularization and diagnostic checking.

## 7. Inverse Tempered Stable Subordinator and Related Processes

The inverse subordinator $E_{\alpha,\lambda}(t)$ is essential for modeling random time changes in non-Markovian models:

- **Laplace Transform and Series**: The joint density $h_{\alpha,\lambda}(u,t)$ admits series and contour integral representations involving the generalized Mittag–Leffler function [1410.1632, 2405.09966].
- **Asymptotic Regimes**: For $t \rightarrow 0$, $E[E_{\alpha,\lambda}(t)] \sim t^\alpha/\Gamma(1+\alpha)$, while for $t\to\infty$ the mean grows linearly, reflecting the exponentially tempered scaling of the parent process.
- **Governing PDEs**: Governing equations for the inverse density emerge as time-fractional evolution (Caputo–tempered) equations and as space-fractional PDEs in particular cases.

The inverse process governs "subordination" in fractional Fokker–Planck and "tempered" Kolmogorov equations, as well as in time-changed Hawkes processes and branched continuous-time random walks [2405.09966, 1410.1632].

---

Collectively, the tempered stable subordinator and its inverse form the probabilistic and analytic backbone of a wide range of "semi-heavy-tailed" models, retaining the rich structure of stable processes with exponentially suppressed extremes, robust path properties, and analytic/numerical tractability across applied probability and mathematical finance [2303.07060, 2405.09966, 1110.2868, 1912.12190, 2409.09909, 2303.11964, 2009.02583, 1905.00192, 1203.0892, 1410.1632, 2004.02267, 1909.07139].

Source: https://www.emergentmind.com/topics/tempered-stable-subordinator