The Donsker delta function and local time for McKean-Vlasov processes and applications

The purpose of this paper is to establish a stochastic diﬀerential equation for the Donsker delta measure of the solution of a McKean-Vlasov (mean-ﬁeld) stochastic diﬀerential equation. If the Donsker delta measure is absolutely continuous with respect to Lebesgue measure, then its Radon-Nikodym derivative is called the Donsker delta function. In that case it can be proved that the local time of such a process is simply the integral with respect to time of the Donsker delta function. Therefore we also get an equation for the local time of such a process. For some particular McKean-Vlasov processes, we ﬁnd explicit expressions for their Donsker delta functions and hence for their local times.


Introduction
The Donsker delta function of a random variable or a stochastic process arises in many studies, including quantum mechanical particles on a circle [7], financial markets with insider trading as in [11], and in [3] for financial markets with singular drift.It has also been used as a tool to determine explicit formulae for replicating portfolios in complete and incomplete markets, see [10].
Moreover, the Donsker delta function is also of interest because it can be regarded as a time-derivative of the local time.Therefore, explicit expressions for the Donsker delta function lead to explicit formulae of the local time.
For example, if we let B be a Brownian motion defined on a filtered probability space (Ω, F , F = {F t } t≥0 , P ), then the Donsker delta function δ B(t) (x) of a Brownian motion B at the point x can be regarded as the time derivative of the local time L t (x) of B. More precisely, we have Such an integral exists as an element of the Hida space (S * ) of stochastic distributions.See Section 2. 3.
In [6] the authors use white noise theory to obtain an explicit solution formula for a general stochastic differential equation (SDE), and this is used to find an expression for the Donsker delta function for the solution of an SDE.Subsequently this was also extended to SDEs driven by Lévy noise in [9].
The main result of the current paper is that the Donsker delta measure of a McKean-Vlasov process (see below) always satisfies a certain Fokker-Planck type SPDE in the sense of distributions.Moreover, we use this to find explicit formulae for the Donsker delta functions for McKean-Vlasov processes, and hence their local times, in specific cases.
We call X a McKean-Vlason process.
Here the σ-algebra F = {F t } t≥0 denotes the filtration generated by Z and B(•), Z is a random variable which is independent of the σ-algebra generated by B(•) and such that We denote by µ t = µ t (ω) = µ t (dx, ω) the conditional law of X(t) given the filtration F t generated by the Brownian motion B.More precisely, we consider the following model: Since we consider only a one-dimensional Brownian motion B(t) ∈ R, we will show that the regular conditional distribution of X(t) given the filtration {F t } t≥0 can be identified with the Donsker delta measure in the sense of distribution.See details in Section 3.1

Preliminaries
In this section we review some basic notions and results that will be used throughout this work.

Radon measures
A Radon measure on R d is a Borel measure which is finite on compact sets, outer regular on all Borel sets and inner regular on all open sets.In particular, all Borel probability measures on R d are Radon measures.In the following, we let • M 0 be the set of deterministic Radon measures.
• C 0 (R d ) be the uniform closure of the space C c (R d ) of continuous functions with compact support.
If we equip M 0 with the total variation norm ||µ|| := |µ|(R d ), then M 0 becomes a Banach space, and it is the dual of C 0 (R d ).See Chapter 7 in Folland [5] for more information.
If µ ∈ M 0 is a finite measure, we define to be the Fourier transform of µ at y.In particular, if µ(dx) is absolutely continuous with respect to Lebesgue measure dx with Radon-Nikodym-derivative m(x) = µ(dx) dx , so that µ(dx) = m(x)dx with m ∈ L 1 (R d ), we define the Fourier transform of m at y, denoted by m(y) or F [m](y), by We let M denote the set of all random measures µ(dx, ω); ω ∈ Ω such that µ(dx, ω) ∈ M 0 for each given ω ∈ Ω.

The Schwartz space of tempered distributions
We recall now some notions from white noise analysis.
• S = S(R d ) be the Schwartz space of rapidly decreasing smooth real functions on R d .

The Hida space (S) * of stochastic distributions
We restrict ourselves to the white noise probability space (Ω = S ′ , F = B, P ), where B is the Borel σ-algebra and the probability P is the probability measure on S ′ defined in virtue of the Bochner-Minlos-Sazonov theorem).Let J denote the set of all finite multi-indices α = (α 1 , α 2 , . . ., α m ), m = 1, 2, . .., of non-negative integers α i .
• ((S) k ) k∈R is the Hilbert space consisting of all f = α∈J c α H α ∈ L 2 (P ) such that • The space (S) = k∈R (S) k equipped with the projective topology is the Hida space of stochastic test functions.
• ((S) −k ) k∈R is the Hilbert space consisting of all formal sums F = α∈J c α H α equipped with the norm • The space (S) * = k∈R (S) −k equipped with the inductive topology is the Hida space of stochastic distributions.It can be regarded as the dual of (S).

The Donsker delta function
We now recall some basic definitions: Definition 2.1 Let Y : Ω → R be a random variable which also belongs to the Hida space (S) * of stochastic distributions.Then a continuous function for all (measurable) g : R → R such that the integral converges in (S) * .
The Donsker delta function is related to the regular conditional distribution.The connection is the following: The regular conditional distribution with respect to the σ-algebra F of a given real random variable Y , denoted by µ Y (dy) = µ Y (dy, ω); ω ∈ Ω, is defined by the following properties: • For each fixed ω ∈ Ω, µ Y (dy, ω) is a probability measure on the Borel subsets of R.
It is well-known that such a regular conditional distribution always exists.See e.g.[4], p.79.From the required properties of µ Y (dy, ω), we get the following formula: (2.7) We will prove in Theorem 3.3 that the Donsker delta function can be regarded as a stochastic distribution in S ′ , satisfying a Fokker-Planck type SPDE in the sense of distributions.It can also be represented as an element of the Hida stochastic distribution space (S) * , and as such it can in some cases be expressed explicitly in terms of Wick calculus.For example, if Y (t) = B(t), we have where ⋄ denotes Wick multiplication and exp ⋄ denotes Wick exponential.Note that even though the Donsker delta function can only be represented as a distribution, its conditional expectation can be a real valued stochastic process.For example, for t < T we have (2.9) For more examples, we refer to e.g.[1] or [10].
3 The Donsker delta equation for McKean-Vlasov processes

The general multidimensional Fokker-Planck equation
To explain the background for this section, let us recall the general multidimensional situation studied in [2], where X(t) ∈ R d is a McKean-Vlasov diffusion, of the form (using matrix notation), where B is a multi-dimensional Brownian motion.
Here Z is a random variable which is independent of the σ-algebra generated by B(•) and such that Define the σ-algebra F = {F t } t≥0 to be the filtration generated by Z and B(•).
Let M denote the set of all Borel measures on R d .We assume that the coefficients b(t, x, µ) bounded and F-predictable processes for all x, µ, and that b and σ are continuous with respect to t and x for all µ.
One can check that under some assumptions, such as Lipschitz and linear growth conditions, there exists a unique solution of equation (3.1).
We denote by µ t = µ t (ω) = µ t (dx, ω) the conditional law of X(t) given the filtration F t generated by the Brownian motion B.More precisely, we consider the following model: Definition 3.1 Fix one of the Brownian motions, say B 1 = B 1 (t, ω), with filtration {F (1) t } t≥0 .We define µ t = µ X(t) (ω, dx) to be regular conditional distribution of X(t) given F (1) t .This means that µ t (ω, dx) is a Borel probability measure on R d for all t ∈ [0, T ], ω ∈ Ω and for all functions g such that The following version of the stochastic Fokker-Planck integro-differential equation for the conditional law for McKean-Vlasov jump diffusions was proved by Agram and Øksendal [2].For simplicity we consider only the case without jumps here.Theorem 3.2 (Conditional stochastic Fokker-Planck equation [2]) Let X(t) be as in (3.1) with m ≥ 2 and let µ t := µ X(t) (dx, ω) be the regular conditional distribution of X(t) given F (1) t .Then for a.a.ω ∈ Ω the conditional law µ t ∈ S ′ and it satisfies the following SPDE (in the sense of distributions): Here A * 0 , A * 1 are the integro-differential operator and the differential operator which are given respectively by: and In the above D j , D n,j denote ∂ ∂x j and ∂ 2 ∂xn∂x j respectively, in the sense of distributions.

The Fokker-Planck equation for the Donsker measure
In [2] the theorem above was proved under the assumption that m ≥ 2. However, the proof also works if m = 1 and F (1) t = F t .Note that in this case, since X(t) is F t -measurable, the identity (3.2) states that for all functions g such that R d |g(x)|µ t (dx, ω) < ∞.
In particular, if we choose d = m = 1 in the above we get that the conditional law coincides with the Donsker measure, i.e. (3.7) Therefore we get the following Fokker-Planck equation for the Donsker measure: Then the Donsker delta measure µ t = δ X(t) (dx, ω) satisfies the following equation (in the sense of distribution) where D = ∂ ∂x and

Local time
In this section we first recall the definition of local time of a stochastic process Y (•): where λ denotes Lebesgue measure on R and the limit is in L 2 (P ).
In the white noise context the local time can be represented as the integral of the Donsker delta function.More precisely, we have the following result: Theorem 4.2 The local time L t (x) of X at the point x and the time t is given by where the integration takes place in (S) * (or in S ′ for each ω).
Proof.For completeness we give the proof.By definition of the local time and the Donsker delta function, we have because the function y → δ Y (s) (y) is continuous in (S) * (and in S ′ ).
Remark 4.3 Note that even though we in general can only say that δ X(t) (x) ∈ (S) * , L t (x) usually exists as a real-valued stochastic process.

Explicit solutions
In this Section, we find explicitly the Donsker delta function for some particular McKean-Vlasov processes and accordingly their local time.Suppose that µ t is absolutely continuous i.e.

Brownian motion
Consider the special case when X(t) = B(t); B(0) = Z.Then b = 0 and σ = 1 and equation (3.8) becomes We can easily verify by Wick calculus that a solution in (S) * of equation (5.5) is which is in agreement with (2.8).The details are as follows: Try Then and and Collecting the terms we see that satisfies the Fokker-Planck equation (5.5) for the conditional law of B(t).
From white noise theory we know that for all random variables X, Y with a finite expectation (independent or not).From this we see that which has a singularity at x = z.

Coefficients not depending on x
The next result shows that, under some conditions, the Donsker delta function can be an ordinary function if the initial value X(0) has a density: Theorem 5.2 Assume that X(t) is the solution of the following McKean-Vlasov equation: where the coefficients α(t, x, µ) = α(t, µ) and β(t, x, µ) = β(t, µ) do not depend on x.Suppose that X(0) = Z is a random variable (independent of B) with density (5.11) (5.12) (5.13) Then Y (t, x) is the Donsker delta function of X(t).
2. The solution X(t) of (5.10) is given by 1. We show that Y (t, x) satisfies equation (5.3).By the Ito formula we have we see that the equation (5.15) can be written which is the same as equation (5.3).Since Y (0, x) = h(x) = m(0, x) we conclude by uniqueness that Y (t, x) = m(t, x) for all t.
2. This follows from the definition of the Donsker delta function.

Constant coefficients
As a special case of the case above, suppose that where α and β are constants.Then by Theorem 5.2 the Donsker delta function is (5.19)

Mean-field geometric Brownian motion
Suppose that X(t) is a McKean-Vlasov process of the form We call this a mean-field geometric Brownian motion.For such processes we have: The Donsker delta function m t (x) for the mean-field geometric Brownian motion X(t) is where (ii) The solution X(t) of the mean-field geometric Brownian motion equation (5.20) can be written (i) The corresponding Fokker-Planck equation for the Donsker delta function This is a stochastic partial differential equation in m t (x).It seems difficult to find directly an explicit solution of this equation.However, we can find the solution m t (x) = δ X(t) (x) by proceeding as follows: The solution of (5.20) is where By Theorem 5.2 we know that where By definition we have With g(y) = exp(y) this gives R exp(y)δ Yt (y)dy = exp(Y t ) = X t .
Hence, substituting exp(y) = x, From this we deduce that is the Donsker delta function of X t .
(ii) This part follows by the definition of the Donsker delta function.

An example related to the Burgers equation
then we see that the Burgers equation (5.27) becomes the following equation in ψ: (5.30) Integrating with respect to x this gives for some constant γ.Then in terms of ϕ the above equation gets the form .33)This simplifies to the equation for ϕ reduces to the (linear) stochastic heat equation ϕ(0, x) = k(x) (to be determined), (5.37) or, using Ito differential notation, (5.39) To find an expression for the solution of (5.36), define an auxiliary process R(t) = R (x) (t) by where B is an auxiliary Brownian motion with law P and independent of B. Then by the Feynman-Kac formula where Ẽ denotes expectation with respect to P and k(z) = ϕ(0, z), solves equation (5.36).
Going back to m we get In particular, setting t = 0 we get We summarize what we have proved as follows: Theorem 5.4 1.The Donsker delta function m(t, x) = δ X(t) (x) for the solution X(t) of the McKean-Vlasov equation (5.26) is given by where  For comparison, recall that the density of Brownian motion at t, x (when starting at x 0 ) is ).

Acknowledgments
We are grateful to Frank Proske for helpful comments.

Definition 4 . 1
The local time L t (y) of Y (•) at the point y and at time t is defined by .31) Now define the function ϕ = ϕ(t, x) by ψ = γ ln ϕ, (5.32)