Technopedia

The technologist’s encyclopedia

Lévy processes

A Lévy process is a stochastic process with stationary, independent increments starting at zero: the law of Xt−XsX_t - X_s depends only on t−st - s, and increments over disjoint intervals are independent. It is the canonical model for random motion that combines smooth drift, diffusion, and sudden jumps — the single framework that unifies Brownian motion, Poisson processes, and everything in between.

When it is relevant. Whenever a quantity accumulates by many small, independent contributions over time — a stock price with both continuous noise and discrete crashes, the total claims in an insurance portfolio, the net displacement of a particle in a turbulent medium — a Lévy process is the natural model. Its power is that one object, the Lévy–Khintchine triplet, encodes the entire law.

A Lévy process is a continuous-time Markov process, and its compensated version is a martingale. It builds on Random variables and expectation and Limit theorems.

Authoritative sources

Definition (cheatsheet)

A common source of confusion. The drift parameter bb is not an absolute number: it depends on the truncation function hh you choose. Changing hh (for instance to h(x)=x 1∣x∣≤1h(x) = x\,\mathbf{1}_{|x|\le 1} or h(x)=x/(1+x2)h(x) = x/(1+x^2)) shifts bb by a compensating constant so that the triplet always describes the same process. Only the combination of bb and ν\nu is invariant, not bb alone.

The Lévy–Khintchine triplet

Definition. A Lévy process is completely determined by its triplet (b,σ2,ν)(b, \sigma^2, \nu): a real drift bb, a diffusion coefficient σ2≥0\sigma^2 \ge 0, and a Lévy measure ν\nu (a measure on R∖{0}\mathbb{R}\setminus\{0\} with ∫(1∧x2) ν(dx)<∞\int (1 \wedge x^2)\,\nu(dx) < \infty). The characteristic exponent is

ψ(θ)=ibθ−12σ2θ2+∫R∖{0}(eiθx−1−iθx h(x)) ν(dx)\psi(\theta) = i b \theta - \tfrac{1}{2}\sigma^2\theta^2 + \int_{\mathbb{R}\setminus\{0\}} (e^{i\theta x} - 1 - i\theta x\, h(x))\,\nu(dx).

Details

The whole theory hangs on one fact: stationary independent increments.

  1. One increment is the whole process. Stationary independent increments mean you only ever need to know the law of one increment, say X1X_1; the rest follows by convolution. The characteristic function turns convolution into multiplication, so E[eiθXt]=etψ(θ)\mathbb{E}[e^{i\theta X_t}] = e^{t\psi(\theta)}: the whole process is one function ψ\psi.
  2. The Lévy–Khintchine theorem says every such ψ\psi has exactly the triplet form above: a linear drift term bb, a Gaussian term σ2\sigma^2, and an integral against the Lévy measure ν\nu for jumps. Conversely, any choice of (b,σ2,ν)(b, \sigma^2, \nu) with ∫(1∧x2) ν(dx)<∞\int (1 \wedge x^2)\,\nu(dx) < \infty produces a valid Lévy process.
  3. Finite vs infinite activity. If ν(R)<∞\nu(\mathbb{R}) < \infty, there are finitely many jumps per unit time (a compound Poisson part). If ν(R)=∞\nu(\mathbb{R}) = \infty but ∫∣x∣<1∣x∣ ν(dx)<∞\int_{|x|<1} |x|\,\nu(dx) < \infty, jumps occur infinitely often but remain summable — the process still has finite variation. Dropping the last condition gives infinite-variation processes, where the “jumps” are so numerous and small that they blur into a continuous-looking drift.
  4. Examples. The Poisson process has ν=δ1\nu = \delta_1 (unit jumps at rate 1); Brownian motion has ν=0\nu = 0; the Cauchy process has a heavy-tailed ν\nu with infinite activity and infinite variation.

The truncation function h(x)=x 1∣x∣<1h(x) = x\,\mathbf{1}_{|x|<1} is what makes the integral converge. For large xx the term eiθx−1e^{i\theta x} - 1 is harmless, but near x=0x = 0 the expansion eiθx−1≈iθxe^{i\theta x} - 1 \approx i\theta x would make ∫iθx ν(dx)\int i\theta x\,\nu(dx) diverge when there are infinitely many small jumps. Subtracting iθx h(x)i\theta x\, h(x) cancels that divergence, and the leftover is absorbed into the drift bb — which is precisely why bb depends on hh.

Example problem

A compound Poisson process has jump rate λ=3\lambda = 3 per unit time and jump sizes YY uniformly distributed on {1,2}\{1, 2\} (so P(Y=1)=P(Y=2)=1/2P(Y=1) = P(Y=2) = 1/2). Write down its Lévy–Khintchine triplet and compute E[Xt]E[X_t] and Var(Xt)\mathrm{Var}(X_t).

Solution

The Lévy measure is ν=λ (12δ1+12δ2)=32δ1+32δ2\nu = \lambda\,(\tfrac{1}{2}\delta_1 + \tfrac{1}{2}\delta_2) = \tfrac{3}{2}\delta_1 + \tfrac{3}{2}\delta_2, there is no diffusion (σ2=0\sigma^2 = 0), and with the standard truncation h(x)=x 1∣x∣<1h(x) = x\,\mathbf{1}_{|x|<1} (which vanishes on {1,2}\{1,2\}) the drift is b=0b = 0. So the triplet is (0,0,32δ1+32δ2)(0, 0, \tfrac{3}{2}\delta_1 + \tfrac{3}{2}\delta_2).

For a compound Poisson process, E[Xt]=λt E[Y]E[X_t] = \lambda t\, E[Y] and Var(Xt)=λt E[Y2]\mathrm{Var}(X_t) = \lambda t\, E[Y^2]. Here E[Y]=32E[Y] = \tfrac{3}{2} and E[Y2]=12(1)+12(4)=52E[Y^2] = \tfrac{1}{2}(1) + \tfrac{1}{2}(4) = \tfrac{5}{2}, so

E[Xt]=3t⋅32=92 t,Var(Xt)=3t⋅52=152 t.E[X_t] = 3t \cdot \tfrac{3}{2} = \tfrac{9}{2}\,t, \qquad \mathrm{Var}(X_t) = 3t \cdot \tfrac{5}{2} = \tfrac{15}{2}\,t.

Sample paths

Definition. The Lévy–Itô decomposition writes every Lévy process as Xt=b t+σWt+JtX_t = b\,t + \sigma W_t + J_t: a deterministic drift, a Brownian part σWt\sigma W_t, and a pure-jump part JtJ_t built from a Poisson random measure. The three terms are independent, and the path you see is their sum.

Details

Each ingredient has a very different look, and the decomposition is exactly what the four sample paths below show:

The four panels below are one realization of each extreme of the family: pure diffusion, a unit-jump Poisson process, a compound Poisson process with random jump sizes, and pure drift.

Brownian motion (σ > 0, no jumps) Poisson process (ν = δ₁) Compound Poisson (finite activity, random sizes) Deterministic drift (b t) time t

The interactive simulation below lets you turn each ingredient on and off and watch the resulting path.

Example problem

A process has triplet (b,σ2,ν)=(1,4,2 δ0.5)(b, \sigma^2, \nu) = (1, 4, 2\,\delta_{0.5}). Describe its path and compute E[Xt]E[X_t] and Var(Xt)\mathrm{Var}(X_t).

Solution

The path is a drift of slope 11, plus a Brownian motion with variance parameter σ2=4\sigma^2 = 4 (so σ=2\sigma = 2), plus a compound Poisson process of rate 22 with all jumps equal to 0.50.5.

For the drift, E=tE = t and Var=0\mathrm{Var} = 0. For the Brownian part, E=0E = 0 and Var=σ2t=4t\mathrm{Var} = \sigma^2 t = 4t. For the compound Poisson part with jump size c=0.5c = 0.5 and rate λ=2\lambda = 2: E=λt c=2t⋅0.5=tE = \lambda t\, c = 2t \cdot 0.5 = t, and Var=λt c2=2t⋅0.25=0.5 t\mathrm{Var} = \lambda t\, c^2 = 2t \cdot 0.25 = 0.5\,t. The three parts are independent, so expectations add and variances add:

E[Xt]=t+0+t=2t,Var(Xt)=0+4t+0.5 t=4.5 t.E[X_t] = t + 0 + t = 2t, \qquad \mathrm{Var}(X_t) = 0 + 4t + 0.5\,t = 4.5\,t.

Members of the family

PageDriftDiffusionJumps
Deterministic driftyesnono
Brownian motionoptionalyesno
Poisson processoptionalnounit jumps
Compound Poisson processoptionalnofinite activity, arbitrary sizes
Variance-gamma processoptionaloptionalinfinite activity
CGMY processoptionalnoinfinite activity, no diffusion
Stable processoptionaloptionalinfinite activity, self-similar

The first three rows are the building blocks; the last three are infinite-activity processes whose small jumps are too numerous to count individually. The stable processes are the self-similar members of the family, which ties into the Subordinators section below.

Subordinators

Definition. A subordinator is a non-decreasing Lévy process — a Lévy process with only positive jumps and no Brownian part, so its triplet is (b,0,ν)(b, 0, \nu) with ν\nu supported on (0,∞)(0, \infty). It is a “random clock”: a non-negative process (Tt)(T_t) that can be used to reparameterize (subordinate) another process.

Details

A subordinator is what you get when you strip a Lévy process down to its jump part and forbid negative jumps. Because it only moves upward, it can be read as a time change: given any process (Ys)(Y_s), the composition YTtY_{T_t} is called the subordination of YY by TT.

The simplest subordinator is a Poisson process (unit positive jumps). A stable subordinator of index α∈(0,1)\alpha \in (0,1) has Lévy measure ν(dx)=c x−1−α dx\nu(dx) = c\, x^{-1-\alpha}\,dx on (0,∞)(0,\infty) and is self-similar: Tct=dc1/αTtT_{ct} \stackrel{d}{=} c^{1/\alpha} T_t. Subordinators are the bridge between the Lévy family and the self-similarity theme: a stable process is a Lévy process that is self-similar, and a stable subordinator is the time-change that produces it.

Example. Let (Nt)(N_t) be a Poisson process of rate λ\lambda and (Ws)(W_s) a standard Brownian motion, independent. The subordinated process WNtW_{N_t} is a Lévy process: it is a compound Poisson process whose jumps are the Brownian increments between successive Poisson arrivals.

Example problem

Is Tt=2t+NtT_t = 2t + N_t (where (Nt)(N_t) is a rate-1 Poisson process) a subordinator? If so, write its triplet.

Solution

Yes. It is non-decreasing (the drift 2t2t and the counting process NtN_t both only increase), it has stationary independent increments, and it starts at 00. Its triplet is (b,σ2,ν)=(2,0,δ1)(b, \sigma^2, \nu) = (2, 0, \delta_1): drift b=2b = 2, no diffusion, and unit jumps at rate 11.

Applications

Quizz

Try it yourself

A Lévy process has characteristic exponent

ψ(θ)=i θ−12 θ2+∫R∖{0}(eiθx−1−iθx 1∣x∣<1) 12δ3(dx)\psi(\theta) = i\,\theta - \tfrac{1}{2}\,\theta^2 + \int_{\mathbb{R}\setminus\{0\}} (e^{i\theta x} - 1 - i\theta x\,\mathbf{1}_{|x|<1})\,\tfrac{1}{2}\delta_{3}(dx).

  1. Read off the triplet (b,σ2,ν)(b, \sigma^2, \nu).
  2. Is the jump part finite- or infinite-activity?
  3. Compute E[X1]E[X_1] and Var(X1)\mathrm{Var}(X_1).
Solution

1. Comparing with ψ(θ)=ibθ−12σ2θ2+∫(eiθx−1−iθx h(x)) ν(dx)\psi(\theta) = i b \theta - \tfrac{1}{2}\sigma^2\theta^2 + \int (e^{i\theta x} - 1 - i\theta x\, h(x))\,\nu(dx) with h(x)=x 1∣x∣<1h(x) = x\,\mathbf{1}_{|x|<1}: the drift is b=1b = 1, the diffusion is σ2=1\sigma^2 = 1, and the Lévy measure is ν=12δ3\nu = \tfrac{1}{2}\delta_3.

2. ν(R)=12<∞\nu(\mathbb{R}) = \tfrac{1}{2} < \infty, so the jump part is finite-activity (a compound Poisson process of rate 12\tfrac{1}{2} with jumps of size 33).

3. The three parts are independent. Drift: E=1E = 1 and Var=0\mathrm{Var} = 0. Brownian (σ2=1\sigma^2 = 1): E=0E = 0 and Var=1\mathrm{Var} = 1. Compound Poisson (rate λ=12\lambda = \tfrac{1}{2}, jump c=3c = 3): E=λc=12⋅3=32E = \lambda c = \tfrac{1}{2}\cdot 3 = \tfrac{3}{2} and Var=λc2=12⋅9=92\mathrm{Var} = \lambda c^2 = \tfrac{1}{2}\cdot 9 = \tfrac{9}{2}. Adding:

E[X1]=1+0+32=52,Var(X1)=0+1+92=112.E[X_1] = 1 + 0 + \tfrac{3}{2} = \tfrac{5}{2}, \qquad \mathrm{Var}(X_1) = 0 + 1 + \tfrac{9}{2} = \tfrac{11}{2}.