---
title: Malliavin Calculus & BEL Formula
url: https://www.emergentmind.com/topics/malliavin-calculus-and-bel-formula
type: topic
---

# Malliavin Calculus & BEL Formula

Malliavin calculus is an infinite-dimensional differential calculus on Wiener space that enables rigorous analysis of derivatives of random variables defined on probability spaces generated by Brownian motion. One of its central achievements is the derivation of integration-by-parts and sensitivity formulas for functionals of stochastic differential equations (SDEs), both regular and degenerate. Central to probabilistic analysis and quantitative finance is the Bismut–Elworthy–Li (BEL) formula, which uses Malliavin calculus to provide representations of derivatives (Greeks) of expectations for SDEs, semigroups, and associated PDE solutions—often in the form of unbiased Monte Carlo estimators. The BEL formula is a fundamental tool in deriving gradient estimates, Harnack inequalities, sensitivity representations for SPDEs, mean-field SDEs, hypoelliptic diffusions, systems with jumps, rough volatility models, and diffusions with killing or boundaries.

## 1. Mathematical Setting of Malliavin Calculus and BEL Formula

Malliavin calculus constructs a differential structure on functionals \(F = f(W)\) of Brownian motion \(W\), with the Malliavin derivative operator \(D\) acting as an infinite-dimensional gradient and its adjoint, the divergence operator \(\delta\) (Skorokhod integral), generalizing stochastic integration. On an SDE
\[
dX_t = b(X_t)dt + \sigma(X_t)dW_t,\quad X_0 = x
\]
under appropriate regularity and nondegeneracy (\(\sigma\sigma^\top \geq \varepsilon I\)), the Malliavin derivative captures how changes in the underlying Brownian path affect \(X_t\), facilitating a probabilistic version of classical differentiation under the expectation.

The BEL formula expresses the gradient of the semigroup \(P_tf(x) = \mathbb{E} f(X_t^x)\) as
\[
\nabla_x P_t f(x) = \mathbb{E}\left[ f(X_t^x) \int_0^t \left(\sigma^{-1}(X_s) Y_s\right)^\top dW_s \right]
\]
where \(Y_s = \partial_x X_s^x\) is the Jacobian (the first-variation) process, and the stochastic integral is interpreted as a Skorokhod or Itô integral depending on adaptedness [2106.04901].

If the dynamics are degenerate or hypoelliptic,
\[
\begin{cases}
dX_t^{(1)} = Z^{(1)}(X_t)dt \\
dX_t^{(2)} = Z^{(2)}(X_t)dt + \sigma dB_t
\end{cases}
\]
the representation is more intricate and requires constructing a suitable control process \(h\) via a linear control problem. The resulting representation is again
\[
\nabla_v P_T f(x) = \mathbb{E}\left[ f(X_T^x) \delta(h(v)) \right]
\]
with \(\delta(h)\) the divergence operator, and the path \(h\) constructed to ensure \(D_h X_T^x = \nabla_v X_T^x\). This underpins explicit gradient formulas, dimension-free Harnack inequalities, and integrability estimates [1107.0096].

## 2. Malliavin Calculus: Core Constructions and Integration by Parts

Malliavin’s approach introduces the Malliavin derivative \(D\), which for a smooth random variable \(F = f(W(h_1),...,W(h_n))\) has
\[
D F = \sum_{i=1}^n \frac{\partial f}{\partial x_i}(W(h_1),...,W(h_n)) h_i
\]
and the divergence operator \(\delta\) satisfies the duality
\[
\mathbb{E}\left[\langle D F, u\rangle_H\right] = \mathbb{E}\left[ F \delta(u) \right]
\]
for suitable \(u\) [2312.00405].

The BEL formula arises directly from this integration-by-parts principle, transferring derivatives from the payoff/composite functional onto a stochastic “weight” process independent of the payoff:
\[
\mathbb{E}\left[\nabla f(X_T)\cdot Y_T\right]
= \mathbb{E}\left[f(X_T)\, \delta(U) \right]
\]
with \(U\) constructed from the inverse diffusion and Jacobian. In nondegenerate cases \(U = \sigma^{-1}(X)Y\), but in degenerate chains, more complex control-theoretic constructions for \(h\) are essential [1107.0096].

Discontinuous payoffs are handled by localizing or approximating the payoff, yielding robust estimators without loss of convergence, as in the computation of Greeks for spread options or digital payoffs [2106.04901, 1510.06961, 2312.00405].

## 3. Generalizations: Degenerate, Mean-Field, and Infinite-Dimensional Systems

### 3.1 Degenerate and Hypoelliptic Diffusions

For SDEs with degenerate noise, the BEL formula persists if the missing directions can be recovered by the underlying Lie bracket structure (Hörmander’s condition) or explicit nondegeneracy of a control matrix, as in
\[
\nabla^{(2)} Z^{(1)} = B_0 + B(x),\quad \text{$B_0$ full rank}
\]
Hamiltonian structures and linearized ODE/PDE flows (the control problem for \(g_t\)) ensure the solvability of the Malliavin weight [1107.0096].

### 3.2 Mean-Field SDEs

If the coefficients depend on the law of the solution (mean-field, McKean–Vlasov systems), the first-variation process incorporates flow terms from distribution-dependent drift/diffusion. The adapted BEL formula is
\[
\nabla_x v(x) = \mathbb{E}\Big[ \Phi(X_T^x) \int_0^T a(s)[\sigma^{-1}(s,X_s^x,\pi_s^x)Y_su(T)]^* \delta W_s \Big]^*
\]
where \(u(T)\) resolves mean-field contributions, and the Skorokhod integral may require correction to yield an adapted Itô form [1510.06961].

### 3.3 SPDEs and Infinite Dimensions

The Malliavin framework and the BEL formula extend to infinite-dimensional systems (e.g., SPDEs with multiplicative or subordinate noise), under appropriate trace-class and invertibility assumptions for the noise operator. For a semigroup \(P_T\) associated to the solution of an SPDE,
\[
dX_t = A X_t dt + b_t(X_t) dt + \sigma_t(X_t) dW_t
\]
the integration by parts yields
\[
\nabla_v P_T f(x) = \mathbb{E}\left[ f(X_T(x)) M_T(x;v) \right]
\]
with \(M_T\) explicitly decomposed into stochastic integrals and trace terms involving the Jacobian and the noise. In finite-dimensional reductions, this recovers classical BEL [1610.02778, 1601.01733].

## 4. Applications: Sensitivities, Greeks, and Quantitative Bounds

The primary application of the Malliavin–BEL formalism is the computation of sensitivities—Greeks—in financial mathematics, where the formula yields expressions for delta, gamma, vega, etc., in models with multiple assets, stochastic and rough volatility, mean-field effects, and more. For the spread option in two dimensions with correlated geometric Brownian motion,
\[
\Delta_1 = \frac{1}{x_1\sigma_1\sqrt{1-\rho^2}T}(W_1(T) - \rho W_2(T))
\]
so that
\[
\partial_{x_1} u = e^{-rT} \mathbb{E}\left[ \phi(S_1(T),S_2(T))\,\Delta_1 \right]
\]
Monte Carlo estimation proceeds by simulating the underlying dynamics, accumulating the stochastic integrals (Malliavin weights), and estimating the expectation for any payoff function, with the weight process fully independent of the payoff [2106.04901].

Analogous approaches yield Greek representations for rough Volterra models, including rough Stein–Stein, SABR, and Bergomi models. Under mild regularity, Malliavin-weighted Monte Carlo demonstrates efficient and stable numerics even for discontinuous payoffs or rough volatility [2312.00405].

In the elliptic and hypoelliptic settings, the BEL formula underlies the derivation of Harnack inequalities and explicit Gaussian and non-Gaussian upper/lower bounds on the heat kernel, connecting Malliavin calculus to analysis on diffusion semigroups and degenerate PDEs [1107.0096, 1610.02778, 1601.01733].

Tables of BEL Representations in Distinct Settings:

| Setting               | BEL/Malliavin Formula Structure                                      | Reference          |
|-----------------------|---------------------------------------------------------------------|--------------------|
| Elliptic SDE          | $\mathbb{E}[f(X_T)\int_0^T (\sigma^{-1}Y_s)dW_s]$                  | [2106.04901]       |
| Degenerate/Hypoelliptic | $\mathbb{E}[f(X_T)\,\delta(h)]$, $h$ via control problem            | [1107.0096]        |
| Mean-field SDE        | $\mathbb{E}[\Phi(X_T)\int_0^T K_s u(T)^\ast \delta W_s]$            | [1510.06961]       |
| SPDEs                 | $\mathbb{E}[f(X_T(x)) M_T(x;v)]$, $M_T$ includes traces of Jacobians | [1610.02778]       |
| SDEs with jumps       | $\mathbb{E}[f(X_T) M_T^v]$, $M_T^v$ via time-changed Brownian motion | [1601.01733]       |
| Rough volatility      | $\mathbb{E}[f(S_T)\delta(\frac{\partial_\theta S_T}{\int D_s S_T ds})]$ | [2312.00405]   |
| Killed processes      | Unbiased Markov-chain based weighted sums (reflection, discrete IBP) | [1908.04550]       |

## 5. Extensions: Jumps, Killing, Boundaries, and Rough Dynamics

The BEL construction admits robust extensions:

- **Jump processes**: For SDEs/SPDEs with jump-type noise (subordinators, Lévy processes), finite-jump approximations yield integration-by-parts and BEL formulas, with weight processes reflecting the time-changed nature of the noise [1601.01733].
- **Killed processes/Boundaries**: Approximations based on Poisson time-embedding and novel discrete Malliavin constructions permit BEL-type formulas for diffusions killed at boundaries, enabling unbiased simulation and unbiased estimation of boundary-corrected derivatives [1908.04550].
- **Rough and Volterra dynamics**: In rough volatility models, the BEL formula generalizes to non-semimartingale settings, provided the Volterra kernel admits a suitable regularization and the required Malliavin differentiability holds; explicit Malliavin-weighted formulas for Greeks are derived and analyzed for numerical convergence [2312.00405].

## 6. Implications, Limitations, and Further Directions

The Malliavin–BEL framework offers several advantages:

- Malliavin weights are independent of the payoff; once realized, the same trajectory provides unbiased estimators for arbitrary payoffs.
- Sensitivity estimations for discontinuous payoffs remain robust—avoiding the bias and instability of finite-difference approximations.
- The formalism provides explicit dimension-free gradient bounds and Harnack inequalities for non-degenerate and even degenerate semigroups, with Gaussian-type or stable-type heat kernel estimates depending on the noise structure [1107.0096, 1610.02778, 1601.01733].

However, the approach relies on invertibility or nondegeneracy conditions on the diffusion/coefficient matrices (or their control-theoretic generalizations). Heavy degeneracy may render the Malliavin matrix singular, in which case classical BEL formulas can fail. For high-variance or rough paths (e.g., under exponential waiting times in killed chains), variance explosion can occur, necessitating advanced renewal schemes or variance-reduction techniques [1908.04550].

Future directions include further generalizations to Lévy-driven SDEs, SPDEs with non-trivial boundary interactions, stochastic control with mean-field interactions, and rougher structures with possibly distributional coefficients.

The extensive development of Malliavin calculus for infinite-dimensional systems and singular noise (including rough paths and jumps) continues to yield both theoretical and computational advances, as evidenced especially in the domain of quantitative finance and modern stochastic analysis [1107.0096, 2106.04901, 1510.06961, 1610.02778, 1601.01733, 2312.00405, 1908.04550].

Source: https://www.emergentmind.com/topics/malliavin-calculus-and-bel-formula