Malliavin calculus

Last updated December 14, 2024

In probability theory and related fields, Malliavin calculus is a set of mathematical techniques and ideas that extend the mathematical field of calculus of variations from deterministic functions to stochastic processes. In particular, it allows the computation of derivatives of random variables. Malliavin calculus is also called the stochastic calculus of variations. P. Malliavin first initiated the calculus on infinite dimensional space. Then, the significant contributors such as S. Kusuoka, D. Stroock, J-M. Bismut, Shinzo Watanabe, I. Shigekawa, and so on finally completed the foundations.

Malliavin calculus is named after Paul Malliavin whose ideas led to a proof that Hörmander's condition implies the existence and smoothness of a density for the solution of a stochastic differential equation; Hörmander's original proof was based on the theory of partial differential equations. The calculus has been applied to stochastic partial differential equations as well.

The calculus allows integration by parts with random variables; this operation is used in mathematical finance to compute the sensitivities of financial derivatives. The calculus has applications in, for example, stochastic filtering.

Overview and history

Malliavin introduced Malliavin calculus to provide a stochastic proof that Hörmander's condition implies the existence of a density for the solution of a stochastic differential equation; Hörmander's original proof was based on the theory of partial differential equations. His calculus enabled Malliavin to prove regularity bounds for the solution's density. The calculus has been applied to stochastic partial differential equations.

Gaussian probability space

Consider a Wiener functional $F$ (a functional from the classical Wiener space) and consider the task of finding a derivative for it. The natural idea would be to use the Gateaux derivative

D_{g}F:=\left.{\frac {d}{d\tau }}F[f+\tau g]\right|_{\tau =0},

however this does not always exist. Therefore it does make sense to find a new differential calculus for such spaces by limiting the directions.

The toy model of Malliavin calculus is an irreducible Gaussian probability space $X=(\Omega ,{\mathcal {F}},P,{\mathcal {H}})$ . This is a (complete) probability space $(\Omega ,{\mathcal {F}},P)$ together with a closed subspace ${\mathcal {H}}\subset L^{2}(\Omega ,{\mathcal {F}},P)$ such that all $H\in {\mathcal {H}}$ are mean zero Gaussian variables and ${\mathcal {F}}=\sigma (H:H\in {\mathcal {H}})$ . If one chooses a basis for ${\mathcal {H}}$ then one calls $X$ a numerical model. On the other hand, for any separable Hilbert space ${\mathcal {G}}$ exists a canonical irreducible Gaussian probability space $\operatorname {Seg} ({\mathcal {G}})$ named the Segal model (named after Irving Segal) having ${\mathcal {G}}$ as its Gaussian subspace. In this case for a $g\in {\mathcal {G}}$ one notates the associated random variable in $\operatorname {Seg} ({\mathcal {G}})$ as $W(g)$ .

Properties of a Gaussian probability space that do not depend on the particular choice of basis are called intrinsic and such that do depend on the choice extrensic.^[1] We denote the countably infinite product of real spaces as $\mathbb {R} ^{\mathbb {N} }=\prod \limits _{i=1}^{\infty }\mathbb {R}$ .

Recall the modern version of the Cameron-Martin theorem

Consider a locally convex vector space $E$ with a cylindrical Gaussian measure $\gamma$ on it. For an element in the topological dual $f\in E'$ define the distance to the mean

t_{\gamma }(f):=f-\int _{E}f(x)\gamma (\mathrm {d} x),

which is a map $t_{\gamma }\colon E'\to L^{2}(E,\gamma )$ , and denote the closure in $L^{2}(E,\gamma )$ as

E_{\gamma }':=\operatorname {clos} \left\{t_{\gamma }(f)\colon \ f\in E'\right\}

Let $\gamma _{m}:=\gamma (\cdot -m)$ denote the translation by $m\in E$ . Then $E_{\gamma }'$ respectively the covariance operator $R_{\gamma }:E_{\gamma }'\to (E_{\gamma }')^{*}$ on it induces a reproducing kernel Hilbert space $R$ called the Cameron-Martin space such that for any $m\in R$ there is equivalence $\gamma _{m}\sim \gamma$ .^[2]

In fact one can use here the Feldman–Hájek theorem to find that for any other $h\not \in R$ such measure would be singular.

Let $\gamma$ be the canonical Gaussian measure, by transferring the Cameron-Martin theorem from $(\mathbb {R} ^{\mathbb {N} },{\mathcal {B}}(\mathbb {R} ^{\mathbb {N} }),\gamma ^{\mathbb {N} }=\otimes _{n\in \mathbb {N} }\gamma )$ into a numerical model $X$ , the additive group of ${\mathcal {H}}$ will define a quasi-automorphism group on $\Omega$ . A construction can be done as follows: choose an orthonormal basis in ${\mathcal {H}}$ , let $\tau _{\alpha }(x)=x+\alpha$ denote the translation on $\mathbb {R} ^{\mathbb {N} }$ by $\alpha$ , denote the map into the Cameron-Martin space by $j:{\mathcal {H}}\to \ell ^{2}$ , denote

L^{\infty -0}(\Omega ,{\mathcal {F}},P)=\bigcap \limits _{p<\infty }L^{p}(\Omega ,{\mathcal {F}},P)\quad

and

\quad q:L^{\infty -0}(\mathbb {R} ^{\mathbb {N} },{\mathcal {B}}(\mathbb {R} ^{\mathbb {N} }),\gamma ^{\mathbb {N} })\to L^{\infty -0}(\Omega ,{\mathcal {F}},P),

we get a canonical representation of the additive group ${\displaystyle \rho$ acting on the endomorphisms by defining

\rho (h)=q\circ \tau _{j(h)}\circ q^{-1}.

One can show that the action of $\rho$ is extrinsic meaning it does not depend on the choice of basis for ${\mathcal {H}}$ , further $\rho (h+h')=\rho (h)\rho (h')$ for $h,h'\in {\mathcal {H}}$ and for the infinitesimal generator of $(\rho (h))_{h}$ that

\lim \limits _{\varepsilon \to 0}{\frac {\rho (\varepsilon h)-I}{\varepsilon }}=M_{h}

where $I$ is the identity operator and $M_{h}$ denotes the multiplication operator by the random variable $h\in {\mathcal {H}}$ (acting on the endomorphisms). In the case of an arbitrary Hilbert space ${\mathcal {G}}$ and the Segal model $\operatorname {Seg} ({\mathcal {G}})$ one has $j:{\mathcal {G}}\to \ell ^{2}$ (and thus ${\displaystyle \rho$ . Then the limit above becomes the multiplication operator by the random variable $W(g)$ associated to $g\in {\mathcal {G}}$ .^[3]

For $F\in L^{\infty -0}(\Omega ,{\mathcal {F}},P)$ and $h\in {\mathcal {H}}$ one now defines the directional derivative

\langle DF,h\rangle =D_{h}F=\lim \limits _{\varepsilon \to 0}{\frac {\left(\rho (\varepsilon h)-I\right)F}{\varepsilon }}.

Given a Hilbert space $H$ and a Segal model $\operatorname {Seg} (H)$ with its Gaussian space ${\mathcal {H}}=\{W(h):h\in H\}$ . One can now deduce for $F\in L^{\infty -0}(\Omega ,{\mathcal {F}},P)$ the integration by parts formula

\mathbb {E} [D_{h}F]=\mathbb {E} [M_{W(h)}F]=\mathbb {E} [W(h)F]

.^[4]

Invariance principle

The usual invariance principle for Lebesgue integration over the whole real line is that, for any real number ε and integrable function f, the following holds

\int _{-\infty }^{\infty }f(x)\,d\lambda (x)=\int _{-\infty }^{\infty }f(x+\varepsilon )\,d\lambda (x)

and hence

\int _{-\infty }^{\infty }f'(x)\,d\lambda (x)=0.

This can be used to derive the integration by parts formula since, setting f = gh, it implies

0=\int _{-\infty }^{\infty }f'\,d\lambda =\int _{-\infty }^{\infty }(gh)'\,d\lambda =\int _{-\infty }^{\infty }gh'\,d\lambda +\int _{-\infty }^{\infty }g'h\,d\lambda .

A similar idea can be applied in stochastic analysis for the differentiation along a Cameron-Martin-Girsanov direction. Indeed, let $h_{s}$ be a square-integrable predictable process and set

\varphi (t)=\int _{0}^{t}h_{s}\,ds.

If $X$ is a Wiener process, the Girsanov theorem then yields the following analogue of the invariance principle:

E(F(X+\varepsilon \varphi ))=E\left[F(X)\exp \left(\varepsilon \int _{0}^{1}h_{s}\,dX_{s}-{\frac {1}{2}}\varepsilon ^{2}\int _{0}^{1}h_{s}^{2}\,ds\right)\right].

Differentiating with respect to ε on both sides and evaluating at ε=0, one obtains the following integration by parts formula:

E(\langle DF(X),\varphi \rangle )=E{\Bigl [}F(X)\int _{0}^{1}h_{s}\,dX_{s}{\Bigr ]}.

Here, the left-hand side is the Malliavin derivative of the random variable $F$ in the direction $\varphi$ and the integral appearing on the right hand side should be interpreted as an Itô integral.

Clark–Ocone formula

One of the most useful results from Malliavin calculus is the Clark–Ocone theorem, which allows the process in the martingale representation theorem to be identified explicitly. A simplified version of this theorem is as follows:

Consider the standard Wiener measure on the canonical space $C[0,1]$ , equipped with its canonical filtration. For $F:C[0,1]\to \mathbb {R}$ satisfying $E(F(X)^{2})<\infty$ which is Lipschitz and such that F has a strong derivative kernel, in the sense that for $\varphi$ in C[0,1]

\lim _{\varepsilon \to 0}{\frac {1}{\varepsilon }}(F(X+\varepsilon \varphi )-F(X))=\int _{0}^{1}F'(X,dt)\varphi (t)\ \mathrm {a.e.} \ X

then

F(X)=E(F(X))+\int _{0}^{1}H_{t}\,dX_{t},

where H is the previsible projection of F'(x, (t,1]) which may be viewed as the derivative of the function F with respect to a suitable parallel shift of the process X over the portion (t,1] of its domain.

This may be more concisely expressed by

F(X)=E(F(X))+\int _{0}^{1}E(D_{t}F\mid {\mathcal {F}}_{t})\,dX_{t}.

Much of the work in the formal development of the Malliavin calculus involves extending this result to the largest possible class of functionals F by replacing the derivative kernel used above by the "Malliavin derivative" denoted $D_{t}$ in the above statement of the result. ^{[ citation needed ]}

Skorokhod integral

The Skorokhod integral operator which is conventionally denoted δ is defined as the adjoint of the Malliavin derivative in the white noise case when the Hilbert space is an $L^{2}$ space, thus for u in the domain of the operator which is a subset of $L^{2}([0,\infty )\times \Omega )$ , for F in the domain of the Malliavin derivative, we require

E(\langle DF,u\rangle )=E(F\delta (u)),

where the inner product is that on $L^{2}[0,\infty )$ viz

\langle f,g\rangle =\int _{0}^{\infty }f(s)g(s)\,ds.

The existence of this adjoint follows from the Riesz representation theorem for linear operators on Hilbert spaces.

It can be shown that if u is adapted then

\delta (u)=\int _{0}^{\infty }u_{t}\,dW_{t},

where the integral is to be understood in the Itô sense. Thus this provides a method of extending the Itô integral to non adapted integrands.

Applications

The calculus allows integration by parts with random variables; this operation is used in mathematical finance to compute the sensitivities of financial derivatives. The calculus has applications for example in stochastic filtering.

Related Research Articles

In fractional calculus, an area of mathematical analysis, the differintegral is a combined differentiation/integration operator. Applied to a function ƒ, the q-differintegral of f, here denoted by

In vector calculus, Green's theorem relates a line integral around a simple closed curve $C$ to a double integral over the plane region $D$ bounded by $C$ . It is the two-dimensional special case of Stokes' theorem. In one dimension, it is equivalent to the fundamental theorem of calculus. In three dimensions, it is equivalent to the divergence theorem.

In mathematics, a Lie algebroid is a vector bundle $together with a Lie bracket on its space of sections and a vector bundle morphism, satisfying a Leibniz rule. A Lie algebroid can thus be thought of as a "many-object generalisation" of a Lie algebra.$

In mathematics, and especially differential geometry and gauge theory, a connection on a fiber bundle is a device that defines a notion of parallel transport on the bundle; that is, a way to "connect" or identify fibers over nearby points. The most common case is that of a linear connection on a vector bundle, for which the notion of parallel transport must be linear. A linear connection is equivalently specified by a covariant derivative, an operator that differentiates sections of the bundle along tangent directions in the base manifold, in such a way that parallel sections have derivative zero. Linear connections generalize, to arbitrary vector bundles, the Levi-Civita connection on the tangent bundle of a pseudo-Riemannian manifold, which gives a standard way to differentiate vector fields. Nonlinear connections generalize this concept to bundles whose fibers are not necessarily linear.

In the theory of stochastic processes, the Karhunen–Loève theorem, also known as the Kosambi–Karhunen–Loève theorem states that a stochastic process can be represented as an infinite linear combination of orthogonal functions, analogous to a Fourier series representation of a function on a bounded interval. The transformation is also known as Hotelling transform and eigenvector transform, and is closely related to principal component analysis (PCA) technique widely used in image processing and in data analysis in many fields.

In theoretical physics and mathematics, a Wess–Zumino–Witten (WZW) model, also called a Wess–Zumino–Novikov–Witten model, is a type of two-dimensional conformal field theory named after Julius Wess, Bruno Zumino, Sergei Novikov and Edward Witten. A WZW model is associated to a Lie group, and its symmetry algebra is the affine Lie algebra built from the corresponding Lie algebra. By extension, the name WZW model is sometimes used for any conformal field theory whose symmetry algebra is an affine Lie algebra.

In mathematics, in particular in algebraic geometry and differential geometry, Dolbeault cohomology (named after Pierre Dolbeault) is an analog of de Rham cohomology for complex manifolds. Let M be a complex manifold. Then the Dolbeault cohomology groups $depend on a pair of integers p and q and are realized as a subquotient of the space of complex differential forms of degree (p, q).$

In mathematics, the Clark–Ocone theorem is a theorem of stochastic analysis. It expresses the value of some function F defined on the classical Wiener space of continuous paths starting at the origin as the sum of its mean value and an Itô integral with respect to that path. It is named after the contributions of mathematicians J.M.C. Clark (1970), Daniel Ocone (1984) and U.G. Haussmann (1978).

<span class="mw-page-title-main">Classical Wiener space</span>

In mathematics, classical Wiener space is the collection of all continuous functions on a given domain, taking values in a metric space. Classical Wiener space is useful in the study of stochastic processes whose sample paths are continuous functions. It is named after the American mathematician Norbert Wiener.

In mathematics, Weyl's lemma, named after Hermann Weyl, states that every weak solution of Laplace's equation is a smooth solution. This contrasts with the wave equation, for example, which has weak solutions that are not smooth solutions. Weyl's lemma is a special case of elliptic or hypoelliptic regularity.

In mathematics, uniform integrability is an important concept in real analysis, functional analysis and measure theory, and plays a vital role in the theory of martingales.

In mathematics, the Möbius energy of a knot is a particular knot energy, i.e., a functional on the space of knots. It was discovered by Jun O'Hara, who demonstrated that the energy blows up as the knot's strands get close to one another. This is a useful property because it prevents self-intersection and ensures the result under gradient descent is of the same knot type.

In mathematics, the Skorokhod integral, also named Hitsuda–Skorokhod integral, often denoted $, is an operator of great importance in the theory of stochastic processes. It is named after the Ukrainian mathematician Anatoliy Skorokhod and Japanese mathematician Masuyuki Hitsuda. Part of its importance is that it unifies several concepts:$

In mathematics, the Kodaira–Spencer map, introduced by Kunihiko Kodaira and Donald C. Spencer, is a map associated to a deformation of a scheme or complex manifold X, taking a tangent space of a point of the deformation space to the first cohomology group of the sheaf of vector fields on X.

In stochastic analysis, a rough path is a generalization of the notion of smooth path allowing to construct a robust solution theory for controlled differential equations driven by classically irregular signals, for example a Wiener process. The theory was developed in the 1990s by Terry Lyons. Several accounts of the theory are available.

Quantum stochastic calculus is a generalization of stochastic calculus to noncommuting variables. The tools provided by quantum stochastic calculus are of great use for modeling the random evolution of systems undergoing measurement, as in quantum trajectories. Just as the Lindblad master equation provides a quantum generalization to the Fokker–Planck equation, quantum stochastic calculus allows for the derivation of quantum stochastic differential equations (QSDE) that are analogous to classical Langevin equations.

The theory of causal fermion systems is an approach to describe fundamental physics. It provides a unification of the weak, the strong and the electromagnetic forces with gravity at the level of classical field theory. Moreover, it gives quantum mechanics as a limiting case and has revealed close connections to quantum field theory. Therefore, it is a candidate for a unified physical theory. Instead of introducing physical objects on a preexisting spacetime manifold, the general concept is to derive spacetime as well as all the objects therein as secondary objects from the structures of an underlying causal fermion system. This concept also makes it possible to generalize notions of differential geometry to the non-smooth setting. In particular, one can describe situations when spacetime no longer has a manifold structure on the microscopic scale. As a result, the theory of causal fermion systems is a proposal for quantum geometry and an approach to quantum gravity.

In mathematics, Katugampola fractional operators are integral operators that generalize the Riemann–Liouville and the Hadamard fractional operators into a unique form. The Katugampola fractional integral generalizes both the Riemann–Liouville fractional integral and the Hadamard fractional integral into a single form and It is also closely related to the Erdelyi–Kober operator that generalizes the Riemann–Liouville fractional integral. Katugampola fractional derivative has been defined using the Katugampola fractional integral and as with any other fractional differential operator, it also extends the possibility of taking real number powers or complex number powers of the integral and differential operators.

In probability theory, a branch of mathematics, white noise analysis, otherwise known as Hida calculus, is a framework for infinite-dimensional and stochastic calculus, based on the Gaussian white noise probability space, to be compared with Malliavin calculus based on the Wiener process. It was initiated by Takeyuki Hida in his 1975 Carleton Mathematical Lecture Notes.

The variational multiscale method (VMS) is a technique used for deriving models and numerical methods for multiscale phenomena. The VMS framework has been mainly applied to design stabilized finite element methods in which stability of the standard Galerkin method is not ensured both in terms of singular perturbation and of compatibility conditions with the finite element spaces.

References

↑ Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. pp. 4–15. ISBN 3-540-57024-1.
↑ Bogachev, Vladimir (1998). Gaussian Measures. Rhode Island: American Mathematical Society.
↑ Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. pp. 20–22. ISBN 3-540-57024-1.
↑ Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. p. 36. ISBN 3-540-57024-1.

Kusuoka, S. and Stroock, D. (1981) "Applications of Malliavin Calculus I", Stochastic Analysis, Proceedings Taniguchi International Symposium Katata and Kyoto 1982, pp 271–306
Kusuoka, S. and Stroock, D. (1985) "Applications of Malliavin Calculus II", J. Faculty Sci. Uni. Tokyo Sect. 1A Math., 32 pp 1–76
Kusuoka, S. and Stroock, D. (1987) "Applications of Malliavin Calculus III", J. Faculty Sci. Univ. Tokyo Sect. 1A Math., 34 pp 391–442
Malliavin, Paul and Thalmaier, Anton. Stochastic Calculus of Variations in Mathematical Finance, Springer 2005, ISBN 3-540-43431-3
Nualart, David (2006). The Malliavin calculus and related topics (Second ed.). Springer-Verlag. ISBN 978-3-540-28328-7.
Bell, Denis. (2007) The Malliavin Calculus, Dover. ISBN 0-486-44994-7; ebook
Sanz-Solé, Marta (2005) Malliavin Calculus, with applications to stochastic partial differential equations. EPFL Press, distributed by CRC Press, Taylor & Francis Group.
Schiller, Alex (2009) Malliavin Calculus for Monte Carlo Simulation with Financial Applications. Thesis, Department of Mathematics, Princeton University
Øksendal, Bernt K.(1997) An Introduction To Malliavin Calculus With Applications To Economics. Lecture Notes, Dept. of Mathematics, University of Oslo (Zip file containing Thesis and addendum)
Di Nunno, Giulia, Øksendal, Bernt, Proske, Frank (2009) "Malliavin Calculus for Lévy Processes with Applications to Finance", Universitext, Springer. ISBN 978-3-540-78571-2

External links

Quotations related to Malliavin calculus at Wikiquote
Friz, Peter K. (2005-04-10). "An Introduction to Malliavin Calculus" (PDF). Archived from the original (PDF) on 2007-04-17. Retrieved 2007-07-23. Lecture Notes, 43 pages

Zhang, H. (2004-11-11). "The Malliavin Calculus" (PDF). Retrieved 2004-11-11. Thesis, 100 pages

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. pp. 4–15. ISBN 3-540-57024-1.

[2] Bogachev, Vladimir (1998). Gaussian Measures. Rhode Island: American Mathematical Society.

[3] Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. pp. 20–22. ISBN 3-540-57024-1.

[4] Malliavin, Paul (1997). Stochastic Analysis. Grundlehren der mathematischen Wissenschaften. Berlin, Heidelberg: Springer. p. 36. ISBN 3-540-57024-1.

[1]

[2]

[3]

[4]