An Introduction to Malliavin Calculus

Total Page:16

File Type:pdf, Size:1020Kb

An Introduction to Malliavin Calculus An Introduction to Malliavin Calculus Lecture Notes SummerTerm 2013 by Markus Kunze Contents Chapter 1. Stochastic Calculus 1 1.1. The Wiener Chaos Decomposition 1 1.2. The Malliavin Derivative 6 1.3. The Divergence Operator 15 1.4. The Ornstein-Uhlenbeck Semigroup 19 1.5. Multiple Wiener Integrals 20 1.6. Stochastic Calculus in the White Noise Case 25 1.7. It^o'sIntegral and the Clark-Ocone Formula 30 Chapter 2. Smoothness of Probability Laws 33 2.1. Absolute Continuity 33 2.2. Smoothness of the Density 39 2.3. Stochastic Differential Equations 44 2.4. H¨ormander'sTheorem 48 Bibliography 57 iii CHAPTER 1 Stochastic Calculus The general setting for Malliavin calculus is a Gaussian probability space, i.e. a proba- bility space (Ω; Σ; P) along with a closed subspace H of L2(Ω; Σ; P) consisting of centered Gaussian random variables. It is often convenient to assume that H is isometric to another Hilbert space H, typically an L2-space over a parameter set T . 1.1. The Wiener Chaos Decomposition We recall that a real-valued random variable X, defined on a probability space (Ω; Σ; P) is called Gaussian (or, more precisely, centered Gaussian) if its characteristic function 'X , q t2 defined by 'X (t) := E exp(itX) is of the form 'X (t) = e 2 for some q ≥ 0. It is well- known that a Gaussian random variable is either constantly zero (in which case q = 0) or, else, q > 0 and the distribution has density 2 1 − x p e 2q 2πq with respect to Lebesgue measure dx. In that case, the random variable has finite moments of all orders, its mean is zero (whence it is called \centered") and its variance is q. A Gaussian random variable is called standard if it has variance 1. A family (Xi)i2I of real-valued random variables is called Gaussian family or jointly Gaussian, if for any n 2 N and any choice i1; : : : ; in of distinct indices in I, the vector Rn (Xi1 ;:::;Xin ) is a Gaussian vector. The latter means that for any α 2 , the real-valued Pn random variable k=1 αkXik is Gaussian. We now introduce isonormal Gaussian processes. Definition 1.1.1. Let H be a real, separable Hilbert space with inner product h · ; · i. An H-isonormal Gaussian process is a family W = fW (h): h 2 Hg of real-valued random variables, defined on a common probability space (Ω; Σ; P), such that W (h) is a Gaussian random variable for all h 2 H and, for h; g 2 H, we have EW (h)W (g) = hh; gi. If the Hilbert space H is clear from the context, we will also speak of isonormal Gaussian processes. Given an isonormal Gaussian process, the probability space on which the random variables are defined will be denoted by (Ω; Σ; P). We now note some properties of isonormal Gaussian processes. Proposition 1.1.2. Let H be a real, separable Hilbert space. (1) There exists an H-isonormal Gaussian process W . Now let W = fW (h): h 2 Hg be any H-isonormal Gaussian process. (2) The map h 7! W (h) is an isometry, in particular, it is linear. (3) W is a Gaussian family. Proof. (2) It suffices to prove that W is linear; in that case, it follows directly from the definition of isonormal process that W is an isometry. Thus, let h; g 2 H and λ, µ 2 R. We have EW (λh + µg) − λW (h) − µW (g)2 = EW (λh − µg)2 − 2λEW (λh + µg)W (h) − 2µEW (λh + µg)W (g) 1 2 1. STOCHASTIC CALCULUS +λ2EW (h)2 + 2λµEW (h)W (g) + µ2EW (g)2 = kλh + µgk2 − 2λhλh + µg; hi − 2µhλh + µg; gi + λ2khk2 + 2λµhg; hi + µ2kgk2 = 0 : This implies that W (λh + µg) = λW (h) + µW (g) almost surely. n (3) For n 2 N, h1; : : : ; hn 2 H and α 2 R , we have, by linearity, n n X X αkW (hk) = W αkhk k=1 k=1 which is Gaussian by assumption. (1) Since H is separable, it has a countable orthonormal basis (hk)k2N (if the basis is finite, the proof is similar). Let γ denote standard Gaussian measure on R and consider the infinite product space 1 1 1 Y O O (Ω; Σ; P) := R; B(R); γ k=1 k=1 k=1 By construction, the random variables Xn, defined by Xn((!k)k2N) := !n are independent and Gaussian. In particular, they form an orthonormal basis of their closed linear span in 2 2 L (Ω; Σ; P). We now define W : H ! L (Ω; Σ; P) to be the isometry that sends hk to Xk. Then W is an H-isonormal Gaussian process as is easy to see. Let us give some examples of isonormal Gaussian processes. Example 1.1.3. (Brownian motion) Consider the Hilbert space H = L2((0; 1); B(0; 1); λ), where λ is Lebesgue measure. By Proposition 1.1.2, there exists an H-isonormal Gaussian process W . Let Bt := W (1(0;t]). Then Bt is a centered Gaussian random variable with variance E 2 E 1 1 1 1 Bt = W ( (0;t])W ( (0;t]) = (0;t]; (0;t] = t 1 1 1 Moreover, given 0 ≤ t0 < t1 : : : < tn = t < t+s, the functions (t0;t1];:::; (tn−1;t]; (t;t+sj are 1 orthogonal in H, hence the random variables Bt1 −Bt0 = W ( (t0;t1]);:::;Bt−Btn−1 ;Bt+s−Bt are orthogonal in L2(Ω), i.e. they are uncorrelated. As these random variables are jointly Gaussian (Proposition 1.1.2), they are independent. This shows that Bt+s−Bt is independent of the σ-algebra Ft := σ(Br : r ≤ t). Thus, the process (Bt) is a Brownian motion | except for the fact that we have not proved that it has continuous paths. However, using Kolmogorov's theorem, it can be proved that Bt has a continuous modification. We refer the reader to the literature for a proof. One often writes Z T f(t) dBt := W (1(0;T )f) 0 R T and calls 0 f(t) dBt the Wiener integral of f over (0;T ). Example 1.1.4. (d-dimensional Brownian motion) Consider the Hilbert space H = L2((0;; 1); B((0; 1); λ; Rd). We denote the canonical d Basis of R by (e1; : : : ; en), i.e. ei is the vector which has a 1 at position i and 0's at all other j 1 1 d positions. If we put Bt := W ( (0;t]ej), then the vector Bt := (Bt ;:::;Bt ) is a d-dimensional j j i Brownian motion, i.e. Bt is a Brownian motion for every j and Bt and Bt are independent for j 6= i. Let us also mention, without going into details, a further example which motivates the abstract setting we consider. 1.1. THE WIENER CHAOS DECOMPOSITION 3 Example 1.1.5. (Fractional Brownian motion) A fractional Brownian motion is a Gaussian process with covariance function 1 c (t; s) = (t2H + s2H − jt − sj2H ) H 2 1 where H 2 (0; 1) is the so-called Hurst parameter. The choice H = yields c 1 (t; s) = 2 2 minft; sg which is the covariance function of Brownian motion. Let E denote the step functions on, say, (0;T ). It can be proved that there exists an inner product h · ; · iH on E such that 1 ; 1 = c (t; s) : (0;t] (0;s] H H We denote by VH the completion of E with respect to the inner product h · ; · iH . Then an VH -isonormal Gaussian process W gives rise, as above, to a fractional Brownian motion with Hurst parameter H. The range of the isonormal process W is the subspace H that was mentioned in the introduction. We now begin to study the structure of that space. This will be done using Hermite polynomials. Definition 1.1.6. For n 2 N0, the n-th Hermite polynomial Hn is defined by H0 ≡ 1 and n 2 n 2 (−1) x d − x H (x) := e 2 e 2 n n! dxn for n ≥ 1. Consider the function F (t; x) := exp(tx − t2=2). Then the Hermite polynomial are the coefficients in the power series expansion of F with respect to t. Indeed, we have x2 1 F (t; x) = exp − (x − t)2 2 2 2 1 n n 2 x X t d − (x−t) = e 2 e 2 n! dtn t=0 n=0 1 n 2 n 2 1 X n (−1) x d − z X n = t e 2 e 2 = t Hn(x) : n! dzn z=x n=0 n=0 We note that the convergence of the series is uniform on compact sets. Now some basic properties of the Hermite polynomials follow easily: Lemma 1.1.7. For n ≥ 1, we have 0 (1) Hn(x) = Hn−1(x); (2) (n + 1)Hn+1(x) = xHn(x) − Hn−1(x) n (3) Hn(−x) = (−1) Hn(x). Proof. Throughout, we set F (t; x) := exp(tx − t2=2). @ P1 n+1 (1) We have @x F (t; x) = tF (t; x) = n=0 t Hn(x). On the other hand, interchanging @ P1 n 0 summation and differentiation, @x F (t; x) = n=0 t Hn(x). Thus (1) follows by equating coefficients. @ P1 n (2) follows similarly, observing that @t F (t; x) = (x − t)F (t; x) = n=0 t xHn(x) − n+1 @ P1 n−1 t Hn(x), while, on the other hand, @t F (t; x) = n=0 nt Hn(x).
Recommended publications
  • Application of Malliavin Calculus and Wiener Chaos to Option Pricing Theory
    Application of Malliavin Calculus and Wiener Chaos to Option Pricing Theory by Eric Ben-Hamou The London School of Economics and Political Science Thesis submitted in partial fulfilment of the requirements of the degree of Doctor of Philosophy at the University of London UMI Number: U615447 All rights reserved INFORMATION TO ALL USERS The quality of this reproduction is dependent upon the quality of the copy submitted. In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if material had to be removed, a note will indicate the deletion. Dissertation Publishing UMI U615447 Published by ProQuest LLC 2014. Copyright in the Dissertation held by the Author. Microform Edition © ProQuest LLC. All rights reserved. This work is protected against unauthorized copying under Title 17, United States Code. ProQuest LLC 789 East Eisenhower Parkway P.O. Box 1346 Ann Arbor, Ml 48106-1346 weseS F 7855 % * A bstract This dissertation provides a contribution to the option pricing literature by means of some recent developments in probability theory, namely the Malliavin Calculus and the Wiener chaos theory. It concentrates on the issue of faster convergence of Monte Carlo and Quasi-Monte Carlo simulations for the Greeks, on the topic of the Asian option as well as on the approximation for convexity adjustment for fixed income derivatives. The first part presents a new method to speed up the convergence of Monte- Carlo and Quasi-Monte Carlo simulations of the Greeks by means of Malliavin weighted schemes. We extend the pioneering works of Fournie et al.
    [Show full text]
  • Malliavin Calculus in the Canonical Levy Process: White Noise Theory and Financial Applications
    Purdue University Purdue e-Pubs Open Access Dissertations Theses and Dissertations January 2015 Malliavin Calculus in the Canonical Levy Process: White Noise Theory and Financial Applications. Rolando Dangnanan Navarro Purdue University Follow this and additional works at: https://docs.lib.purdue.edu/open_access_dissertations Recommended Citation Navarro, Rolando Dangnanan, "Malliavin Calculus in the Canonical Levy Process: White Noise Theory and Financial Applications." (2015). Open Access Dissertations. 1422. https://docs.lib.purdue.edu/open_access_dissertations/1422 This document has been made available through Purdue e-Pubs, a service of the Purdue University Libraries. Please contact [email protected] for additional information. Graduate School Form 30 Updated 1/15/2015 PURDUE UNIVERSITY GRADUATE SCHOOL Thesis/Dissertation Acceptance This is to certify that the thesis/dissertation prepared By Rolando D. Navarro, Jr. Entitled Malliavin Calculus in the Canonical Levy Process: White Noise Theory and Financial Applications For the degree of Doctor of Philosophy Is approved by the final examining committee: Frederi Viens Chair Jonathon Peterson Michael Levine Jose Figueroa-Lopez To the best of my knowledge and as understood by the student in the Thesis/Dissertation Agreement, Publication Delay, and Certification Disclaimer (Graduate School Form 32), this thesis/dissertation adheres to the provisions of Purdue University’s “Policy of Integrity in Research” and the use of copyright material. Approved by Major Professor(s): Frederi Viens Approved by: Jun Xie 11/24/2015 Head of the Departmental Graduate Program Date MALLIAVIN CALCULUS IN THE CANONICAL LEVY´ PROCESS: WHITE NOISE THEORY AND FINANCIAL APPLICATIONS A Dissertation Submitted to the Faculty of Purdue University by Rolando D. Navarro, Jr.
    [Show full text]
  • Weak Approximations for Arithmetic Means of Geometric Brownian Motions and Applications to Basket Options Romain Bompis
    Weak approximations for arithmetic means of geometric Brownian motions and applications to Basket options Romain Bompis To cite this version: Romain Bompis. Weak approximations for arithmetic means of geometric Brownian motions and applications to Basket options. 2017. hal-01502886 HAL Id: hal-01502886 https://hal.archives-ouvertes.fr/hal-01502886 Preprint submitted on 6 Apr 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Weak approximations for arithmetic means of geometric Brownian motions and applications to Basket options ∗ Romain Bompis y Abstract. In this work we derive new analytical weak approximations for arithmetic means of geometric Brownian mo- tions using a scalar log-normal Proxy with an averaged volatility. The key features of the approach are to keep the martingale property for the approximations and to provide new integration by parts formulas for geometric Brownian motions. Besides, we also provide tight error bounds using Malliavin calculus, estimates depending on a suitable dispersion measure for the volatilities and on the maturity. As applications we give new price and im- plied volatility approximation formulas for basket call options. The numerical tests reveal the excellent accuracy of our results and comparison with the other known formulas of the literature show a valuable improvement.
    [Show full text]
  • Stochastic Differential Equations and Hypoelliptic
    . STOCHASTIC DIFFERENTIAL EQUATIONS AND HYPOELLIPTIC OPERATORS Denis R. Bell Department of Mathematics University of North Florida Jacksonville, FL 32224 U. S. A. In Real and Stochastic Analysis, 9-42, Trends Math. ed. M. M. Rao, Birkhauser, Boston. MA, 2004. 1 1. Introduction The first half of the twentieth century saw some remarkable developments in analytic probability theory. Wiener constructed a rigorous mathematical model of Brownian motion. Kolmogorov discovered that the transition probabilities of a diffusion process define a fundamental solution to the associated heat equation. Itˆo developed a stochas- tic calculus which made it possible to represent a diffusion with a given (infinitesimal) generator as the solution of a stochastic differential equation. These developments cre- ated a link between the fields of partial differential equations and stochastic analysis whereby results in the former area could be used to prove results in the latter. d More specifically, let X0,...,Xn denote a collection of smooth vector fields on R , regarded also as first-order differential operators, and define the second-order differen- tial operator n 1 L ≡ X2 + X . (1.1) 2 i 0 i=1 Consider also the Stratonovich stochastic differential equation (sde) n dξt = Xi(ξt) ◦ dwi + X0(ξt)dt (1.2) i=1 where w =(w1,...,wn)isann−dimensional standard Wiener process. Then the solution ξ to (1.2) is a time-homogeneous Markov process with generator L, whose transition probabilities p(t, x, dy) satisfy the following PDE (known as the Kolmogorov forward equation) in the weak sense ∂p = L∗p. ∂t y A differential operator G is said to be hypoelliptic if, whenever Gu is smooth for a distribution u defined on some open subset of the domain of G, then u is smooth.
    [Show full text]
  • An Intuitive Introduction to Fractional and Rough Volatilities
    mathematics Review An Intuitive Introduction to Fractional and Rough Volatilities Elisa Alòs 1,* and Jorge A. León 2 1 Department d’Economia i Empresa, Universitat Pompeu Fabra, and Barcelona GSE c/Ramon Trias Fargas, 25-27, 08005 Barcelona, Spain 2 Departamento de Control Automático, CINVESTAV-IPN, Apartado Postal 14-740, 07000 México City, Mexico; [email protected] * Correspondence: [email protected] Abstract: Here, we review some results of fractional volatility models, where the volatility is driven by fractional Brownian motion (fBm). In these models, the future average volatility is not a process adapted to the underlying filtration, and fBm is not a semimartingale in general. So, we cannot use the classical Itô’s calculus to explain how the memory properties of fBm allow us to describe some empirical findings of the implied volatility surface through Hull and White type formulas. Thus, Malliavin calculus provides a natural approach to deal with the implied volatility without assuming any particular structure of the volatility. The aim of this paper is to provides the basic tools of Malliavin calculus for the study of fractional volatility models. That is, we explain how the long and short memory of fBm improves the description of the implied volatility. In particular, we consider in detail a model that combines the long and short memory properties of fBm as an example of the approach introduced in this paper. The theoretical results are tested with numerical experiments. Keywords: derivative operator in the Malliavin calculus sense; fractional Brownian motion; future average volatility; Hull and White formula; Itô’s formula; Skorohod integral; stochastic volatility models; implied volatility; skews and smiles; rough volatility Citation: Alòs, E.; León, J.A.
    [Show full text]
  • A Minicourse on Stochastic Partial Differential Equations
    Robert C. Dalang Davar Khoshnevisan Carl Mueller David Nualart Yimin Xiao A Minicourse on Stochastic Partial Differential Equations Salt Lake City, Utah, 2006 Editors: D. Khoshnevisan and F. Rassoul-Agha Springer Berlin Heidelberg NewYork Hong Kong London Milan Paris Tokyo Contents A Primer on Stochastic Partial Differential Equations Davar Khoshnevisan:::::::::::::::::::::::::::::::::::::::::::::: 5 1 What is an SPDE? . 5 2 Gaussian random vectors . 6 3 Gaussian processes . 6 4 Regularity of random processes. 12 5 Martingale measures . 18 6 A nonlinear heat equation . 27 7 From chaos to order . 36 References . 40 The Stochastic Wave Equation Robert C. Dalang :::::::::::::::::::::::::::::::::::::::::::::::: 43 1 Introduction . 43 2 The stochastic wave equation . 45 3 Spatially homogeneous Gaussian noise . 51 4 The wave equation in spatial dimension 2. 53 5 A function-valued stochastic integral . 60 6 The wave equation in spatial dimension d ≥ 1.................. 62 7 Spatial regularity of the stochastic integral (d = 3). 65 8 H¨older-continuity in the 3-d wave equation . 74 References . 75 Application of Malliavin Calculus to Stochastic Partial Differential Equations David Nualart ::::::::::::::::::::::::::::::::::::::::::::::::::: 77 1 Introduction . 77 2 Malliavin calculus . 78 3 Application of Malliavin calculus to regularity of probability laws . 88 4 Stochastic heat equation . 96 5 Spatially homogeneous SPDEs . 104 VI Contents References . 113 Some Tools and Results for Parabolic Stochastic Partial Differential Equations Carl Mueller :::::::::::::::::::::::::::::::::::::::::::::::::::: 115 1 Introduction . 115 2 Basic framework . 117 3 Duality . 119 4 Large deviations for SPDEs . 129 5 A comparison theorem . 134 6 Applications . 136 References . 147 Sample Path Properties of Anisotropic Gaussian Random Fields Yimin Xiao ::::::::::::::::::::::::::::::::::::::::::::::::::::: 151 1 Introduction . 151 2 Examples and general assumptions .
    [Show full text]
  • Stochastic Processes and Applications Mongolia 2015
    NATIONAL UNIVERSITY OF MONGOLIA Stochastic Processes and Applications Mongolia 2015 27th July - 7th August 2015 National University of Mongolia Ulaanbaatar Mongolia National University of Mongolia 1 1 Final Report Executive summary • The research school was attended by 120 individuals, the majority of which were students, 54 of which were Mongolians. There was representation from academic institutions in 20 different countries. The lecture room remained full on every one of the 10 working days without numbers waning. • The whole event was generously sponsored by: Centre International de Math´ematiquesPures et Ap- pliqu´ees(CIMPA), Deutscher Akademischer Austauschdienst (DAAD), National University of Mongo- lia (NUM), Mongolian University of Science and Technology (MUST), the State Bank (SB), Mongolian Agricultural Commodities Exchange (MACE), Index Based Livestock Insurance Project (IBLIP) and Tenger Insurance, Mongolia (TIM). There was also a kind cash contribution from Prof. Ga¨etanBorot. The total expenditure of the school came to around 56K EURO. • Feedback indicates the event was an overwhelming scientific success and the school is likely to remain as a landmark event in the history of the Department of Mathematics at the National University of Mongolia. • There was outstanding cultural exchange, with strong social mixing and interaction taking place fol- lowing traditional norms of Mongolian hospitality. Participants experienced cultural activities on two Wednesday afternoons as well as a weekend excursion into Terelj National Park. These include throat singing, ger building, hiking, visiting a Buddhist Temple, museum exhibitions, learning about Mongo- lia's Soviet past, socialising at lunches and dinners with the help of an engineered `buddy system'. • There was significant exchange between senior Mongolian and foreign academics at the meeting re- garding the current and future plight of Mongolian mathematics.
    [Show full text]
  • Neural Stochastic Differential Equations: Deep Latent Gaussian
    Neural Stochastic Differential Equations: Deep Latent Gaussian Models in the Diffusion Limit Belinda Tzen∗ Maxim Raginskyy Abstract In deep latent Gaussian models, the latent variable is generated by a time-inhomogeneous Markov chain, where at each time step we pass the current state through a parametric nonlinear map, such as a feedforward neural net, and add a small independent Gaussian perturbation. This work considers the diffusion limit of such models, where the number of layers tends to infinity, while the step size and the noise variance tend to zero. The limiting latent object is an Itô diffusion process that solves a stochastic differential equation (SDE) whose drift and diffusion coefficient are implemented by neural nets. We develop a variational inference framework for these neural SDEs via stochastic automatic differentiation in Wiener space, where the variational approximations to the posterior are obtained by Girsanov (mean-shift) transformation of the standard Wiener process and the computation of gradients is based on the theory of stochastic flows. This permits the use of black-box SDE solvers and automatic differentiation for end-to-end inference. Experimental results with synthetic data are provided. 1 Introduction Ordinary differential equations (ODEs) and other types of continuous-time flows have always served as convenient abstractions for various deterministic iterative models and algorithms. Recently, however, several authors have started exploring the intriguing possibility of using ODEs for constructing and training very deep neural nets by considering the limiting case of composing a large number of infinitesimal nonlinear transformations (Haber and Ruthotto, 2017; Chen et al., 2018b; Li et al., 2018).
    [Show full text]
  • Martingale Structure of Skorohod Integral Processes
    Martingale Structure of Skorohod Integral Processes Giovanni PECCATI Laboratoire de Statistique Th´eorique et Appliqu´ee Universit´eParis VI, France email: [email protected] Mich`eleTHIEULLEN Laboratoire de Probabilit´eset Mod`elesAl´eatoires Universit´esParis VI and Paris VII, France email: [email protected] Ciprian A. TUDOR∗ Laboratoire de Probabilit´eset Mod`elesAl´eatoires Universit´esParis VI and Paris VII, France email: [email protected] February 3, 2005 Abstract Let the process {Yt, t ∈ [0, 1]}, have the form Yt = δ u1[0,t] , where δ stands for a Skorohod integral with respect to Brownian motion, and u is a measurable process verifying some suitable regularity conditions. We use a recent result by Tudor (2004), to prove that Yt can be represented as the limit of linear combinations of processes that are products of forward and backward Brownian martingales. Such a result is a further step towards the connection between the theory of continuous- time (semi)martingales, and that of anticipating stochastic integration. We establish an explicit link between our results and the classic characterization, due to Duc and Nualart (1990), of the chaotic decomposition of Skorohod integral processes. We also explore the case of Skorohod integral processes that are time-reversed Brownian martingales, and provide an “anticipating” counterpart to the classic Optional Sampling Theorem for Itˆostochastic integrals. Key words – Malliavin calculus; Anticipating stochastic integration; Martingale theory; Stopping times. AMS 2000 classification – 60G15; 60G40; 60G44; 60H05; 60H07 Running title –Martingale structure of integrals 1 Introduction Let C[0,1], C, P = (Ω, F, P) be the canonical space, where P is the law of a standard Brownian motion started from zero, and write X = {Xt : t ∈ [0, 1]} for the coordinate process.
    [Show full text]
  • Itô Calculus and Malliavin Calculus
    ItˆoCalculus and Malliavin Calculus Shinzo Watanabe 527-10, Chaya-cho, Higashiyama-ku, Kyoto, 605-0931, Japan. Dedicated to Professor Kiyosi Itˆoon his 90th birthday 1 Introduction. Since the Wiener space was established by N. Wiener as a mathematical model of Brownian motion in 1923, a rigorous theory of integrations on a function space started. In these almost eighty years, it has been providing us with important methods in stochastic analysis and its applications. Around 1942, R. Feynman ([F 1],[F 2]) had an epoch making idea of rep- resenting the propagators for Schr¨odinger equations by a path integral over trajectories of quantum mechanical particles. M. Kac noticed that its coun- terpart could be discussed rigorously on Wiener space and thus found the Feynman-Kac formula. He also applied probabilistic representations of heat kernels by Wiener functional expectations to study asymptotics of spectra of Schr¨odinger operators ([Ka 1],[Ka 2],[Ka 3]). This study was further developed in a fundamental paper by McKean and Singer [MS], which may be regarded as an origin of the heat equation methods in the analysis of manifold. The approach by McKean and Singer is based on PDE theory, the method of parametrix for heat kernels, in particular. If we would give a similar prob- abilitistic approach as Kac in the problems of McKean and Singer, we have several difficulties to overcome. In the case of Kac, the second order term of the Schr¨odinger operator is Laplacian in Euclidean space so that a use of Wiener process and pinned Wiener process is sufficient, which could be easily set up on a Wiener space.
    [Show full text]
  • Malliavin Calculus and Its Applications
    B.Sc. Thesis Malliavin calculus and its applications Ad´amGyenge´ Supervisor: Dr. Tam´asSzabados Associate professor Institute of Mathematics Department of Stochastics Budapest University of Technology and Economics 2010 Contents 1 Introduction 1 1.1 The structure of the thesis . 1 2 Theory of Malliavin calculus 3 2.1 It^ointegral . 3 2.1.1 Brownian motion . 3 2.1.2 Construction of the It^ointegral . 4 2.2 Isonormal Gaussian processes . 5 2.3 Wiener chaos expansion . 6 2.3.1 Iterated It^o-and multiple stochastic integrals . 8 2.3.2 Wiener chaos and iterated It^ointegrals . 9 2.4 The derivative operator . 11 2.4.1 Closability and other properties of the derivative operator . 12 2.4.2 The derivative in the white noise case . 14 2.5 The divergence operator and the Skorohod integral . 16 2.5.1 The Skorohod integral . 17 2.6 The Clark-Ocone formula . 20 3 Applications 22 3.1 Stochastic differential equations . 22 3.1.1 Absolute continuity of distributions . 22 3.1.2 Lie bracket and H¨ormander'scondition . 24 3.1.3 Absolute continuity under H¨ormander'scondition . 25 3.2 Relation to Stein's method . 28 3.2.1 Stein's lemma . 28 3.2.2 Ornstein-Uhlenbeck operators . 29 3.2.3 Bounds on the Kolmogorov distance . 30 3.3 Financial mathematics . 31 3.3.1 Investments . 31 3.3.2 Black-Scholes model . 34 4 Summary and conclusions 36 i Chapter 1 Introduction The Malliavin calculus (or the stochastic calculus of variations) is an infinite dimensional differential calculus on the Wiener space.
    [Show full text]
  • The Malliavin Calculus
    THE MALLIAVIN CALCULUS Han Zhang Supervisor: Professor Tony Dooley School of Mathematics, The University of New South Wales. December 2004 Submitted in partial fulfillment of the requirements of the degree of Bachelor of Science with Honours Acknowledgements I enjoyed the help and support from numerous friends and faculty members during whilst writing this thesis. My greatest debt however, goes to my supervisor, Pro- fessor Tony Dooley, who initiated me to the fantastic field of Malliavin calculus. With his incredible breadth and depth of knowledge and intuition, he guided me to structure and clarify my thought and suggested valuable insightful comments. This thesis would not have been possible to write without his help and support. I am also indebted to express my sincere thankfulness to the honours coordi- nators, Dr Ian Doust and Dr Brian Jefferies, who were very kind and supportive to us during the year. Dr Ben Goldys also deserves a special mention, for the rec- ommendation of a number of useful references, and I very much enjoyed the casual discussions I had with him. Finally, I would like to thank everyone in the School of Mathematics, and the Department of Actuarial Studies for providing me with such a wonderful four years at UNSW. Han Zhang, November 2004. i History and Introduction The Malliavin calculus, also known as the stochastic calculus of variations, is an infinite dimensional differential calculus on the Wiener space. Much of the theory builds on from Itˆo’s stochastic calculus, and aims to investigate the structure and also regularity laws of spaces of Wiener functionals.
    [Show full text]