Malliavin calculus

Last updated September 29, 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.

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 having ${\mathcal {G}}$ as its Gaussian subspace. 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 and denote the closure of the topological dual to the mean

E_{\gamma }^{*}:=\operatorname {clos} \left\{f-\int _{E}f(x)\gamma (\mathrm {d} x),\ f\in E^{*}\right\}.

Let $\gamma _{m}:=\gamma (\cdot -m)$ denote the translation by $m\in E$ . Then $E_{\gamma }^{*}$ 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 H$ such measure would be singular.

We denote the countably infinite product of real spaces as $\mathbb {R} ^{\mathbb {N} }=\prod \limits _{i=1}^{\infty }\mathbb {R}$ .

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 on $\Omega$ associated to $h\in {\mathcal {H}}$ (acting on the endomorphisms).^[3]

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

<span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

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, 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 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 the mathematical field of dynamical systems, a random dynamical system is a dynamical system in which the equations of motion have an element of randomness to them. Random dynamical systems are characterized by a state space S, a set of maps $from S into itself that can be thought of as the set of all possible equations of motion, and a probability distribution Q on the set that represents the random choice of map. Motion in a random dynamical system can be informally thought of as a state evolving according to a succession of maps randomly chosen according to the distribution Q .$

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, 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:$

The narrow escape problem is a ubiquitous problem in biology, biophysics and cellular biology.

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.

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.

In stochastic calculus, the Ogawa integral, also called the non-causal stochastic integral, is a stochastic integral for non-adapted processes as integrands. The corresponding calculus is called non-causal calculus which distinguishes it from the anticipating calculus of the Skorokhod integral. The term causality refers to the adaptation to the natural filtration of the integrator.

In probability theory particularly in the Malliavin calculus, a Gaussian probability space is a probability space together with a Hilbert space of mean zero, real-valued Gaussian random variables. Important examples include the classical or abstract Wiener space with some suitable collection of Gaussian random variables.

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.

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.

[1]

[2]

[3]