---
title: Skellam Point Process Overview
url: https://www.emergentmind.com/topics/skellam-point-process
type: topic
---

# Skellam Point Process Overview

Searching arXiv for recent and foundational papers on Skellam point processes and related variants.
A Skellam point process is a signed counting process built from differences of Poisson counting mechanisms. In its classical form, it is defined by
\[
S(t)=N_1(t)-N_2(t),
\]
where \(N_1\) and \(N_2\) are independent homogeneous Poisson processes with intensities \(\lambda_1,\lambda_2>0\). At fixed time \(t\), \(S(t)\) has the Skellam distribution, while as a process it is a pure-jump Lévy process with upward and downward jumps corresponding to the two Poisson drivers. Subsequent work has extended this construction to order-\(k\) and arbitrary-jump models, non-homogeneous and multivariate settings, inverse-subordinator and stable-subordinator time changes, spatial and multiparameter random fields, and biologically motivated reset mechanisms for neural spike trains [1806.00277], [2504.07672], [2206.09903].

## 1. Classical construction and basic probabilistic structure

The classical Skellam process is the difference of two independent Poisson processes,
\[
S(t)=N_1(t)-N_2(t), \qquad t\ge 0.
\]
Its one-dimensional law is
\[
s_k(t)=\Pr\{S(t)=k\} = e^{-t(\lambda_1+\lambda_2)} \left(\frac{\lambda_1}{\lambda_2}\right)^{k/2} I_k\!\left(2t\sqrt{\lambda_1\lambda_2}\right), \qquad k\in\mathbb Z,
\]
where \(I_k\) is the modified Bessel function of the first kind. The process has mean and variance determined by the two Poisson rates, and its infinitesimal behavior is that a jump of \(+1\) occurs at rate \(\lambda_1\) and a jump of \(-1\) occurs at rate \(\lambda_2\) [1806.00277].

This process is also a compound Poisson process. One representation is
\[
S_{\underline{\lambda}}(t)\stackrel{d}{=}\sum_{k=1}^{N_{\lambda_1+\lambda_2}(t)}X_k,
\]
where \(X_k\in\{1,-1\}\) with
\[
P(X_k=1)=\frac{\lambda_1}{\lambda_1+\lambda_2}, \qquad P(X_k=-1)=\frac{\lambda_2}{\lambda_1+\lambda_2}.
\]
That representation makes explicit that the Skellam process is a signed counting process rather than merely a static distribution on \(\mathbb Z\) [2208.06376].

Its state probabilities satisfy the forward system
\[
\frac{d}{dt}s_k(t) = \lambda_1\big(s_{k-1}(t)-s_k(t)\big)-\lambda_2\big(s_k(t)-s_{k+1}(t)\big), \qquad k\in\mathbb{Z},
\]
with \(s_0(0)=1\) and \(s_k(0)=0\) for \(k\neq 0\). The corresponding Lévy measure is concentrated at \(\pm1\), so the classical model is the two-sided unit-jump member of a much wider family of signed jump processes [2003.09471].

## 2. Time changes, fractional dynamics, and nonlocal evolution

A major line of development replaces calendar time by a random operational time. For an inverse subordinator
\[
Y_f(t)=\inf\{s\ge 0: H_f(s)>t\},
\]
the time-changed process
\[
Z(t)=S(Y_f(t))
\]
has marginal probabilities
\[
r_k^f(t)=\Pr\{Z(t)=k\} = \int_0^\infty s_k(u)\,l_f(t,u)\,du,
\]
where \(l_f(t,u)\) is the density of the inverse subordinator. These marginals satisfy the nonlocal system
\[
{}^{f}\!D_t r_k^f(t) = \lambda_1\big(r_{k-1}^f(t)-r_k^f(t)\big) - \lambda_2\big(r_k^f(t)-r_{k+1}^f(t)\big), \qquad k\in\mathbb Z,
\]
with the generalized Caputo-Djrbashian convolution derivative \( {}^{f}\!D_t\). When \(f(x)=x^\alpha\), this reduces to the known fractional Skellam equation [1806.00277].

The literature distinguishes two fractional Skellam constructions. In type 1, one takes the difference of two independent fractional Poisson processes,
\[
Y_{\underline{\nu},\underline{\lambda}}(t)=N_{\nu_1,\lambda_1}(t)-N_{\nu_2,\lambda_2}(t).
\]
In type 2, one time-changes the Skellam process itself,
\[
Z_{\nu,\underline{\lambda}}(t)=S_{\underline{\lambda}}(L_\nu(t)),
\]
with \(L_\nu\) the inverse stable subordinator. The two models have different moderate deviation behavior, and the paper on noncentral moderate deviations emphasizes that type 2 is “less random” in a sense because there is only one inverse-stable randomization instead of two [2208.06376].

Further generalization leads to the generalized fractional Skellam process
\[
S_\alpha(t)=S(Y_\alpha(t)),
\]
where \(S\) is a generalized Skellam process and \(Y_\alpha\) is an independent inverse stable subordinator. Its state probabilities satisfy a Caputo-fractional system, its probability generating function has Mittag-Leffler form, and the one-dimensional distributions are not infinitely divisible. The same paper derives long-range dependence for the generalized fractional Skellam process and short-range dependence for its increment process [2107.08307].

For order-\(k\) models, the fractional Skellam process of order \(k\),
\[
S_k^\alpha(t)=S_k(Y_\alpha(t)),
\]
inherits grouped jump sizes \(1,\dots,k\) while replacing the ordinary time derivative in the state equations by a Caputo derivative. Its correlations decay as \(t^{-\alpha}\), which establishes long-range dependence [2103.09187].

## 3. Generalized, order-\(k\), and non-homogeneous families

A broad generalization replaces the two jump sizes \(\pm1\) by an arbitrary fixed jump alphabet. One formulation is
\[
S(t)=\sum_{i\in I} i\,N_i(t), \qquad t\ge 0,
\]
where \(I\subset \mathbb R\setminus\{0\}\) and the \(N_i\) are independent non-homogeneous Poisson processes with cumulative rates
\[
\Lambda_i(t)=\int_0^t \lambda_i(s)\,ds.
\]
The probability generating function is
\[
G_S(t;u)=\mathbb{E}[u^{S(t)}] =\exp\!\left(-\sum_{i\in I}\Lambda_i(t)\,(1-u^i)\right).
\]
This family is closed under scaling and summation, admits law-of-large-numbers and central-limit-type limits, and in the homogeneous case has a compound Poisson representation with Lévy measure concentrated on the jump sizes \(I\) [2504.07672].

The order-\(k\) construction specializes this idea to symmetric positive and negative jumps up to size \(k\). The Skellam process of order \(k\) is
\[
S^k(t)=N_1^k(t)-N_2^k(t),
\]
where \(N_1^k\) and \(N_2^k\) are independent Poisson processes of order \(k\). Its marginal law retains the Skellam-Bessel form but with time scale multiplied by \(k\), and its Lévy measure is supported on \(\pm1,\dots,\pm k\). This model generalizes the classical Skellam process by allowing jump magnitudes beyond one [2003.09471].

A different generalization uses generalized counting processes with jump sizes \(1,\dots,k\) and possibly time-dependent rates. The non-homogeneous generalized Skellam

Source: https://www.emergentmind.com/topics/skellam-point-process