--- title: Time-Changed Compound Poisson-Gamma Subordinators url: https://www.emergentmind.com/topics/time-changed-compound-poisson-gamma-subordinators type: topic --- # Time-Changed Compound Poisson-Gamma Subordinators A time-changed compound Poisson–Gamma subordinator is a specific class of stochastic process constructed by employing the Gamma subordinator with Poissonian randomization, and then utilizing this process to randomize the clock (i.e., the time argument) of another process, most often a Lévy or Poisson process. This framework unifies finite-activity jump processes with subordinated dynamics and is essential in modeling bursty or clustered phenomena in a range of probabilistic systems, with deep connections to anomalous diffusion, aggregated stochastic modeling, and non-local governing equations [1707.00523], [1806.03833], [2106.12201]. ## 1. Definition and Construction Let \( N_\lambda(t) \) be a Poisson process with rate \( \lambda > 0 \). For a Gamma subordinator \( G_{\alpha,\beta}(s) \) — with shape parameter \( \alpha > 0 \) and rate \( \beta > 0 \) — the compound Poisson–Gamma subordinator is defined as \[ S(t) = \sum_{n=1}^{N_\lambda(t)} G_n \equiv G(N_\lambda(t)), \quad S(0) = 0, \] where \( \{G_n\} \) are i.i.d. Gamma(\(\alpha,\beta\)) random variables [1707.00523, Sec. 3]. The resulting process \( S(t) \) is a pure-jump subordinator with finite activity. When this subordinator is used to direct (“time-change”) another independent stochastic process — often a Poisson process \( N_{\lambda_1} \) — one obtains processes such as \( X(t) = N_{\lambda_1}(S(t)) \). More generally, one can replace the outer process with any Lévy process [1707.00523], [1806.03833]. ## 2. Analytical Properties ### Laplace Exponent and Lévy Measure The Lévy–Khintchine exponent of the compound Poisson–Gamma subordinator is \[ \psi(\theta) = -\log \mathbb{E} [e^{-\theta S(1)}] = \lambda \beta^\alpha\left[\beta^{-\alpha} - (\beta+\theta)^{-\alpha}\right]. \] Its Lévy measure is \[ \nu(du) = \lambda \frac{\beta^\alpha}{\Gamma(\alpha)} u^{\alpha-1} e^{-\beta u}\, du. \] For \(\alpha=1\), \( S(t) \) becomes a compound Poisson process with exponential jumps [1707.00523, (3.1)–(3.2)]. ### Marginal Distribution For \( t > 0 \), \[ P\{S(t) = 0\} = e^{-\lambda t}, \quad P\{S(t) \in ds\} = e^{-\lambda t} \delta_0(ds) + e^{-\lambda t - \beta s} s^{-1} \Phi(\alpha, 0; \lambda t (\beta s)^\alpha)\, ds \] where \( \Phi \) is the Wright function [1707.00523, (3.2)]. Alternatively, \[ f_{S(t)}(s) = e^{-\lambda t - \beta s} \sum_{n=1}^\infty \frac{(\lambda t \beta^\alpha)^n}{n! \Gamma(\alpha n)} s^{\alpha n - 1}. \] ## 3. Time-Changed Poisson Processes ### Distributional Law Letting \( X(t) = N_{\lambda_1}(S(t)) \), the marginal law is given by \[ p_k(t) = P\{X(t) = k\} = e^{-\lambda t} \frac{k! \lambda_1^k}{(\lambda_1+\beta)^k} \sum_{n=1}^\infty \frac{(\lambda t \beta^\alpha)^n\,\Gamma(\alpha n + k)}{n! \Gamma(\alpha n) (\lambda_1+\beta)^{\alpha n}}, \quad k \ge 1, \] with \[ p_0(t) = \exp\left\{ -\lambda t \left[1 - \left(\frac{\beta}{\lambda_1+\beta}\right)^{\!\alpha}\,\right]\right\} \] [1707.00523, (4.1)–(4.2)], [1806.03833]. ### Moment Generating Function The mgf of \( X(t) \) is given by \[ M(u; t) = \mathbb{E}[e^{u X(t)}] = \exp\left\{-\lambda t \left[1 - \beta^\alpha (\beta + \lambda_1(1-e^u))^{-\alpha}\right]\right\}. \] This expression leverages the subordinator's explicit Laplace exponent and Bochner subordination [1707.00523, (4.12)]. ### Hitting Times and Ruin Probabilities Hitting times, or first passage times to state \( m \), have the distribution: \[ P\{T_m<\infty\} = \frac{\lambda_1^m}{m!(\lambda_1+\beta)^m \left[1-\left(\frac{\beta}{\lambda_1+\beta}\right)^\alpha\right]} \sum_{n=1}^\infty \left(\frac{\beta}{\lambda_1+\beta}\right)^{\alpha n} \frac{\Gamma(\alpha n + m)}{\Gamma(\alpha n)} \] indicating that the process may never reach state \( m \) with strictly positive probability [1806.03833]. ## 4. Generalizations and Iterations ### Extensions A linear deterministic drift can be incorporated: \[ Y(t) = G_{\alpha,\beta}(N_\lambda(t) + a t), \] with Laplace exponent \[ \phi_Y(u) = \alpha a \ln\!\left(1+\frac{u}{\beta}\right) + \lambda\left[1 - (1+u/\beta)^{-\alpha}\right]. \] Time-changing a Poisson process by \( Y(t) \) yields distributions analogous to the basic time-changed model, with parameters modified by the drift \( a \) [1806.03833]. ### Iterated Compositions One may construct \[ G_1(N(t)), \quad G_2(N(G_1(N(t)))), \ldots \] For \( \alpha = 1 \), a “semigroup” property holds: iterated compound Poisson–exponential subordinators again yield a process within the same class, with effective parameters composed algebraically. For general \( \alpha \), Laplace exponents add hierarchically, leading to stable or quasi-stable limiting behavior in some regimes [1806.03833]. ## 5. Inverse Processes and Skellam-Type Constructions ### Inverse Subordinators For \( \alpha = 1 \), the inverse of a compound Poisson–exponential process, \( Y(u) = \inf\{t: S(t)>u\} \), has explicit density and Laplace transform given by Bessel and exponential terms. This structure produces time-changed Poisson processes \( N_1(Y(t)) \) whose distributions are described by generalized Mittag-Leffler functions [1707.00523, Section 5]. ### Skellam-Type Differences Two variants arise: - **Type I (Common Clock):** \( S_I(t) = N_1(S(t)) - N_2(S(t)) \) - **Type II (Independent Clocks):** \( S_{II}(t) = N_1(X_1(t)) - N_2(X_2(t)) \), with \( X_i(t) \) iid subordinators Explicit pmfs and mgfs are obtained using integration against the compound Poisson–Gamma laws, with forms involving Wright and Bessel functions [1707.00523, eqs. (4.20), (4.33)]. ## 6. Connections to Anomalous Diffusion and Fractional Calculus ### Lower-Incomplete Gamma Subordinator A related process is defined via the lower-incomplete gamma function: for \( a \in (0,1] \), \[ \phi(\eta) = \alpha \gamma(\alpha; \eta), \quad \gamma(\alpha;\eta) = \int_0^\eta w^{\alpha-1}e^{-w} dw \] [2106.12201, (3.1)]. As \( \alpha \to 1 \), this recovers the standard Poisson process; for \( \alpha < 1 \), it serves as a finite-activity approximation to the stable subordinator. ### Governing Equations Time-changing Markov semigroups by the compound Poisson–Gamma or incomplete gamma subordinator yields non-local integro-differential equations, serving as a discrete or mesoscopic analog to fractional-derivative evolution equations \[ \frac{\partial}{\partial t}g(x,t) = \frac{1}{\Gamma(1-\alpha)} \int_\varepsilon^\infty (T_s g(\cdot, t) - g(\cdot, t))(s-\varepsilon)^{-\alpha-1} ds \] [2106.12201, §6]. ### Sub-Diffusive Scaling For a Brownian motion time-changed by the lower-incomplete gamma subordinator, only fractional moments exist; e.g., for \( 0