CONTENTS

Lévy Processes

Motivation: why this matters in quant finance

A Brownian motion model is powerful because it has independent and stationary increments. Those two assumptions make returns over non-overlapping periods easy to combine, and they are the reason characteristic functions, generators, and option-pricing equations stay tractable.

Continuity is a separate assumption. Lawler begins the jump-process chapter by asking what happens if we keep independent and stationary increments but stop requiring continuous paths. The resulting class is the Lévy processes. Brownian motion is still inside the class, but it is no longer the whole story.

For quantitative finance, this is the clean conceptual step from Black-Scholes toward jump models. A Lévy process is not "a fancy jump distribution"; it is the process-level structure that lets Brownian motion, Poisson processes, and compound Poisson shocks be handled with one set of tools.

The informal idea

A Lévy process has increments that are:

  1. Independent: what happens after time ss is independent of the past once you look at the increment.
  2. Stationary: the distribution of an increment depends on its length, not on calendar time.

Brownian motion adds continuous paths. Lévy processes do not require that. If continuity is dropped, the process may jump, while retaining the same clean time-additive structure that makes Brownian motion mathematically useful.

Lawler's first examples make this contrast concrete:

  • Brownian motion is the continuous-path Lévy process.
  • The Poisson process is the basic integer-valued jump Lévy process.
  • A compound Poisson process allows random jump sizes.
  • Hitting-time and Cauchy-process examples show that not every Lévy process looks Gaussian or Poisson.

Formal definitions

A stochastic process XtX_t is a Lévy process if it has independent and stationary increments: for s,t>0s,t>0,
Xs+tXsX_{s+t}-X_s

is independent of {Xr:rs}\{X_r:r\le s\} and has the same distribution as XtX0X_t-X_0.

Lawler also emphasises the path convention needed for a precise continuous-time object. The process is taken to be right-continuous with left limits, often called cadlag:
Xt=Xt+:=limqtXq,Xt:=limqtXq.X_t=X_{t+}:=\lim_{q\downarrow t}X_q, \qquad X_{t-}:=\lim_{q\uparrow t}X_q.

At a jump time, XtXtX_t\ne X_{t-} and the jump size is

ΔXt=XtXt.\Delta X_t=X_t-X_{t-}.

Infinitely divisible increments

If XtX_t is a Lévy process, then X1X_1 can be decomposed into nn i.i.d. increments:

X1=j=1n(Xj/nX(j1)/n).X_1=\sum_{j=1}^n \left(X_{j/n}-X_{(j-1)/n}\right).
So X1X_1 has an infinitely divisible distribution. Conversely, Lawler notes that every infinitely divisible distribution can arise as the time-one law of a Lévy process.

Two core examples are:

  • If XN(m,σ2)X\sim N(m,\sigma^2), then XX is the sum of nn independent N(m/n,σ2/n)N(m/n,\sigma^2/n) variables.
  • If XX is Poisson with mean λ\lambda, then XX is the sum of nn independent Poisson variables with mean λ/n\lambda/n.

This explains why normal and Poisson laws are the first two pillars of the Lévy family.

Characteristic exponents

Lawler uses characteristic exponents to express the additivity of independent increments. For a random variable XX,

E[eisX]=eΨX(s),ΨX(0)=0.\mathbb{E}[e^{isX}]=e^{\Psi_X(s)}, \qquad \Psi_X(0)=0.

If XX and YY are independent, then

ΨX+Y(s)=ΨX(s)+ΨY(s).\Psi_{X+Y}(s)=\Psi_X(s)+\Psi_Y(s).

For a Lévy process starting at the origin, stationary independent increments imply

ΨXt(s)=tΨX1(s).\Psi_{X_t}(s)=t\Psi_{X_1}(s).

This is the main tractability property. The distribution over time tt is encoded by multiplying the one-period exponent by tt.

Poisson and compound Poisson examples

For a Poisson random variable with mean λ\lambda,

Ψ(s)=λ(eis1).\Psi(s)=\lambda(e^{is}-1).

For a compound Poisson process with jump-size distribution μ#\mu^\# and rate λ\lambda, define the finite measure

μ=λμ#.\mu=\lambda\mu^\#.

Then the characteristic exponent is

Ψ(s)=R(eisy1)dμ(y).\Psi(s)=\int_{\mathbb{R}}(e^{isy}-1)\,d\mu(y).

This formula is the grounded version of the finance slogan "jumps add a tail term". The jump-size distribution enters through the integral, while the jump frequency is the total mass of μ\mu.

Generators

A Lévy process is Markov, so it has a generator

Lf(x)=limt0E[f(Xt)X0=x]f(x)t.Lf(x)=\lim_{t\downarrow 0}\frac{\mathbb{E}[f(X_t)\mid X_0=x]-f(x)}{t}.

For Brownian motion with drift mm and variance parameter σ2\sigma^2,

Lf(x)=mf(x)+12σ2f(x).Lf(x)=mf'(x)+\frac{1}{2}\sigma^2f''(x).

For a Poisson process with rate λ\lambda,

Lf(x)=λ[f(x+1)f(x)].Lf(x)=\lambda[f(x+1)-f(x)].

For a compound Poisson process with Lévy measure μ\mu,

Lf(x)=R[f(x+y)f(x)]dμ(y).Lf(x)=\int_{\mathbb{R}}[f(x+y)-f(x)]\,d\mu(y).

If two independent Lévy processes are added, their generators add. For a Brownian motion with drift plus an independent unit-jump Poisson process,

Lf(x)=mf(x)+12σ2f(x)+λ[f(x+1)f(x)].Lf(x)=mf'(x)+\frac{1}{2}\sigma^2f''(x)+\lambda[f(x+1)-f(x)].

Worked examples

Example 1: Brownian plus Poisson exponent

Let

Xt=mt+σWt+Nt,X_t=mt+\sigma W_t+N_t,

where NtN_t is independent Poisson with rate λ\lambda. The characteristic exponent of X1X_1 is

Ψ(s)=ims12σ2s2+λ(eis1).\Psi(s)=ims-\frac{1}{2}\sigma^2s^2+\lambda(e^{is}-1).

The first two terms come from the Gaussian component; the last term comes from unit jumps.

Example 2: compound Poisson mean and martingale compensation

Let Jt=i=1NtYiJ_t=\sum_{i=1}^{N_t}Y_i with finite first moment and measure μ=λμ#\mu=\lambda\mu^\#. Then

mJ=Rydμ(y)=λE[Y1],m_J=\int_{\mathbb{R}}y\,d\mu(y)=\lambda\mathbb{E}[Y_1],

and

Mt=JttmJM_t=J_t-tm_J

is a martingale. This compensation is central when jump models are used for P&L: it separates predictable jump drift from unexpected jump variation.

Example 3: jump quadratic variation

For a compound Poisson process,

Jt=st(ΔJs)2.\langle J\rangle_t=\sum_{s\le t}(\Delta J_s)^2.

So a path with two jumps of sizes 0.03-0.03 and 0.050.05 has jump quadratic variation

0.032+0.052=0.0034.0.03^2+0.05^2=0.0034.

That is a realised path quantity, not a deterministic clock.

Common confusions and pitfalls

"Lévy process means jump process." No. Brownian motion is a Lévy process. Jumps enter when the continuous-path assumption is dropped.
"Stationary increments mean the process itself is stationary." No. The increment law depends only on the time length. The level XtX_t usually changes distribution as tt grows.
"Every Lévy process has continuous paths except at rare jumps." Not from the definition alone. Lawler introduces cadlag paths so discontinuities can be handled, and later theory distinguishes finite-activity and more general jump behaviour.
"The Lévy measure in a compound Poisson process is a probability measure." In Lawler's compound Poisson setup, μ=λμ#\mu=\lambda\mu^\# is a finite measure whose total mass is the jump rate. The probability distribution is μ#\mu^\#.
"Characteristic exponents are decorative notation." They carry the main additivity. Independent components add exponents, and time scales them by tt.

Where this goes next

References

  • Lawler, G. F. (2023). Stochastic Calculus: An Introduction with Applications. Ch. 6 §6.1 (Lévy processes), §6.2 (Poisson process), §6.3 (Compound Poisson process).

title: Lévy Processes order: 7 verified: false latest_update_date: 26/04/2026 progress: completed difficulty: advanced prerequisites:

  • learn/math/probability/stochastic-processes/brownian-motion
  • learn/math/probability/stochastic-processes/jump-diffusion-processes
  • learn/math/probability/foundations/random-variables

Lévy Processes

Motivation: why this matters in quant finance

Lévy processes provide the standard framework for continuous-time models with independent and stationary increments. They include Brownian Motion, Jump-Diffusion Processes, Poisson processes, compound Poisson processes, and infinite-activity jump models under one language.

In quantitative finance, the payoff is practical: Lévy models can keep the tractable independent-increment structure of Brownian motion while allowing discontinuities, skew, and heavier tails. They form the mathematical backbone of exponential Lévy option models and clarify exactly what is added when Black-Scholes is extended with jumps.

Lawler's construction frames the contrast sharply: Brownian motion is the continuous-path member of the Lévy family, while Poisson and compound Poisson processes are the basic discontinuous members. The lesson's job is to understand that common structure before moving to generators, characteristic exponents, or jump-pricing equations.

Definition and Basic Properties

Formal Definition

A stochastic process {Xt}t0\{X_t\}_{t \geq 0} is a Lévy process if:
  1. Initial condition: X0=0X_0 = 0 almost surely
  2. Independent increments: For 0t1<t2<<tn0 \leq t_1 < t_2 < \cdots < t_n, the random variables Xt1,Xt2Xt1,,XtnXtn1X_{t_1}, X_{t_2} - X_{t_1}, \ldots, X_{t_n} - X_{t_{n-1}} are mutually independent
  3. Stationary increments: Xt+hXtX_{t+h} - X_t has the same distribution as XhX_h for all t,h>0t, h > 0
  4. Stochastic continuity: limh0P(Xt+hXt>ε)=0\lim_{h \to 0} \mathbb{P}(|X_{t+h} - X_t| > \varepsilon) = 0 for all ε>0\varepsilon > 0 and t0t \geq 0
  5. Càdlàg paths: Sample paths are right-continuous with left limits

Examples

  • Brownian Motion: Xt=σWt+μtX_t = \sigma W_t + \mu t
  • Compound Poisson Process: Xt=i=1NtYiX_t = \sum_{i=1}^{N_t} Y_i
  • Gamma Process: Increasing Lévy process with Gamma-distributed increments
  • Inverse Gaussian Process: Subordinator used in time-changed models
  • Variance Gamma Process: Time-changed Brownian motion

The Lévy-Khintchine Formula

Characteristic Function

The fundamental theorem characterizing Lévy processes states that the characteristic function of XtX_t has the form:

E[eiuXt]=etψ(u)\mathbb{E}[e^{iu X_t}] = e^{t\psi(u)}
where the Lévy exponent ψ(u)\psi(u) is given by:
ψ(u)=iγu12σ2u2+R{0}(eiuz1iuz1z<1)ν(dz)\psi(u) = i\gamma u - \frac{1}{2}\sigma^2 u^2 + \int_{\mathbb{R} \setminus \{0\}} \left(e^{iuz} - 1 - iuz\mathbf{1}_{|z|<1}\right) \nu(dz)

Lévy Triplet

Every Lévy process is uniquely determined by its Lévy triplet (γ,σ2,ν)(\gamma, \sigma^2, \nu):
  • γR\gamma \in \mathbb{R}: Drift coefficient
  • σ20\sigma^2 \geq 0: Gaussian (diffusion) coefficient
  • ν\nu: Lévy measure on R{0}\mathbb{R} \setminus \{0\} satisfying R(1z2)ν(dz)<\int_{\mathbb{R}} (1 \wedge |z|^2) \nu(dz) < \infty

Lévy Measure

The Lévy measure ν\nu encodes the jump structure:
  • ν({0})=0\nu(\{0\}) = 0 (no jumps of size zero)
  • ν(A)\nu(A) represents the expected number of jumps in set AA per unit time
  • z>1ν(dz)<\int_{|z|>1} \nu(dz) < \infty ensures finite jump activity for large jumps
  • z1z2ν(dz)<\int_{|z| \leq 1} |z|^2 \nu(dz) < \infty ensures finite quadratic variation from small jumps

Lévy-Itô Decomposition

Decomposition Formula

Any Lévy process can be written as:

Xt=γt+σWt+limε0[0<st,ΔXs>εΔXs+0<st,ΔXsε(ΔXsE[ΔXs])]X_t = \gamma t + \sigma W_t + \lim_{\varepsilon \to 0} \left[\sum_{0 < s \leq t, |\Delta X_s| > \varepsilon} \Delta X_s + \sum_{0 < s \leq t, |\Delta X_s| \leq \varepsilon} (\Delta X_s - \mathbb{E}[\Delta X_s])\right]

More precisely:

Xt=γt+σWt+Xt(1)+Xt(2)X_t = \gamma t + \sigma W_t + X_t^{(1)} + X_t^{(2)}

where:

  • γt\gamma t: Linear drift
  • σWt\sigma W_t: Brownian motion component
  • Xt(1)X_t^{(1)}: Compensated sum of small jumps
  • Xt(2)X_t^{(2)}: Sum of large jumps

Small and Large Jumps

  • Small jumps (z1|z| \leq 1): Infinite activity, compensated to have zero mean
  • Large jumps (z>1|z| > 1): Finite activity, uncompensated

This decomposition shows that Lévy processes are the most general processes with independent increments.

Special Cases and Examples

1. Subordinators

A Lévy process with non-decreasing paths. Examples:

  • Gamma process: ΓtGamma(αt,β)\Gamma_t \sim \text{Gamma}(\alpha t, \beta)
  • Inverse Gaussian process: Used in variance gamma models
  • Stable subordinator: Heavy-tailed increasing process

2. Symmetric α-Stable Processes

Characterized by stability under addition and heavy tails:

ψ(u)=cuα,0<α2\psi(u) = -c|u|^\alpha, \quad 0 < \alpha \leq 2

When α=2\alpha = 2, this reduces to Brownian motion.

3. Variance Gamma Process

Xt=θGt+σWGtX_t = \theta G_t + \sigma W_{G_t}

where GtG_t is a Gamma subordinator. This creates infinite activity pure-jump processes popular in option pricing.

4. Normal Inverse Gaussian (NIG)

Has Lévy density:

ν(x)=αδπxK1(αx)eβx\nu(x) = \frac{\alpha\delta}{\pi|x|} K_1(\alpha|x|) e^{\beta x}

where K1K_1 is the modified Bessel function of the first kind.

Time Changes and Subordination

Definition

A subordinator is an increasing Lévy process GtG_t with G0=0G_0 = 0. Time-changing a process YtY_t by a subordinator creates:
Xt=YGtX_t = Y_{G_t}

Examples

  1. Variance Gamma: Xt=WGtX_t = W_{G_t} where GtG_t is Gamma
  2. Normal Inverse Gaussian: Uses inverse Gaussian subordinator
  3. CGMY: Uses a different subordination scheme

Financial Interpretation

Time changes model varying market activity:

  • Business time: Markets are more active during certain periods
  • Transaction time: Time measured in number of trades
  • Information time: Time measured by information arrival

Infinitely Divisible Distributions

Definition

A random variable XX is infinitely divisible if for every n1n \geq 1, there exist i.i.d. random variables X1(n),,Xn(n)X_1^{(n)}, \ldots, X_n^{(n)} such that:
X=dX1(n)++Xn(n)X \stackrel{d}{=} X_1^{(n)} + \cdots + X_n^{(n)}

Connection to Lévy processes

  • The distribution of XtX_t for any Lévy process is infinitely divisible.
  • Every infinitely divisible distribution corresponds to a Lévy process.
  • Normal and Poisson laws are the canonical examples: Brownian motion decomposes normal increments, while Poisson processes decompose event counts.

This is the structural reason Lévy processes are not an arbitrary model class. They are exactly the continuous-time processes whose fixed-horizon increments can be split consistently into any number of smaller independent increments.

Option Pricing with Lévy Processes

Exponential Lévy Models

Asset prices follow:

St=S0eXtS_t = S_0 e^{X_t}

where XtX_t is a Lévy process with appropriate drift adjustment.

Risk-Neutral Pricing

Under the risk-neutral measure, the discounted price must be a martingale:
ertSt=ertS0eXte^{-rt}S_t = e^{-rt}S_0 e^{X_t}

This requires the Lévy exponent to satisfy:

ψ(i)=r\psi(-i) = r

Fourier Methods

Option prices can be computed using characteristic functions:

C=erT0(STK)+fST(s)dsC = e^{-rT} \int_0^\infty (S_T - K)^+ f_{S_T}(s) ds

where fSTf_{S_T} is recovered via inverse Fourier transform.

  1. Black-Scholes: Geometric Brownian motion
  2. Merton Jump-Diffusion: Brownian motion + compound Poisson
  3. Variance Gamma: Pure jump model
  4. CGMY: Infinite activity model
  5. NIG: Solvable in terms of special functions

Simulation Methods

Direct Simulation

  1. Simulate large jumps: Poisson process for ν({z>ε})\nu(\{|z| > \varepsilon\})
  2. Simulate small jumps: Use series representation or approximation
  3. Add Gaussian component: Standard Brownian motion
  4. Add drift: Deterministic term

Series Representation

For infinite activity processes:

Xt=σWt+i=1zi1UitX_t = \sigma W_t + \sum_{i=1}^{\infty} z_i \mathbf{1}_{U_i \leq t}

where UiU_i are arrival times and ziz_i are jump sizes.

Acceptance-Rejection

For processes with known Lévy densities, use acceptance-rejection to simulate jump sizes.

Estimation from Data

Maximum Likelihood

For discretely observed Lévy processes, maximize:

L(θ)=i=1nfΔ(XtiXti1;θ)L(\theta) = \prod_{i=1}^n f_{\Delta}(X_{t_i} - X_{t_{i-1}}; \theta)

where fΔf_\Delta is the transition density over interval Δ\Delta.

Method of Moments

Match sample moments to theoretical moments derived from the characteristic function.

Generalized Method of Moments (GMM)

Use multiple moment conditions based on the Lévy-Khintchine formula.

Applications in Risk Management

Value at Risk (VaR)

Lévy models better capture:

  • Heavy tails: More accurate extreme quantile estimation
  • Asymmetry: Different behavior for gains vs. losses
  • Jump risk: Sudden large moves

Expected Shortfall (ES)

ESα=E[XX>VaRα]\text{ES}_\alpha = \mathbb{E}[X | X > \text{VaR}_\alpha]

Lévy models provide more realistic tail risk measures.

Operational Risk

Compound Poisson processes model operational loss frequency and severity.

Advanced Topics

Fluctuation Theory

Studies the behavior of suprema, infima, and first passage times of Lévy processes.

Lévy Copulas

Extend the concept to multivariate settings with dependent marginal Lévy processes.

Path-Dependent Options

Barrier options, lookback options, and Asian options under Lévy dynamics.

Credit Risk Models

Jump-to-default models use Lévy processes to model credit events.

Common confusions and pitfalls

"Lévy process means jump process." No. Brownian motion is a Lévy process. Jumps enter only when the continuous-path assumption is dropped.
"Stationary increments mean the process is stationary." No. The increments have time-homogeneous laws; the level XtX_t usually changes distribution with tt.
"Infinite activity means infinitely many large jumps." No. Infinite-activity models can have infinitely many small jumps on finite intervals while still having finite large-jump activity.
"A Lévy model is automatically risk-neutral." No. It is a probabilistic model class. Risk-neutrality requires a measure under which discounted prices are martingales.
"The Lévy measure is an ordinary probability distribution." Not generally. It measures jump intensity by size, can have infinite total mass near zero, and is constrained by integrability conditions rather than total probability one.

Where this goes next

Lévy processes unify many concepts in stochastic finance:

References

  • Lawler, G. F. (2023). Stochastic Calculus: An Introduction with Applications. Ch. 6 §6.1 (Lévy processes), §6.2 (Poisson process), §6.3 (Compound Poisson process).