---
title: Poissonian Malliavin Calculus
url: https://www.emergentmind.com/topics/poissonian-malliavin-calculus
type: topic
---

# Poissonian Malliavin Calculus

Poissonian Malliavin calculus is the stochastic differential calculus of variations formulated for Poisson random measures, Lévy processes with jumps, and general pure-jump processes. Its central objects are discrete difference operators and their duals, which contrast fundamentally with the derivations and structure underlying classical Malliavin calculus on Wiener (Gaussian) space. The resulting calculus is fundamental for quantitative limit theorems, stochastic analysis of SDEs with jumps, sensitivity formulas (“Greeks”) in finance, stochastic partial differential equations driven by Lévy noise, absolute continuity results, and the Malliavin–Stein program for Gaussian approximations.

## 1. Construction of the Poisson Space and Chaos Decomposition

Given a probability space \((\Omega, \mathcal{F}, \mathbb{P})\) carrying a Poisson random measure \(N\) on \([0,\infty)\times \mathbb{R}\) with Lévy measure \(\nu\), the compensator is \(\tilde N(dt,dx) = N(dt,dx) - dt\,\nu(dx)\) [1712.04098]. Every \(F \in L^2(\Omega)\) admits a unique Wiener–Itô chaos expansion:
\[
F = \sum_{q=0}^\infty I_q(f_q),
\]
where \(I_q(f_q)\) are multiple compensated Poisson integrals with symmetric, square-integrable kernels \(f_q\) over \((\mathbb{R}_+\times \mathbb{R})^q\) and orthogonal decomposition [1712.04098, 1703.07259, 1605.07413]. Isometries and product/contracting formulas for these multiple integrals provide a Fock-space structure analogous to the Gaussian case [2110.04877].

## 2. Fundamental Operators: Difference, Divergence, and OU Generator

### Gradient (Difference/Malliavin) Operator \(D\)

For measurable \(F(\eta)\), the discrete gradient is
\[
D_z F := F(\eta + \delta_z) - F(\eta),
\]
acting as a closed, unbounded operator from \(L^2(\Omega)\) into \(L^2(\Omega \times [0, \infty) \times \mathbb{R})\), with domain those \(F\) whose chaos kernels \(f_q\) satisfy \(\sum_{q=1}^\infty q \, q! \|f_q\|^2 < \infty\). On the chaos expansion, \(D_z F = \sum_{q=1}^\infty q\, I_{q-1}(f_q(z, \cdot))\) [1712.04098, 1111.2140, 1703.07259].

### Difference Rules

Unlike the Gaussian case, the product and chain rules incorporate jump terms: for \(F, G \in \mathrm{Dom}\, D\), 
\[
D_z(FG) = G D_z F + F D_z G + x D_z F D_z G,
\]
and for a smooth \(\varphi\),
\[
D_z \varphi(F) = \varphi(F + D_z F) - \varphi(F) \quad \text{in the pure jump case} [2502.00793].
\]

### Divergence (Skorohod) Operator \(\delta\)

\(\delta\) is the formal adjoint of \(D\). For a process \(u \in L^2(\Omega \times [0,T] \times \mathbb{R})\), if \(F \in \mathbb{D}^{1,2}\),
\[
\E\left[F \delta(u)\right] = \E\left[\langle DF, u \rangle_{L^2(dt\,\nu(dx))}\right],
\]
and if \(u(z) = \sum_{q \ge 0} I_q(f_q(z,\cdot))\) with \(\sum (q+1)! \|f_q\|^2 < \infty\), then \(\delta(u) = \sum_{q \ge 0} I_{q+1}(f_q)\) [1712.04098, 1703.07259].

### Ornstein–Uhlenbeck Generator \(L\) and Pseudo-Inverse \(L^{-1}\)

\(L\) acts on chaos expansions via
\[
L F = -\sum_{q \ge 1} q I_q(f_q),
\]
with domain determined by \(\sum q^2 q! \|f_q\|^2 < \infty\). The fundamental commutation relation \(-\delta D = L\) holds on a suitable domain, and the pseudo-inverse is given by \(L^{-1} F = -\sum_{q \ge 1} (1/q) I_q(f_q)\), so \(L L^{-1}F = F - \E[F]\) [1712.04098, 1703.07259].

## 3. Advanced Calculus: Chain Rule, Energy Bracket, Gamma Calculus

The Poissonian chain rule is substantially more complex than in Gaussian analysis. For a multivariate smooth function \(\varphi(F_1,\ldots,F_d)\) and points \(z=(t,x)\),
\[
D_z(\varphi(F)) = \sum_{i=1}^d \partial_i \varphi(F) D_z F_i + \frac{1}{2} \sum_{i,j=1}^d \partial^2_{ij} \varphi(F) D_z F_i D_z F_j + \frac{x}{2} \sum_{i,j,k=1}^d \partial^3_{ijk} \varphi(F + \theta x D_z F) D_z F_i D_z F_j D_z F_k,
\]
for some \(\theta \in (0,1)\), with the last term reflecting jump compensation absent in Wiener analysis [1712.04098].

The carré du champ and energy bracket are extended via the Poisson version of the Gamma calculus. For \(F, G\), 
\[
\Gamma(F,G) = \frac{1}{2}(L(FG) - F LG - G LF)
\]
and
\[
[u,v]_\Gamma := \frac{1}{2} \int u(z) \otimes v(z) \nu(dz) + \frac{1}{2} \int (1-D^-_z)u(z)\otimes(1-D^-_z)v(z) \eta(dz),
\]
provide higher-order variance/covariance control in limiting theorems [1910.10181, 2110.04877].

## 4. Limit Theorems, Quantitative Approximation, and Stein–Malliavin Bounds

The “Nourdin–Peccati (NP) bound” for the normal approximation in Poisson space states that for centered, unit variance \(F\in\mathrm{Dom}\,D\), the Wasserstein distance satisfies
\[
d_W(F, N) \leq C \sqrt{\mathrm{Var}( \langle DF, -D L^{-1} F \rangle ) },
\]
with explicit expressions for the remainder involving higher-order derivatives and the jump terms [1712.04098]. A second-order Poincaré inequality further refines this to
\[
d_W(F, N) \leq C \left( \E[\|DF\|^4]^{1/4} \E[\|D^2F\|^4_{\text{op}}]^{1/4} + \E[\|D^2F\|^2_{\text{HS}}] \right),
\]
providing the explicit asymptotic rate of convergence, often \(O(n^{-1/2})\) or \(O(n^{-1})\) [1712.04098, 1111.2140, 1940.10181].

In the Hilbert space–valued setting, functional limit theorems and fourth-moment bounds admit clean formulations involving kernel contractions and Gamma operators, extending classical “fourth-moment phenomena” to the Poisson framework [2110.04877].

## 5. Applications: SDEs with Jumps, SPDEs, Sensitivity, and Hawkes Processes

### SDEs and SPDEs with Poisson Noise

The Poisson Malliavin calculus supports rigorous sensitivity analysis and error bounds for both finite-dimensional SDEs with jumps and infinite-dimensional SPDEs (particularly with α-stable drivers) [1703.07259, 2502.00793, 2510.05300]. The integration by parts and commutation relations 
\[
D_{t,x}\int_0^T\int_U \Phi(s,y) \tilde N(ds,dy) = \int_0^T\int_U D_{t,x}\Phi(s,y) \tilde N(ds,dy) + \Phi(t,x)
\]
enable sharp error estimates in weak and strong approximations. The Poisson–Alekseev–Gröbner formula provides an explicit stochastic error representation between solutions of SDEs differing in coefficients, under minimal regularity [2510.05300].

### Sensitivity (Greeks) in Financial Mathematics

Poissonian Malliavin calculus, via Skorokhod integration and explicit chaos expansions, yields efficient representations for Greeks in mean-field jump-diffusion models and path-dependent discontinuous payoff functions, paralleling the Bismut–Elworthy–Li formula but for pure jump processes [2502.00793]. Formulae for Delta involve the computation of predictable Skorohod weights solving integration-by-parts equations, with theoretical and computational advantages over naive finite-difference estimators.

### Poisson–Malliavin–Stein Program for Gaussian Approximation

Key results connect the Malliavin covariance \(\langle D F, -D L^{-1} F \rangle\) to optimal bounds in normal approximation for functionals of Poisson process, Poisson U-statistics, and functionals of Hawkes processes [1111.2140, 2104.01583, 2109.07749, 2110.04877]. The multidimensional and functional functional CLTs employ explicit contraction/kernels–based bounds and are central in stochastic geometry, random graphs, and non-diffusive stochastic models.

## 6. Absolute Continuity, Malliavin-Sobolev Spaces, and Regularity

Sharp conditions for absolute continuity of the law of Poisson functionals arise naturally via Malliavin calculus. The criterion
\[
\int_0^T\int_{\R_0} |D_{t,x} F|^2 \nu(dx) dt > 0 \quad \text{a.s.}
\]
on an event ensures absolute continuity of the law of \(F\) restricted to that event [1210.1156]. Appropriate weighted Lebesgue space characterizations of differentiability and fractional differentiability also follow, with the key result that for \(\mathcal{F}_A\)-measurable \(F\),
\[
F \in D^{1,p} \iff F (N(A)+1)^{1/p} \in L^p
\]
and for fractional spaces via real interpolation,
\[
F \in (L^p, D^{1,p})_{\theta,p} \iff F (N(A)+1)^\theta \in L^p
\]
demonstrating the explicit relationship between jump-count and regularity on Poisson spaces [1605.07413].

## 7. Comparison: Poisson, Wiener, and Mixed Wiener–Poisson Calculi

The Poissonian Malliavin calculus is part of a general unified stochastic calculus of variations, admitting Gaussian (Wiener), Poissonian (jump), and mixed Wiener–Poisson frameworks [1712.04098, 2502.00793]. In the Wiener case, the chain rule is global and the chaos is Hermitian. In the Poissonian setting, extra compensation terms and different commutation/measurability properties arise, and the semigroup structure is more intricate. Mixed spaces combine both differentials, support joint limit theorems, and exhibit error rates inherited from both the Gaussian and Poissonian cases.

---

**References**:  
[1712.04098] Normal Convergence Using Malliavin Calculus With Applications and Examples  
[1703.07259] Poisson Malliavin calculus in Hilbert space with an application to SPDE  
[1605.07413] A note on Malliavin smoothness on the Lévy space  
[1210.1156] Local Malliavin Calculus for Lévy Processes and Applications  
[1111.2140] Malliavin-Stein method for multi-dimensional U-statistics of Poisson point processes  
[2110.04877] Functional Gaussian approximations on Hilbert-Poisson spaces  
[1910.10181] Stable limit theorems on the Poisson space  
[2502.00793] Sensitivity Analysis for Mean-Field SDEs With Jump By Malliavin Calculus: Chaos Expansion Approach  
[2510.05300] A Poisson-Alekseev-Gröbner formula through Malliavin calculus for Poisson random integrals  
[2104.01583] The Malliavin-Stein method for Hawkes functionals  
[2109.07749] Malliavin-Stein method for the multivariate compound Hawkes process

Source: https://www.emergentmind.com/topics/poissonian-malliavin-calculus