Translation Invariant Diffusions and Stochastic Partial
Differential Equations in S′
B.Rajeev
Email: [email protected]
Abstract
In this article we show that the ordinary stochastic differential
equations of K.Itô maybe considered as part of a larger class of
second order stochastic PDE’s that are quasi linear and have the
property of translation invariance. We show using the ‘monotonicity
inequality’ and the Lipshitz continuity of the coefficients
σij and bi, existence and uniqueness of strong
solutions for these stochastic PDE’s. Using pathwise uniqueness, we
prove the strong Markov property.
Keywords :S′ valued process, diffusion
processes, Hermite-Sobolev space, Strong solution, quasi linear
SPDE, Monotonicity inequality, Translation invariance
Subject
classification :[2010]60G51, 60H10, 60H15
1 Introduction
The notion of an ordinary stochastic differential equation (SDE) was
introduced by K.Itô in [23] and since then has become the
main tool for modelling diffusion phenomena as a random process (see
for example [34]). The approach to diffusions as a random
process goes back to the works of A.N.Kolmogorov [26], and was
studied by N.Weiner [48], W.Feller[16], J.L.Doob [13]
and P.Lévy [31]. The theory was extended further by
D.W.Stroock and S.R.S.Varadhan in their well known ‘weak
formulation’ or ‘martingale formulation’ [43]. The subject of
stochastic partial differential equations (SPDE) on the other hand,
is of more recent vintage ([47],[30] ). It extends the
logic of perturbing an ordinary differential equation by noise,
inherent in the Itô approach, to partial differential equations.
Although the underlying probabilistic logic is the same, the
mathematics of these two models can be vastly different, the latter
more often than not involving the tools and techniques of function
space analysis (see for example [21], [12]). On the other
hand, one of the fundamental features that continues to sustain
interest in the Itô approach, both in applications and theory,
is the connection with other areas of mathematics like partial
differential equations and potential theory
([29],[3],[42]); more recent examples are the notion
of ‘viscosity solutions’ related to the Hamilton-Jacobi-Bellman
equation ([32]), and backward SDE’s ([35]). In this paper,
we show that the two approaches viz. the SDE and the SPDE approaches
can be unified into a single framework, in which the SDE approach
(with an extra parameter) is equivalent to the SPDE approach and
mathematically speaking both maybe viewed as part of a single
structure. Our method may be considered a variant of the well known
‘ method of characteristics’ in PDE, that constructs solutions of
PDE’s from the ordinary differential equations satisfied by the
‘characteristic curves’ associated with the PDE (see [15],
Chapter 3, and [30], Chapter 6, for the stochastic case). The
difference in our approach lies in the treatment of non linearities
i.e. in the manner in which the coefficients in the PDE or SPDE are
allowed to depend on the solutions. In this paper, we first
construct the solutions of the SPDE and then deduce the solutions of
the corresponding SDE. The reverse construction of solutions of
SPDE’s from that of the associated SDE’s was already done, via the
Itô formula, in [36],[39]. It turns out that the
solutions of Itô’s SDE’s correspond to rather singular solutions
of the associated (quasi linear) SPDE in a manner analogous to the
way in which ‘fundamental solutions’ are associated to certain
second order partial differential equations. The solutions of the
SPDE so constructed, arise, in a unique fashion, as translations of
the initial condition of the SPDE by the solution of the
‘characteristic’ SDE starting at the origin.
In more detail, we construct in this paper a general method of
solving the stochastic partial differential equation (SPDE) driven
by an n-dimensional Brownian motion (Bt) in the form
[TABLE]
Here L and
A=(A1,⋯,An) are non linear partial differential operators
of the second and first order respectively on the space of tempered
distributions S′→S′ on Rd
given by equations (2) and (3) below. The initial condition y
is an arbitrary tempered distribution whose regularity maybe
measured on a decreasing scale of Hilbert spaces Sp,p∈R. In particular y∈Sp for some p∈R. The operators L and Ai are quasi-linear i.e. they are
constant coefficient differential operators once the value of the
coefficients σij,bi:Sp→R are
fixed, L is of order two and the Ai’s of order one.
Consequently, L,Ai:Sp→Sq,q≤p−1,i=1,⋯,n (see Section 2). Thus the domain and range of
the operators L and Ai differ, leading to what K.Itô in
[24] refers to as a ‘Type 2’ equation. This also introduces
the principal difficulty in solving our SPDE, since there is no
obvious way of using techniques such as ‘Picard iteration’. However
by assuming a Lipschitz condition on σij and bi with
respect to the norm ∥⋅∥q,q≤p−1 (equation (4)) and
exploiting the quasi linear structure of the operators L,Ai we
implement a modified form of Picard iteration to solve the above
SPDE. The solutions of the above equation have the property that
they are translation invariant i.e. they can be written as Yt(y)=τZt(y)y, where τx:S′→S′
are the translation operators and (Zt(y)) is a finite dimensional
process that depends on the initial value y. This has the
consequence that the solution corresponding to the translate
τxy is the translate of y by the process (x+Zt(τxy)). Note that the Sp themselves are invariant
under translations i.e. τx:Sp→Sp
([37]). The action of the translation operators on y gives
rise to finite dimensional coefficients σˉij,bˉi,i=1,⋯,d,j=1,⋯n by σˉij(z):=σij(τzy),bˉi(z):=bi(τzy),z∈Rd. It turns out that Xtx:=x+Zt(τxy) solves the
ordinary stochastic differential equation driven by (Bt) with
coefficients σˉij,bˉi and initial value X0x=x. In recent times distribution dependent SDE’s have become an
active area of research (see for example
[2],[11],[27],[1],[10] and
references therein). We refer to Example 6 in Section 6 below, for
some connections between distribution dependent SDE’s and our
results. Our work also relates to the problem of identifying ‘
invariant submanifolds’ of solutions of SPDEs that arise in finance
(see [8],[9],[14],[44]). In effect, the
set of translates {τxy:x∈Rd} serves as an
invariant manifold for the above SPDE with initial distribution y∈Sp, under some smoothness assumptions on y.
Our method relies on three ingredients viz. one, a quasi-linear
extension of linear differential operators by identifying the
coefficients σij(x) as a restriction of the functional
⟨σij,ϕ⟩,ϕ∈S−p=Sp′,p>d,σij∈Sp to the distribution ϕ=δx ;
two, an Itô formula for translations of tempered distributions
by semi-martingales (see [36], [4],[46]); and
finally, the monotonicity inequality (see [5],[20]).
Indeed, this last inequality, whose abstract version has been known
for some time (see [28], [25], [18]), has proved to
be an indispensable tool for proving uniqueness results for SPDE’s
in the framework of a scale of Hilbert spaces of the type discussed
above(see [20],[39]). Our results below show that it
can also be used for proving existence results.
The paper is organised as follows. After the preliminaries in
Section 2, we prove in Section 3, using the monotonicity inequality
(Theorem 3.1), some extensions of the same in Theorems (3.2) and
(3.3) respectively; viz. in the case that the pair of operators
(A,L) have variable coefficients. These inequalities are crucial
for the convergence results in Section 4 which contains the main
existence and uniqueness results in Theorem (4.3). Our proof of
existence is tailored for the infinite dimensional situation and
applies to more general situations. A simpler proof is indicated in
Remark 4.4. In Section 5, we construct the ‘maximal’ solutions upto
an explosion time and prove the strong Markov property (Theorem
(5.6)) using the pathwise uniqueness established in Section 4 . In
section 6, we look at several examples. Examples 1,2 & 3 relate to
finite dimensional diffusions. Example 4 relates solutions of our
SPDE with solutions of the associated martingale problem for L.
Example 5 deals with the stochastic representation of the solutions
of non-linear evolution equation canonically associated with the
operator L. In Example 6, we consider the situation where the
coefficients in the finite dimensional equation depends on the
marginal law of the process. Finally Example 7 deals with extensions
of the operator L, that have a zeroth order term and is related to
the Feynman-Kac formula. In Section 7 we make some remarks on
‘duality’ and invariant measures in the context our SPDE. Some
technical results are in the Appendix. We use well known results on
stochastic calculus for processes with values in a Hilbert space,
for the proofs of which we refer to [12],[18],[33].
2 Preliminaries
Let (Ω,F,{Ft}t≥0,P) be a filtered probability space
satisfying the usual conditions viz. 1) (Ω,F,P) is a complete probability space. 2)
F0 contains all A∈F, such that P(A)=0, and 3) Ft=s>t⋂Fs,t≥0. On this probability space is given a standard
n-dimensional Ft- Brownian motion (Bt)≡(Bt1,…,Btn). We will denote the filtration generated by (Bt)
as (FtB). Let σˉij,bˉi be locally
Lipshitz functions on Rd for i=1,⋯,d,j=1,⋯n. Let σˉ:=(σˉij) (so that
(σˉij(x)),x∈Rd is a d×n
matrix) and bˉ:=(bˉ1,⋯,bˉd) be a vector
field on Rd. We use the notation R^d=Rd∪{∞} for the one point compactification
of Rd.
Theorem 2.1
*Let σˉ,bˉ,(Bt) be as above and x∈Rd. Then
∃ η:Ω→(0,∞],η an (FtB) stopping time and an R^d-valued,
(FtB) adapted process (Xt)t≥0 such that
-
*For all ω∈Ω, X.(ω):[0,η(ω))→Rd, is continuous
and *Xt(ω)=∞, t≥η(ω)
2. 2.
*a.s. (P), η(ω)<∞ implies *t↑η(ω)limXt(ω)=∞.
3. 3.
a.s.(P),
[TABLE]
for 0≤t<η(ω).
The solution (Xt,η) is (pathwise)
unique i.e. if (Xt1,η1) is another solution then P{Xt=Xt1,0≤t<η∧η1}=1.
Proof : We refer to [22], Chapter IV,Theorem 2.3 and
Theorem 3.1 for the proofs (with appropriate modifications for the
case d=r) of existence and uniqueness
respectively.\hfill□
Let α,β∈Z+d:={(x1,⋯,xd):xi≥0,xi integer}. Let xα be the product
xα:=x1α1…xdαd∈R and
∂β:=∂1β1…∂dβd, the differential operator of order β1+⋯βd
corresponding to the monomial xβ.
For a multi index α, we use the notation ∣α∣:=i=1∑dαi. Let S denote the space
of rapidly decreasing smooth real functions on Rd with
the topology given by the family of semi norms
{∧α,β}, defined for f∈S and
multi indices α,β by ∧α,β(f):=xsup ∣xα∂βf(x)∣. Then {S,∧α,β:α,β∈Z+d} is a locally convex,
complete, metrisable topological vector space i.e. a Fréchet
space. S′ will denote its continuous dual. The duality
between S and S′ will be denoted by ⟨ψ,ϕ⟩ for ϕ∈S and ψ∈S′.
For x∈Rd the translation operators τx:S→S are defined as τxf(y):=f(y−x) for
f∈S and then for ϕ∈S′ by duality : ⟨τxϕ,f⟩:=⟨ϕ,τ−xf⟩.
Let {hk;k∈Z+d} be the orthonormal basis in
the real Hilbert space L2(Rd,dx)⊃S
consisting of the Hermite functions (see for eg. [45]); here
dx denotes Lebesgue measure, where the dependence on the dimension
is suppressed whenever there is no risk of confusion. Let ⟨⋅,⋅⟩0 be the inner product in L2(Rd,dx). For f∈S and p∈R define the
inner product ⟨f,g⟩p on S as follows :
[TABLE]
The corresponding norm will be denoted by ∥⋅∥p. We define the Hilbert space Sp as the
completion of S with respect to the norm ∥⋅∥p
over the field of real numbers. The following basic relations hold
between the Sp spaces (see for eg. [24], [25]):
For 0<q<p, S⊂Sp⊂Sq⊂L2=S0⊂S−q⊂S−p⊂S′. Further,
S′=p∈R⋃Sp and
p∈R⋂Sp=S. If
{hkp:k∈Z+d} denotes the orthonormal basis in
Sp consisting of the (normalised) Hermite functions hkp:=(2∣k∣+d)−phk, then the dual space Sp′ may be
identified with S−p, via the basis {hk−p:k∈Z+d} of S−p. For ϕ∈S and
ψ∈S′ the bilinear form (ψ,ϕ)→⟨ψ,ϕ⟩ also gives the duality between Sp(⊃S) and S−p(⊂S′). It
is also well known that ∂i:Sp→Sp−21 are bounded linear operators for every p∈R and i=1,⋯,d. For Banach spaces Xand Y,L(X,Y) will denote the Banach space of bounded linear
operators from X into Y.
Let p∈R and let σij,bi:Sp→R,i=1,⋯,d,j=1,⋯,n. We consider
the (non-linear) operators A:=(A1,⋯,An):Sp→L(Rn,Sp−21),
from Sp to the space of linear operators from Rn to Sp−21, defined by
[TABLE]
and the non-linear operator L:Sp→Sp−1 defined as follows :
[TABLE]
where
aij(ϕ):=(σ(ϕ)σ(ϕ)t)ij and the
superscript ‘t’ denotes matrix transpose. Clearly if
σij(ϕ) and bi(ϕ) are bounded on the set {ϕ∈Sp:∥ϕ∥p≤λ} for some λ>0, then ∃ C=C(λ)>0 such that if q≤p−1 and {ei:i=1,⋯,n} is the standard
orthonormal basis in Rn, then
[TABLE]
for ϕ∈Sp with ∥ϕ∥p≤λ. In the above equalities and in
what follows we use the notation A(ϕ)⋅h:=i=1∑nAi(ϕ)⋅hi,h∈Rn. The
subscript ‘HS’ refers to the Hilbert-Schmidt norm. The above
inequalities follow from the boundedness of the operators
∂i:Sp→Sp−21
and the assumed (local) bounds on the coefficients σij and
bi.
3 The Monotonicity Inequality
In this section we will prove the ‘monotonicity inequality’
involving the pair (L,A) defined in equations (2) and (3) and
which we will use in the proof of existence and uniqueness of the
SPDE (18). The constant coefficient case was proved in [20].
Using techniques developed in [5] we prove the corresponding
inequality and a variant of the same in the variable coefficient
case, in theorems (3.2) and (3.3) below.
Let σ=(σij)∈Rdn, h=(h1…hn)∈Rn and ϕ∈Sp. Then we define
A0:Rdn×Sp→L(Rn,Sq), a bilinear map, as follows:
[TABLE]
Note that the symbols σij,σij(⋅) have different meanings, the latter being a
function on Sp and the former an element of R. A
similar remark holds for bi and bi(⋅). For ψ,ϕ∈Sp, and q≤p−1 we can write (from (2))
[TABLE]
Similarly
we can write
[TABLE]
where L1,L2 are Sq valued, bilinear maps on Rd×Sp and Rd2×Sp,
respectively, given as follows : Let (b,ϕ)∈Rd×Sp,b:=(b1,⋯,bd). Then
[TABLE]
and to define
L2, let (σ,ϕ)∈Rd2×Sp,σ:=(σij). Then,
[TABLE]
Note that for ϕ∈Sp, we have the d×d
matrix a(ϕ)≡σσt(ϕ):=((σ(ϕ)σt(ϕ))ij) and the d-dimensional vector b(ϕ):=(b1(ϕ),⋯,bd(ϕ)). For λ>0, define the
constant K1(λ) as follows :
[TABLE]
We then have the
following restatement of the Monotonicity inequality for constant
coefficient operators [20].
Theorem 3.1
Let p∈R,q≤p−1. Suppose that σij(⋅),bi(⋅),
i=1,…,d, j=1,…n are bounded on the set
Bp(0,λ):={ϕ∈Sp:∥ϕ∥p≤λ} for every λ>0. Then for every λ>0, ∃
a constant C=C(n,d,p,K1(λ))>0 such that
[TABLE]
for all ψ∈Sp and for all ϕ∈Bp(0,λ).
Proof: It follows from the Monotonicity inequality (see
[20],[5]) that the inequalities in the statement of
the theorem holds for fixed ψ,ϕ∈Sp with a
constant C′ that depends quadratically on the numbers max{∣σij(ϕ)∣:i=1,…,d,j=1,…,n} and linearly
on max{∣bi(ϕ)∣:i=1,…,d}. Taking supremum over
∥ϕ∥p≤λ, we get the required constants.
\hfill□
We now prove the Monotonicity
inequality in the form required to obtain uniqueness of solutions to
our stochastic partial differential equation (18) below.
Theorem 3.2
Let p∈R,q≤p−1. Let
σij,bi:Sp→R,
i=1,…,d, j=1,…n. Suppose that for λ>0, ∃ K(λ)>0 such that
[TABLE]
for ϕ,ψ∈ Bp(0,λ). Then ∃
a constant C=C(n,d,p,q,λ,K(λ),K1(λ)) such that
[TABLE]
for all ϕ,ψ∈Bp(0,λ).
Proof: Using the notation established in the discussion
preceding the statement of Theorem (3.1),
[TABLE]
From Theorem
(3.1), ∃ C1=C1(n,d,p,q,K1(λ))>0 such that for
all ϕ∈Bp(0,λ)
[TABLE]
[TABLE]
for all ψ∈Sp.
Using the Lipschitz continuity of σij and bi , the fact
that products of locally Lipschitz continuous functions are again
locally Lipschitz continuous, and the boundedness of ∂i:Sp→Sq,q≤p−21, we have
[TABLE]
and for q≤p−1
[TABLE]
for some constants C2=C2(n,d,p,q,λ,K(λ))>0 and C3=C3(n,d,p,q,λ,K(λ))>0 and for all ϕ,ψ∈Bp(0,λ). Similarly ∃ C4=C4(r,d,p,q,λ,K(λ))>0 such that
[TABLE]
We
now show that ∃ C5=C5(r,d,p,q,λ,K(λ),K1(λ))>0
[TABLE]
for ϕ,ψ∈Bp(0,λ). Consequently, the inequality (5) in the statement
now follows from equality (6) and the inequalities (7) - (12), with
the constant C in (5) given by C=2C1+C2+C3+C4+C5. To prove (12), we note that from the definition of A0 that
[TABLE]
Clearly it suffices to
show that ∃ C6′:=C6′(n,d,p,q,λ)>0 such
that for all j,k=1,…d and ϕ,ψ∈Bp(0,λ)⋂S
[TABLE]
where C6:=ψ∈Bp(0,λ)⋂Ssup ∥ψ∥q+1 C6′. But this
is an immediate consequence of the representation of the adjoint
∂j∗ of ∂j:S⊂Sq→Sq. Indeed it was shown in [5] that for
f∈S, we have,∂j∗f=−∂jf+Tjf, where Tj:Sq→Sq
is a bounded operator. In particular, for ϕ,ψ∈S,
[TABLE]
and (13) follows. This completes
the proof of Theorem (3.2).\hfill□
The following variant of the monotonicity inequality will be needed
in the proof of existence of translation invariant diffusions.
Before stating the result, we introduce some notation.
Let aij(ϕ),σij(ϕ),bi(ϕ) be as in Section 2.
For i=1,⋯,n we define Ai:Sp×Sp→Sq and L:Sp×Sp→Sq as follows :
[TABLE]
Note that
Ai(ϕ,ϕ)=Ai(ϕ) where the non-linear operator
Ai(ϕ) (of a single argument) is given by equation (2). Note
also that Ai(ϕ,ψ) is linear in the second variable and is
given in terms of the operator A0(⋅,⋅) defined in the
beginning of Section 3 as Ai(ϕ,ψ)=A0(σ(ϕ),ψ)⋅ei. Similar remarks hold for the
operator L(ϕ,ψ).
Theorem 3.3
Let p∈R,q≤p−1. Let σij(.),bi(.) be as
in Theorem (3.2). Then there exists a positive constant C1=C1(n,d,p,q,λ,K(λ),K1(λ)) such that
[TABLE]
for all ϕ1,ϕ2,ϕ3∈Bp(0,λ).
Proof: Let ϕi,i=1,2,3∈Bp(0,λ).
The left hand side of (16)
[TABLE]
The 1st term + the 3rd term in the right hand side of (17)
[TABLE]
where C1′=C1′(n,d,p,q,K1(λ)).
Using the Lipshitz continuity of the coefficients σij,bi
(see (4)), the 2nd term in right hand side of (17)
[TABLE]
where C3′′,C2′′,C3′,C2′ are positive constants depending only on
n,d,p,q,λ,K(λ) and K1(λ).
Similarly, the 4th term in right hand side of (17)
[TABLE]
where C4=C4(n,d,p,q,K(λ),K1(λ)).
Finally in the same manner as in the proof of Theorem (3.2), the 5th
term in right hand side of (17)
[TABLE]
where C5=C5(n,d,p,q,λ,K(λ),K1(λ)). The
proof of the theorem follows by summing up the terms in the RHS of
(17), using the above inequalities.\hfill□
4 Existence and Uniqueness of Solutions of SPDE’s.
Let p∈R. Let σij,bi:Sp→R i=1,⋯d,j=1,⋯,n be locally bounded functions on
Sp. Let (Bt) be a given n-dimensional Ft-
Brownian motion on (Ω,F,P) as in Section 2. Let
Ai,i=1,⋯,n and L be partial differential operators as
defined in equations (2) and (3). We now consider a stochastic
partial differential equation in S′ driven by the Brownian
motion (Bt) and ‘coefficients ’given by the differential
operators Ai,i=1,⋯,n and L defined above and initial
condition y∈Sp viz.
[TABLE]
where Y:Ω→Sp.
Note that if (Yt) is an Sp valued, locally bounded,
(Ft) adapted process then Ai(Ys),i=1,⋯,n and
L(Ys) are Sp−1 valued, adapted , locally bounded
processes and hence for i=1,⋯n, the stochastic integrals
∫0tAi(Ys) dBsi and ∫0tL(Ys) ds are well defined
Sp−1 valued, continuous Ft-adapted processes
and in addition, the former processes are Ft local
martingales . We then have the following definition of a ‘local’
strong solution of equation (18).
Definition 4.1
Let p∈R. Let σij,bi:Sp→R,i=1,⋯,d,j=1,⋯,n
be locally bounded functions, {Bt,Ft} a given standard
n-dimensional (Ft) Brownian motion and Y:Ω→Sp an F0 measurable random variable
independent of the filtration (FtB). Let δ be an
arbitrary state,viewed as an isolated point of Sp^:=Sp∪{δ}. By an Sp^ valued,
strong, (local) solution of equation (18), we mean a pair (Yt,η) where η:Ω→(0,∞] is an FtB-stopping time and (Yt) an
Sp^ valued (FtB)
adapted process such that
-
For all ω∈Ω, Y.(ω):[0,η(ω))→Sp is a continuous map
and
Yt(ω)=δ,t≥η(ω)
2. 2.
a.s. (P) the following equation holds in Sp−1 for 0≤t<η(ω),
[TABLE]
We note that equation (19) also holds in Sq for any q≤p−1. To prove the existence of solutions to equation (18) we
need a few well known facts. Let
(σˉij(s,ω)),(bˉi(s,ω)),j=1,⋯,n,i=1,⋯,d be locally bounded, Ft-adapted processes
and let Zt:=(Zt1,⋯,Ztd) be a d-dimensional (Ft) semi-martingale defined as follows :
[TABLE]
For p∈R, define the operator valued adapted
processes Lˉ(s,ω),Aˉi(s,ω):[0,∞)×Ω→L(Sp,Sp−1) i=1,⋯,n as follows : For ϕ∈Sp
[TABLE]
and for i=1,⋯,n,
[TABLE]
where aˉij(s,ω):=(σˉ(s,ω)σˉt(s,ω))ij,i,j=1,⋯,d. Let Y:Ω→Sp. Note that the
Sp valued process τZt(Y) has the Sp
valued trajectories t→τZt(ω)(Y(ω)).
We then have the following Lemma.
Lemma 4.2
Let p∈R . Let Yˉt:=τZt(Y) where Y:Ω→Sp an F0 measurable random variable independent of the filtration
(FtB), and (Zt) as above.
Suppose (σˉij(s,ω)),(bˉi(s,ω)),j=1,⋯,n,i=1,⋯,d are FtB-adapted locally
bounded processes. Then (Yˉt) is an Sp-valued
continuous FtB-adapted process which is the unique
solution of the following linear equation in Sq,q≤p−1 : almost surely,
[TABLE]
for every t≥0.
Let (Xt) be an Sp-valued progressively measurable process
which is uniformly bounded i.e. ∃K>0 such that ∥Xt(ω)∥p≤K, ∀(t,ω). Let
σˉij(s,ω):=σij(Xs(ω)),bˉi(s,ω):=bi(Xs(ω)) where σij,bi are
as in Defn.(4.1). Let (Zt),(Yˉt) be as defined above. Then
[TABLE]
where C=C(d,n,K,t) is a constant.
Proof: (a) The proof of the existence part of (a) for Y=y∈Sp fixed, is an immediate consequence of Itô’s
formula and we refer to [36] for the details. For Y arbitrary
but independent of the Brownian motion B, the result follows by a
conditioning argument. The proof of uniqueness
follows from the results in [20].
(b) It is sufficient to consider the case E∥Y∥p2<∞.
From the results of [37], we have
[TABLE]
where P(x) is a polynomial in x∈R with nonnegative
coefficients and degree m depending on ∣p∣. Now the result
follows by an application of the Burkholder-Davis-Gundy inequality
to each term in P(∣Zt∣), using the boundedness assumption on
(Xt) and the local boundedness of the coefficients σij
and bi,i=1,⋯,d,j=1,⋯,n. \hfill□
We now come to the existence and uniqueness of solutions to equation
(18). Recall that Bp(y,r) is the ball in Sp with centre
y and radius r>0.
Theorem 4.3
Let p∈R,q≤p−1. Let
σij,bi:Sp→R,
i=1,…,d, j=1,…n. Suppose that for every λ>0, ∃ K(λ)>0 such that
[TABLE]
for ϕ,ψ∈ Bp(0,λ)={η∈Sp:∥η∥p≤λ}. Then for every Y:Ω→Sp which is F0 measurable and independent of B
and for every r>0 there exists a strictly positive (FtB) stopping time ηr, and a Sp⋂Bq(Y,r)
valued,continuous, FtB-adapted process (Ytr)
satisfying equation (19) on [0,ηr), almost surely. If
Y1,Y2 are two F0 measurable Sp-valued random
variables with P{Y1=Y2}>0, then the corresponding solutions
(Yt1,r,η1,r),(Yt2,r,η2,r) satisfy :
η1,r=η2,r and Yt1,r=Yt2,r,0≤t<η1,r on the set {Y1=Y2}.
In particular, (19) has an S^p valued local, strong
solution. The solution is unique in the sense that if
(Yt1,η1) and (Yt2,η2) are any two solutions of
equation (18) with initial condition Y , then P[Yt1=Yt2,0≤t<η1∧η2]=1.
Proof: We will first prove uniqueness.
Uniqueness : Let (Yt1,η1) and (Yt2,η2) be any
two local solutions of equation (18) with initial condition Y∈Sp. Let λ>0. Let Yt=Yt1−Yt2. Let η0(λ):=inf{t:Yt1 or Yt2∈/Bp(0,λ)} and η≡η(λ):=η0(λ)∧η1∧η2. Then
by using the identity
[TABLE]
where {hk,q} is an
ortho-normal basis for Sq and expanding ⟨Yt,hk,q⟩q2 using Itô’s formula, we see that the
following equation holds a.s., for all 0≤t<η:
[TABLE]
where (Mt) is a continuous local martingale. Now
using inequality (5) of Theorem (3.2) , the Gronwall inequality and
a localisation argument (see for example [19]), we get for
each λ>0, almost surely , Yt1=Yt2, 0≤t<η(λ) . Letting λ↑∞ the result
follows.
Existence: To prove existence, we first consider the case
ω∈Ωsup∥Y(ω)∥p2<∞.
Recall the operator maps L(.,.) and Ai(.,.),i=1,⋯,n
defined in the paragraph prior to Lemma (4.2). Let for each k≥1,(Xsk) be an Sp-valued process. We define the operator
valued process Lk,Aik:[0,∞)×Ω→L(Sp,Sp−1),i=1,⋯,n, whose action
on ϕ∈Sp is given by,
[TABLE]
and for i=1,⋯,n,
[TABLE]
We define a sequence of FtB-adapted, Sp-valued
processes (Ytk), inductively, using operator valued processes
Lk(s,ω) and Aik(s,ω) as follows:
[TABLE]
where Y∈Sp is the given initial value of equation
(18). If (Ytk−1) is defined, then (Ytk) is defined as the
unique (FtB)-adapted solution of the linear equation
[TABLE]
where Lk(s,ω) and
Aik(s,ω) are defined as above, with Xsk(ω):=Ys∧ηk−1k−1(ω). Here, ηk−1 is an
FtB-stopping time, defined inductively, as follows : Let
r>0 be as in the statement of the theorem.
[TABLE]
j=1,…k−1. For notational convenience, in what follows, we often suppress
the dependence on r when there is no ambiguity. For notational
clarity, we note that σi denotes stopping times, whereas
σij denotes the coefficients in the operators L,A.
The existence and uniqueness of solutions of equation (20), is a
consequence of Lemma (4.2) with (Zt),(Yˉt) there taken to
be the processes (Ztk−1),(Ytk) respectively, where
(Ztk−1) is defined as follows :
[TABLE]
and
Ytk:=τZtk−1(Y).
Define η:=k→∞limηk. We note
that η≡ηr depends on r. We will show below that
η>0 almost surely. We now show that for each t≥0 the
sequence {Yt∧ηk} converges in L1(Ω→Sq) for q≤p−1. We have as in the proof
of uniqueness,
[TABLE]
where (Mtk) is a local martingale. Then using (14) and (15)
we have Lk(Ysk)=L(Ysk−1,Ysk),Aik(Ysk)=Ai(Yk−1,Ysk) with similar expressions for
Lk−1(Ysk−1) and Aik−1(Ysk−1) involving the
processes (Ysk−1),(Ysk−2) respectively. From Theorem (3.3)
and using a localisation argument, we can take expectations in the
above expression to get for some constant C1>0 and all k≥1,t>0,
[TABLE]
By the Gronwall inequality (21) now implies
[TABLE]
where K:=C(1+teCt) and C is some
positive constant.
Iterating the above inequality yields, for each t>0,
[TABLE]
where α=t0≤tsupE∥Yt0∧η1−Y∥q2<∞. It follows by the Cauchy-Schwartz inequality that for each
T>0 and 0≤t≤T,
[TABLE]
Define for each t,
[TABLE]
where the series in the right hand side converges in
L1([0,T]×Ω→Sq,dt dP),q≤p−1 and T>0 and defines an (Ft)-progressively
measurable Sq-valued process (Yt). We also note that,
for each t, Yt is an Sq valued random variable such that
E∥Yt−Yt∧ηk∥q→0. Note that for each
t≥0,k≥1, ∥Yt∧ηk−Y∥q≤r almost
surely and by passing to an almost sure convergent subsequence, we
also have, ∥Yt−Y∥q≤r almost surely. Denoting this
subsequence again by Ytk it follows by the bounded convergence
theorem that E∥Yt−Yt∧ηk∥q2→0 for
every t and moreover that E0∫t∥Ys−Ys∧ηk∥q2 ds→0.
We now wish to pass to the limit in equation (20). We note that by
the assumed continuity of the coefficients σij,bi;i=1,⋯d, j=1,…n and the continuity of ∂i:Sq→Sq−1/2,i=1,…,d, ∥Yt−Yt∧ηk∥q→0 almost surely for each t implies,
[TABLE]
for every s≤t∧η and i=1,⋯n, almost surely, where the
convergence takes place in Sq−1. Note also that there exists a
constant K1>0 such that for each 0≤s≤t, almost surely
[TABLE]
and a fortiori,
[TABLE]
holds for each 0≤s≤t, almost surely. It follows from the
above observations that
[TABLE]
for each t≥0, almost surely in Sq−1 and that
[TABLE]
for
each t≥0 in L2(Ω→Sq−1). Hence we can
pass to the limit in Sq−1 in equation (20) with t replaced by
t∧η, to get
[TABLE]
in
Sq−1 for every t>0, almost surely. It then follows, as a
consequence of Lemma (4.2), that if we define
[TABLE]
then (τZt(Y))
is a continuous Sp valued process that satisfies equation (19) in
Sq for any q≤p−1. We denote this process again by (Yt)
i.e. Yt:=τZt(Y). By its very construction the paths of
(Yt) are constant for t>η, and we can redefine this value to be δ to satisfy
definition (4.1) of a ‘ local solution’.
We now show that η>0 almost surely. Recall the stopping times σn,ηn
defined above. Suppose there exists a set A with P(A)>0 such
that if ω∈A then η(ω)=0. Then we claim that
∃ a subset A0⊂A with P(A0)>0 and a
subsequence {nk} such that for ω∈A0, σnk(ω)=ηnk(ω)<ηnk−1(ω),k≥2.
It is intuitively clear that such subsequences as described above
must exist. We give a proof for completeness: Fix t>0. We define a
sequence of integer valued random variables Ki for i≥0 as
follows:
[TABLE]
For i≥1,Ki:=min{j>Ki−1:ηj<ηj−1}.
Note that if ω∈A, then Ki(ω)<∞ for all i≥0 and σKi(ω)=ηKi(ω)<ηKi−1(ω). Further for every i≥0,
[TABLE]
Hence for every i≥1, we can choose an integer mi
satisfying
[TABLE]
Without loss of generality we can take mi>mi−1,i≥2. It
follows that for i sufficiently large
[TABLE]
It follows by the Borel-Cantelli lemma that
[TABLE]
In particular we have the almost sure equality
[TABLE]
where we define the set Bn for n≥1 as
[TABLE]
where the (countable) union in the last equality
is over jℓ≤mℓ for ℓ≥n, and j1<j2<…<jn−1<jn≤mn. Since Bn↑A and
P(A)>0, we choose n so that P(Bn)>0. Since each Bn is a
countable union of sets, each of which corresponds to a sequence {jℓ} and P(Bn)>0, ∃ a sequence {jℓ}, j1<j2<…<jk<… such that
[TABLE]
Define the set A0:={K1=j1,…Kℓ=jℓ,…}∩A. For ω∈A0, σjℓ(ω)=ηjℓ(ω)<ηjℓ−1(ω),ℓ≥1.
Note that P(A0)>0. This proves our claim. We will now work with
the subsequence (Ytnk)k≥1 which we will rename (Ytn) which then has the following property on A0 : If ω∈A0, then σn(ω)=ηn(ω)<ηn−1(ω),n≥2.
Fix t>0. Then we claim that, almost surely along a subsequence,
Yt∧ηkk→Yt∧η in Sp as
k→∞. This can be seen as follows. First note that
ηk defined above, decrease to η almost surely. Next,
observe that Ytk=τZtk−1(Y) where for each k≥1, the process (Ztk) is defined by
[TABLE]
Since the map x→τx(Y):Rd→Sp is continuous, to prove
the claim, it suffices to show that (Zt∧ηkk)
converges to the process (Zt∧η) in probability and
hence almost surely along a subsequence. To see this, note that we
can write
[TABLE]
By stopping (Ysk) at its exit from a suitably large ball
centered at Y∈Sp, we can show that the third and fourth
terms go to zero in probability. Hence using the (Lipschitz)
continuity of the maps σ,b and the convergence of Ys∧ηk to Ys in L2([0,t]×Ω→Sq)
it is easy to see that (Zt∧ηkk) converges to the
process (Zt∧η) in probability and hence almost surely
along a subsequence.
Thus our claim is proved and we have Yt∧ηkk→Yt∧η in Sp, almost surely , along a
subsequence. In particular it converges to Y in Sq on the set
A. For the rest of the proof, we will work with this subsequence,
which, abusing notation, we continue to denote by Yt∧ηkk. In particular we will now assume that the sum of integrals
in the right hand side of equation (20), evaluated at t∧ηk, goes to zero, which is the value of the sum of integrals
in the RHS of (22), almost surely on the set A.
On the other hand, for ω∈A0,σk(ω)=ηk(ω)<ηk−1(ω),n≥2; Hence for k≥1 such that ηk≤t, we have
[TABLE]
and consequently by the continuity of the process Yt∧ηkk in Sq,
[TABLE]
for ω∈A0 and k≥k0(ω) for some
k0(ω)≥1. But this leads to a contradiction to the fact
proved above that the RHS of (20) goes to zero in Sq, almost
surely on A and in particular on the set A0.
To complete the proof of the first part of the theorem, let r>0
and Y1,Y2 be as in the statement of the theorem with P{Y1=Y2}>0. First we assume Y1,Y2 are bounded. We will
denote with a superscript i,i=1,2, the various objects defined
in the construction of the solutions corresponding to the bounded
initial values Y1,Y2 respectively and suppress the dependence
on r, which is fixed. Firstly, we note that since ηi,j=ηi,j−1∧σi,j,i=1,2 and j≥2 and
since uniqueness holds for the linearized equation (20), an
induction argument shows that η1,j=η2,j and, almost
surely, Yt1,j=Yt2,j,0≤t<∞ on {Y1=Y2}. It follows that η1:=j→∞limη1,j=j→∞limη2,j=η2 on the set {Y1=Y2}.
To show that almost surely, Yt1=Yt2, 0≤t<η1 on
{Y1=Y2} we argue as follows. Define Y^ti:=I{Y1=Y2}Yti,i=1,2 and η′:=I{Y1=Y2}η1+I{Y1=Y2}∞. Then, using the quasi linearity of
(19), (Y^t1,η′),(Y^t2,η′) are two solutions of
(19) with initial value I{Y1=Y2}Y1. By uniqueness, our
claim follows.
The existence for the case of a general initial random variable Y
and a fixed r>0 can be reduced to the L2 case by considering
the initial conditions Yn:=YI{∥Y∥p≤n}, as follows.
Denote the solution corresponding to Yn, constructed above, by
(Ytn,ηn) where we have omitted the dependence on r.
Another solution with the same initial condition Yn is given by
(Yˉtn,ηˉn) where we define
[TABLE]
and
[TABLE]
Then by uniqueness we get that almost surely on the
set {∥Y∥p≤n}, Ytn=Yˉtn=Ytn+1, 0≤t<ηn, and ηn=ηn+1. We can now construct
the solution (Ytr,ηr) corresponding to the initial random
variable Y by piecing together the solutions on the sets
{∥Y∥p≤n} as follows : ηr(ω):=ηn(ω)
and Ytr(ω):=Ytn(ω),0≤t<ηn if ω∈{∥Y∥p≤n} and lies outside a suitable null set. That
(Ytr,ηr) is a solution of equation (19) follows from the
fact that (Ytn,ηn) solves (19) with initial condition Y on
the set {∥Y∥p≤n}. That it takes values in Bp(Y,r) on
[0,η) is also clear from the corresponding property for
(Ytn,ηn) on the set {∥Y∥p≤n}.
Let now Y1,Y2 be two F0B measurable, Sp
valued random variables and (Yt1,η1),(Yt2,η2) be the
corresponding solutions constructed above for some fixed r>0. To
show the claimed uniqueness on the set {Y1=Y2}, we define,
for n≥1, the processes
[TABLE]
where η′i:=ηi on {Y1=Y2,∥Y1∥p≤n}; =∞ otherwise. Then, Ytn,i solve (19)
on the interval [0,η′i), with initial values I{Y1=Y2,∥Y1∥p≤n}Yi,i=1,2 respectively; and by uniqueness
for the case of bounded initial random variables discussed above, it
follows that, almost surely, η1=η2 and Yt1=Ytn,1=Ytn,2=Yt2,0≤t<η1 on the set
{Y1=Y2,∥Y1∥p≤n}. Letting n→∞
the uniqueness claim follows. This completes the proof of the
theorem. □
Remark 4.4
A simpler proof of existence can be provided using
the finite dimensional existence results as in Theorem 2.1. In
effect, we fix the initial value y and we define
σˉij(x):=σij(τxy),bˉi(x):=bi(τxy). Then if σij,bi are Lipschitz in Sq,q≤p−21, then one can show (using duality
and the mean value theorem applied to <τxy−y,ψ>,ψ∈S )that σˉij,bˉi are locally Lipshitz on
Rd. If (Zt,η) is the solution of equation (1) with
x=0 then using Ito^’s formula for translations, it is easy
to see that Yt:=τZty,t<η, solves equation (19).
However this proof does not work when, for example, we replace the
operator L in equations (3) and (19) with a perturbation of L
viz. L+cI, where I is the identity and c∈R ( see
also Example 6 of Section 6). **
Remark 4.5
Translation invariance also applies to solutions of an
evolution equation which is a first order quasi linear PDE i.e.
these solutions are translates of the initial condition by the
solution of an appropriate ‘characteristic’ ODE. This follows on
setting the diffusion coefficients in the above calculations to be
equal to zero. These first order systems may also be viewed as the
‘zero noise’ limit of stochastic second order system, a topic of
considerable interest in the last three or four decades (see
[17]), and also volumes 3 & 4 of [7] for the connection
with large deviation theory.**
5 The Strong Markov Property :
In this section we show that the local solution of equation (19)
obtained in Theorem (4.3) can be extended to a maximal interval
[0,η) for a given F0-measurable Y (Theorem (5.3)).
We then show that the solutions (Yt(y),η) obtained when Y≡y∈Sp have a jointly measurable version
(Y(t,ω,y),η(ω,y)) in (t,ω,y) and that the
solutions for arbitrary Y can be represented as Y(t,ω,Y),0≤t<η(ω,y), (Proposition (5.4) and Theorem (5.5)). The
strong Markov property (Theorem (5.7))is then proved as a
consequence of uniqueness in law, which in turn follows from
pathwise uniqueness by a Yamada-Watanabe type argument (Theorem
(5.6)).
Consider now the solution (Ytr,ηr) constructed in Theorem
(4.3) for r>0 and for some F0-measurable random
variable Y:Ω→Sp. From the definition of
ηr in the proof of Theorem (4.3), it follows that, r1<r2
implies ηr1(ω)≤ηr2(ω). Let
η(ω):=r↑∞limηr(ω).
Then by pathwise uniqueness of solutions we have Ytr1=Ytr2,0≤t<ηr1, almost surely. Let rn↑∞. Then ηrn(ω)↑η(ω) ∀ ω. Let Ωn satisfy P(Ωn)=0 and for ω∈/Ωn,Ytrn(ω)=Ytrn+1(ω),0≤t<ηrn. Let Ω0:=n=1⋃∞Ωn. Define, for ω∈/Ω0 and 0≤t<η(ω)
[TABLE]
and let
[TABLE]
For ω∈Ω0 redefine η(ω):=0 and define
[TABLE]
We note the following ‘
maximality’ property of the solution (Yt,η).
Proposition 5.1
Let Y∈Sp,q≤p−1. Then,
a.s., t↑η(ω)lim∥Yt(ω)∥q=∞ on ω∈{η<∞}. In
particular, a.s., t↑η(ω)lim∥Yt(ω)∥p=∞ on ω∈{η<∞}.
Proof: Note that {η<∞}=n⋃{η<n}. It suffices to show that r→∞lim∥Yηr∥=∞ a.s. on {η<∞}. Recall
the approximations (Ytk)≡(Ytk,r) to the solutions
(Ytr,ηr) constructed in the proof of Theorem (4.3). For
fixed r>0 we know that Yt∧ηrk→Yt∧ηrr in L2(Ω,Sq) and in particular for every ϵ>0
[TABLE]
as k→∞. Further if ηk,r are the
approximations to ηr constructed in the proof of Theorem (4.3)
we have for each k≥1,
[TABLE]
Since ηk,r↓ηr and ∥Yuk∥q≤r for t∧ηr≤u≤t∧ηk,r, the first term goes
to zero by the bounded convergence theorem almost surely and the
second term goes to zero in probability as k→∞.
Thus for every ϵ>0,
[TABLE]
It follows from the above that for every ϵ>0,
[TABLE]
Thus for any t>0 and passing to a subsequence {ki} we have
a.s. on {ηr<t},
[TABLE]
We can argue (as in the proof that ηr>0 a.s, in the proof of
Theorem (4.3)), that by passing to a further subsequence, that ηki,r=σki,r and in particular that on
{ηr<t}
[TABLE]
It follows that a.s. on {ηr<t},∥Yηr−Y∥q=r.
Now we take rk↑∞. Then {η<∞}=n=1⋃∞k=1⋂∞{ηrk<n}. In particular a.s. on {η<∞},
[TABLE]
. Since q<p the result
follows. \hfill□
Proposition 5.2
Let E∥Y∥p2<∞ and σij,bi:Sp→R be bounded i.e. ∃ K>0 such that
[TABLE]
for all ϕ∈Sp,i=1,…d,j=1…n. Then η=∞ a.s.
Proof: It suffices to show that for every t>0,P(ηr≤t)→0 as r↑∞. As in the proof of the
previous Proposition, {ηr<t}⊂{ηr<t,∥Yηr−Y∥q=r}⊂{∥Yηr∧t−Y∥q≥r}. Hence,
[TABLE]
Further by the boundedness assumptions on σij,bi and the
monotonicity inequality given in Theorem (3.1) ( applied with ϕ=ψ=Ys), we have
[TABLE]
where C>0 is a constant depending only on r,d,p and K. Using
Gronwal’s inequality we get
[TABLE]
Hence dividing by r2 and letting r↑∞ in the above inequality
we conclude that P(ηr≤t)→0 as r↑∞ for every t>0. \hfill□
Theorem 5.3
Let p∈R,q≤p−1. Let
σij,bi:Sp→R,
i=1,…,d, j=1,…n be as in Theorem (4.3). Then for
every Y:Ω→Sp which is F0
measurable and independent of B, equation (18) has a unique
S^p valued strong (local) solution (Yt,η).
Further η>0 is maximal in the sense that, almost surely,
[TABLE]
Finally, if we define
(Zt) as
[TABLE]
then (Zt)
is a continuous, (FtY)-adapted, Rd-valued
process such that almost surely,
[TABLE]
Proof : The proof follows by ‘patching up’ the solutions
obtained in Theorem (4.3), as described in the beginning of Section
5. This gives us the pair (Yt,η). That it solves equation (19)
follows from Theorem (4.3) and the fact that by construction, for
any r>0, Yt=Ytr,0≤t<ηr. The maximality of the
solution follows from Proposition (5.1). We note that Yt=τZt(Y) follows from Lemma (4.2). \hfill□
To formulate the strong Markov property, we consider the solution
(Yt,η) when the initial value Y is a constant Y≡y∈Sp and denote the corresponding solution by (Yt(y),ηy) or (Yt(ω,y),ηy).
Recall that S^p=Sp∪{δ}. We define
the σ-field B(S^p) on S^p by
A^∈B(S^p) iff A^=A∪{δ}
for some A∈B(Sp). A measurable function f:Sp→R is extended to S^p
by defining f(δ)=0. The resulting extension will also be
denoted by f.
For y∈Sp, we define
[TABLE]
In other words, η(ω)=0 for such y. We now construct
versions of the solution (Yt(y),ηy) constructed in Theorem
(5.3), with initial value Y≡y∈Sp, which are
jointly measurable in (t,ω,y). In the two Propositions below,
we need the approximations (Ytk),k≥1 for initial values Y≡y, constructed in the proof of Theorem (4.3), which we now
denote by (Ytk(y)) or as (Ytr,k(y)), whenever the
dependence on the domain Bp(y,r) needs to be made explicit. The
proofs of the following two results (Proposition (5.4) and Theorem
(5.5) are given in the Appendix).
Proposition 5.4
a). There exists a map Y~:[0,∞)×Ω×Sp→S^p which is B[0,∞)⊗F∞B⊗B(Sp) / B(S^p) measurable and satisfies
[TABLE]
b). There exists a map η~:Ω×Sp→[0,∞] which is F∞B⊗B(Sp) / B[0,∞] measurable and satisfies
[TABLE]
The next result shows that the solution of the SPDE (19) with an
initial random variable Y0 as the composition of Y0 with the
solutions starting at y∈Sp. So let Y0:Ω→Sp be F0-measurable. Define
[TABLE]
Then we have the following
result.
Theorem 5.5
Let (Yt,η) be the solution of equation (19)
with initial r.v. Y0 independent of (Bt). Then for each t≥0, we have η=η~ a.s. and
[TABLE]
We now prove uniqueness in law for equation (19) required to prove
the strong Markov property. This follows from the Yamada-Watanabe
result for SPDE’s of the type (19), which we now state as the next
theorem. We need some preliminaries to deal with the law of
explosive solutions.
We first construct an appropriate path space i.e. a measurable space
(C,C) such that if (Yt,η) is a maximal solution on some
probability space then almost surely, the paths Y⋅∧η(ω) belong to C and the map ω→Y⋅∧η(ω) is measurable. It is clear that the
law of (Y⋅∧η) is essentially determined on
[0,η) where the paths are continuous. However, we need to
distinguish between the cases η<∞ and η=∞.
Further, although our initial conditions and the paths of the
corresponding solutions lie in Sp, we will consider them
as paths in Sq, where the equation holds.
Thus let p∈R,q≤p−1. Let y:[0,∞)→S^q. All such maps that we consider will be B([0,∞))/B(S^q) measurable. Define
ηq(y):=inf{s>0:y(s)∈/Sq}.
Let
[TABLE]
[TABLE]
Define C:=C1⋃C2. For y∈C, define for r>0,τr(y):=inf{t>0:∥y(t)−y(0)∥q>r}. Then y∈C implies y∈C0 and by
continuity of y on [0,ηq(y)), we have τr(y)<∞, and ηq(y)=r↑∞limτr(y). Note that C([0,∞),Sq)={y∈C:ηq(y)=∞}. We now define a sigma field C on C
via the maps Kr:C→C([0,∞),Sq),r≥0,Kr(y):=y(⋅∧τr) as follows.
[TABLE]
Let (Yt,η) be a maximal solution of equation (19) with initial
value Y∈Sp and Brownian motion (Bt), obtained in
Theorem (5.3) on some probability space (Ω,F,P) . Then
recall that this solution is obtained by pasting together the
solutions (Ytr,ηr),r>0, obtained in Theorem (4.3). Then
we have a map Y^:Ω→S^q,
[TABLE]
By Proposition (5.1) it follows
that almost surely, ηq(Y^(ω))=η(ω) and
hence Y^(ω)∈C almost surely. We can redefine Y^
on a null set so that Y^:Ω→C⊂S^q. Since ηr(ω)≤τr(Y^(ω))=:τr(ω)<ηq(Y^(ω))=η(ω) and since
ηr(ω)↑η(ω), we have for a fixed r0>0,
[TABLE]
Since the right hand
side is a limit of measurable maps, we conclude that for each r0≥0, the maps ω→Kr0(Y^(ω)) is
FB measurable and hence ω→Y^(ω) is FB/C measurable.
Let (Yti,Bti,ηi) i=1,2 be two ‘maximal’ solutions of
(19) viz.
[TABLE]
adapted to the n dimensional Fti-Brownian motions (Bti),i=1,2,
with F0i measurable initial values Yi:Ωi→Sp, independent of (Bti), on possibly
different probability spaces (Ωi,Fi,Fti,Pi). Equality in law between two random variables X1,X2
will be denoted by ”X1≜X2”. For r>0 we will denote by ηi,r,r>0, the stopping times upto which
the solution Yti,r lies in a ball of radius r around Yi
in Sq,q≤p−1 (Theorem (4.3)) and ηi=r→∞limηi,r, the explosion time.
Let Pi,i=1,2 be the laws of (Yti) on (C,C). Let Wn:=C([0,∞),Rn).
Theorem 5.6
If Y1≜Y2, then P1=P2.
Proof: Let Pi,y,i=1,2 be the law of the solutions
(Yti,y,ηi,y) of equation (19) corresponding to Yi≡y∈Sq on (C,C). Using Theorem (5.5) and
the independence of Yi and (Bti) and the definition of C we have
[TABLE]
it suffices to show that
for all y∈Sp,P1,y=P2,y on C. From
the definition of C, it suffices to show that for every
r0>0, the laws of (Yti,y,ηi,y) composed with the
map Kr0:C→C([0,∞),Sq) i.e. the laws
of Kr0(Y^i,y) agree on B(C([0,∞),Sq)). We fix y∈Sp. In what follows, we drop the
explicit dependence on y in our notation. Let Pri(A):=P(Y⋅∧ηi,ri,r∈A) where A∈B(C([0,∞),Sq)). Since, as was observed above,
Kr0(Y^i) is the almost sure limit (Yt∧τr0∧ηi,ri,r) as r↑∞, it suffices
to show that for every r>0,Pr1=Pr2 on B(C([0,∞),Sq)).
The proof is basically the same as in the proof of the finite
dimensional Yamada-Watanabe result. In our case the finite
dimensional diffusions are replaced by the infinite dimensional
processes (Yti,r,ηi,r),i=1,2 with a fixed initial
value y. We follow the proof in [22] (Chapter IV, Theorem
(1.1) and its corollary), the only difference in our case being that
the space C([0,∞),Rd) is replaced by the space
C([0,∞),Sq), which is again a Polish space. It
suffices to show then that Pr1(A)=Pr2(A),A∈B(C([0,∞),Sq)). Let Qri(A×B):=P{Y⋅∧η1,ri∈A,W⋅∈B}, A∈B(C([0,∞),Sq)), B∈B(Wn). Let P0 be the
Wiener measure on Wn.
Let Qri(ω,A) be a disintegration of Qri w.r.t.
P0, i.e. for i=1,2,
[TABLE]
A∈B(C([0,∞),Sq)) and B∈B(Wn).
Then as in the proof in [22],Theorem (1.1), using pathwise
uniqueness, there exists a measurable map Fr:W3→C([0,∞),Sq) such that
[TABLE]
for i=1,2. In particular it follows that
[TABLE]
\hfill□
Having defined measurable versions of (Yt(y),ηy) we can now
define the transition probability function P(t,y,A) for t≥0,y∈S^p and A∈B(S^p) in the
usual way:
[TABLE]
From Proposition (5.4), it follows that for fixed A∈B(S^p), the map (t,y)→P(t,y,A) is
jointly measurable and a probability measure for fixed t≥0
and y∈S^p. For y∈Sp we will write
[TABLE]
where
[TABLE]
for 0≤t<ηy. We can
then formulate the strong Markov property of the process (Yt(y))
as follows.
Theorem 5.7
Let T:Ω→[0,∞] be an (FtY) stopping time. Then for each
y∈Sp, a.s. on {T<∞},
[TABLE]
for t≥0,A∈B(S^p).
Proof: We first consider the case y∈Sp and A⊆Sp. Let f:S^p→R be a bounded measurable function, f(∂)=0. Then,
with y∈Sp and η:=ηy
[TABLE]
where
[TABLE]
We note that on {T<∞}
[TABLE]
where (Y^t0(z),η^z)
satisfies equation (26)(with respect to P( . ∣T<∞)) with
y replaced by z and with (Bt) replaced by the Brownian motion
(B^t):=(I(T<∞)(BT+t−BT)); and η^z is
the explosion time for the process Y^t(z):=z+Y^t0(z)
which by its maximality, satisfies η≡ηy=T+η^YT(y) on {ηy>T}. Since the latter
Brownian motion is independent of FT we have (using
Theorem (5.5)),
[TABLE]
on the set {T<η}, and where we have used the uniqueness in
law for equation (26), which follows from the previous theorem.
We now consider the case y=δ. Then both sides of the
equation in the statement of the theorem reduce to IA(δ).
Next let y∈Sp,A={δ}. we have
[TABLE]
where we
have used the independence of the Brownian motion (B^t) and
FTY in the second equality. This completes the proof.
\hfill□
It is clear from the relation Yt=τZt(y),0≤t<ηy, with (Zt) as in equation(23) that the path properties of
the processes (Yt) and (Zt) are closely related, although they
live in different spaces. In particular, as already observed in
Proposition (3.12) of [39] corresponding to the case where
σ,b are given by linear functionals on Sp, the
explosions of (Zt) as t→∞ are related to the
convergence of Yt to zero in the weak topology of S′ and
this correspondence is pathwise. It is easy to see that the result
of Proposition (3.12) of [39] extends to the more general
framework of Theorem (5.3) above. We then have the following result.
Proposition 5.8
Let σi,j,bi be as in Theorem 4.3, with Y≡y=0∈Sp. Let (Yt(y),ηy) be the unique maximal
solution of equation(19),and (Zt(y)) be given by equation (23)
with Yt(y)=τZt(y)(y),0≤t<ηy. Fix ω∈Ω. Then, Zt(ω,y)→∞ as t→ηy(ω) whenever Yt(ω,y)→0
weakly in S′ as t→ηy(ω).
Conversely, suppose one of the following two conditions is satisfied
viz.
-
y∈Lp(Rd),p≥1.**
2. 2.
y* has compact support.*
Then, Yt(ω,y)→0 weakly in S′
whenever Zt(ω,y)→∞ as t→ηy(ω).
Proof : The proof is the same as in Proposition (3.12) of
[39]. The proof for the case y∈Lp,p≥1,p=2
is also the same as the case p=2 with some obvious changes.
\hfill□
Remark 5.9
Note that when ηy<∞ and Zt(ω,y)→∞
then by the above Proposition, Yt(ω,y)→0 weakly
in S′ as t→ηy while by Proposition (5.1),
∥Yt(y)∥p→∞.
6 Some Examples
Our main existence and uniqueness result
viz. Theorem (5.3), applies to a number of different situations. In
this section we give some examples of these applications. In what
follows we use the fact that if p>4d then δx∈S−p (see [38],
Theorem (4.1)) and for such p we also note that ϕ∈S−p is a continuous function.
Example 1
Let p>4d+1. Then note that Sp⊂C2(Rd), the space of two times continuously differentiable
functions on Rd (see Theorem 4.1, [38]). Let
(Yt)0≤t<η be the unique Sp- valued strong
solution of equation (18) with initial condition y∈Sp,
given by Theorem (5.3). Then almost surely, for t<η, it is
given by a C2(Rd) function, say, x→Yt(ω,x) and we also have
[TABLE]
for ∣α∣≤2. In particular acting
on both sides of (19) by δx∈S−p we get for
each t<η,x∈Rd, almost surely,
[TABLE]
where the integrands in the RHS of the above equation are well
defined processes for each x and the stochastic integrals are well
defined. Since Yt=τZt(y),t<η with (Zt) as in
(23), in particular Yt(x)=y(x−Zt),x∈Rd. Since y∈Sp,p>4d+1, it is in C2(Rd)
and the Itô formula applied to y(x−Zt) also yields the RHS of
(27). Thus in this case, (Yt)t≥0≡{Yt(x):t≥0,x∈Rd} gives the unique classical solution of the
SPDE (18) when p>4d+1.
We also note that the Fourier transform f∈Sp→f^∈Sp is a unitary map on the complexified
Hermite-Sobolev spaces Sp(C) (see [45]).
Hence we get from the above that the Fourier transform Y^t of
Yt is given as Y^t=y^ei⟨⋅,Zt⟩, where ⟨.,.⟩ is the inner product in
Rd and the RHS represents the product of the tempered
distribution y^, the Fourier transform of y, with the
bounded C∞ function x→ei⟨x,Zt⟩. Note that for each x∈Rd, Y^t(x)
is a process and it is easily seen that it satisfies a linear SDE
obtained by the Fourier transform of equation (19).**
Example 2
The connection between solutions of equation (19) and the
solutions of the finite dimensional SDE (1) was shown in [39].
Let (Zt) be as in equation (23). Then it follows as in [39]
that the process (Xtx) defined by Xtx:=x+Zt(τxy), 0≤t<η , solves the equation
[TABLE]
where for z∈Rd,σˉij(z):=σij(τz(y)), bˉi(z):=bi(τz(y)),j=1,⋯,n,i=1,⋯,d and y∈Sp acts as a fixed parameter.
Special cases arise when σij,bi:Sp→R are continuous linear functionals on Sp i.e.
they are given by elements in S−p and consequently
[TABLE]
where ⟨⋅,⋅⟩ denotes duality between S−p and Sp. Note that when σij,bi,y
are functions in L2(Rd), then
[TABLE]
where ∗ denotes convolution and y~(z):=y(−z) and
similarly bˉi(z)=bi∗y~(z). When p<−4d, then we can take y=δ0∈Sp,σij,bi∈S−p⊂C(Rd), the
space of real valued continuous functions on Rd, and
σˉij(z)=σij(z), bˉi(z)=bi(z). **
Remark 6.1
The weak existence of solutions to the Itô SDE (28)
can be combined with the pathwise uniqueness of solutions to
equation (18) when σ,b are in Sp to yield pathwise
unique solutions of (28). The weak existence is obtained whenever
the coefficients σˉij, bˉi in (28) are
bounded and continuous. On the other hand any two solutions of (28)
with the same Brownian motion gives rise via Lemma (4.2), to
corresponding solutions of (18) forcing the former solutions to be
the same (see Theorem (3.3) of [40]).**
We can vary the construction in the above example to get strong
solutions in the case of Lipschitz continuous functions. We do this
in the following one dimensional example. The general finite
dimensional case can be handled by considering finitely many
equations like (19).
Example 3
Let d=1,p>41 and σj=σ1j,b=b1:R→R,i=1,j=1,⋯,n be Lipschitz
functions. For x∈R define b^(x,.),σ^j(x,⋅):Sp→R by σ^j(x,y):=σj(y(x)) and b^(x,y):=b(y(x)). Note that under the assumptions on p, the elements of
Sp are continuous functions. Then we note that for
y1,y2∈Sp,
[TABLE]
with a similar inequality for b and where the constant K depends
on x. Let L and A be the operators as in equations (3) and (2)
with σj,b replaced by σ^j(x,⋅) and
b^(x,⋅). Then for any fixed initial value y0∈Sp, equation (18) has a unique Sp valued strong
solution which we denote by (Yt(x,y0)). We then have
Yt(x,y0):=τZt(x,y0)(y0) where (Zt(x,y0)) is
given by (23) with σ and b there replaced with σ^j(x,⋅),b^(x,⋅) respectively and y0∈Sp, as defined above. Then it follows as in Example 2 that
(Zt(x,y0)) solves the ordinary SDE
[TABLE]
let σˉj(z):=σ^j(x,τz(y0))=σj((τzy0)(x))=σ(y0(x−z)) and a similar
expression for bˉ(z). Then Xt(x):=x−Zt(x,y0) solves
[TABLE]
Moreover, in a manner similar to the case of uniqueness discussed
in the Remark (6.1) above, the uniqueness of solutions of (18)
implies that the solution of (30) is unique : Any two (Bt)
adapted solutions of equation (30) will give rise to two solutions
of equation (29), which in turn (via Lemma (4.2)) gives rise to two
solutions of (18). The same arguments also imply that there is local
uniqueness in equation (30), upto a stopping time i.e. local
uniqueness upto a stopping time in equation (18) implies local
uniqueness in equation (30) upto a stopping time. If now we consider
a sequence of elements y0k∈Sp,k≥1, satisfying
[TABLE]
then a localisation argument
implies that the corresponding solutions (Xtk(x)) satisfies
Xtk(x)=Xtk+1(x),t≤τk where τk is the exit
time of (Xtk(x)) from the ball {z:∣z−x∣≤k}. One can
then patch up the solutions (Xtk(x)),k≥1 to obtain the
solution of the equation
[TABLE]
when the coefficients
σj,b,j=1,⋯,n, are given Lipschitz continuous
functions. **
Example 4
In this example we consider martingale problems in the sense of
Stroock and Varadhan, associated with a second order differential
operator Lˉ with coefficients σˉij and bˉi,i,j=1,…d which are bounded and continuous functions on
Rd. If in addition they belong to Sp,p>4d+1 we can solve the SDE (18) with σij and
bi given by the linear functionals on S−p
corresponding to σˉij and bˉi,i,j=1,…d
and initial condition δx∈S−p. In this
situation we have indeed a unique strong solution to the Ito SDE
(1). In case we know only that σˉij and bˉi,i,j=1,…d are bounded and continuous, then since they are
tempered distributions, there exists p>0 such that they belong
to S−p. In this case, we still have strong solutions of
(18) with η=∞ (Proposition (5.2)) for initial conditions
y∈Sp,p>0 and of course δx∈/Sp. Below we show that when yn→δx
weakly in S′ and Yt(yn)=τZn(t)(yn)
are the solutions of (18), then the laws {Pxn} of the
processes (x+Zn(t)), converge weakly to Px, the solution of
the martingale problem for Lˉ starting at x, provided the
latter is well posed.
We have z∈Rd,
[TABLE]
On the other hand consider the SPDE equation (18) with coefficients
σij,bi:Sp→R given by
σij(ϕ)=⟨σˉij,ϕ⟩ and a similar expression for bˉi(ϕ),ϕ∈Sp. Let yn∈Sp∩C(Rd),yn→δx weakly for a fixed x∈Rd.
Let (Ytn),Ytn:=Yt(yn) denote the unique Sp
valued solution to (18) with initial condition Y0n=yn. Then
Ytn=τZtn(yn) where (Ztn) comes from equation (23)
with (Yt) replaced by (Ytn). Let σˉijn(z):=⟨σˉij,τz(yn)⟩,bˉin(z):=⟨bˉi,τz(yn)⟩. Then
[TABLE]
where we have used the notation σˉn=(σˉijn)
and bˉn=(bˉin) for the diffusion and drift
coefficients respectively. Let Pn be the law of (Ztn) on
C([0,∞),Rd) and Zt(ω):=ω(t) the
coordinate process. For φ∈S let
[TABLE]
Let s<t and G be a bounded,continuous and Fs-measurable function of the path ω∈C([0,∞),Rd) depending on finitely many time
coordinates. For f∈S, we have by Itô’s formula
[TABLE]
Suppose now that Pn→P weakly on
C([0,∞),Rd). Let Lˉx be the operator in
(32) wherein σˉij(z),bˉi(z) are replaced with
σˉij(x+z),bˉi(x+z), x∈Rd fixed.
We then have :
[TABLE]
To see this, first note that the integrand is a bounded continuous
function on C([0,∞),Rd). Further, as n→∞, we have σˉijn(z)→σˉij(x+z),bˉin(z)→bˉi(x+z).
Moreover,
[TABLE]
Our claim now follows by using the Skorokhod mapping theorem and the
bounded convergence theorem. In particular, it follows that any weak
limit P of the sequence {Pn} solves the martingale problem
for Lˉx starting at zero. We then have the following
theorem.**
Theorem 6.2
Suppose the martingale problem for Lˉ starting at x has a
unique solution Px. Let yn∈S−p∩C(Rd),yn→δx weakly. Let (Ztn) be
as above and let Pxn be the law of (x+Ztn). Then Pxn→Px weakly.
Proof: Replacing f by τ−xf we see that if P is any
weak limit of the family {Pn,n≥1} where Pn is the law
of (Ztn), then under P,Xtx:=x+Zt solves the martingale
problem for Lˉ starting at x and hence the law of (Xtx)
must be Px. The tightness of the laws {Pn,n≥1} viz.
for every ϵ>0,T>0
[TABLE]
and hence the tightness of {Pxn}, follows easily from Doob’s
maximal inequality, the Burkholder-Davis-Gundy inequalities and the
uniform bounds in (35). □
Example 5
In this example we consider the non-linear evolution equation
[TABLE]
Here y∈Sp for some p∈R and L:Sp→Sp−1 is given by equation (3). By a solution
we mean a pair (Yt,η) where η>0 and (Yt) is a
continuous function t→Yt:[0,η)→Sp
satisfying the following equation in Sp−1
[TABLE]
for 0≤t<η. Suppose (Yt)0≤t<η is an
Sp valued solution. Define the time dependent, linear
operators Aˉt,Lˉt:Sp→Sp−1 as follows:
[TABLE]
where h=(h1,⋯,hn). Note that the coefficients are now
deterministic but time dependent. Define the Rd-valued
process (Zt)0≤t<η by
[TABLE]
for 0≤t<η. Since the integrands are deterministic
(Zt) is a Gaussian process. Let Yˉt:=τZt(y),0≤t<η. Then (Yˉt)0≤t<η is the unique
Sp-valued solution of the equation
[TABLE]
Let φ(t):=EYˉt, 0≤t<η where we note that
E∥Yˉt∥p<∞. Then φ(t) satisfies the linear
evolution equation
[TABLE]
in the interval 0≤t<η. Since Lˉt has constant (in
space) coefficients it satisfies the monotonicity inequality and
hence equation (38) has a unique Sp-valued solution. Hence
we have the following stochastic representation of solutions of
equation (36).**
Theorem 6.3
Let p∈R,y∈Sp, and let σij,bi:Sp→R be
bounded and measurable. Let (Yt)0≤t<η be an Sp-valued solution of equation (36). Then we have,
[TABLE]
where pZt is the density of Zt and ∗-denotes
convolution.
Example 6
The previous example maybe generalised. Consider the following
equation viz.
[TABLE]
of which (37) becomes a special case when there is no dependence
on y in the operator L. However we will make a departure from
the L in (37) by requiring L to act on y as a partial
differential operator with the coefficients σij,bi
depending on Ys in the right hand side above. In other words,
[TABLE]
where y1,y2∈S−p,p∈R. If μ(dz) is
a probability measure, then μ∈S−p,p>4d and we can define the non-linear convolution L(⋅,y2)∘μ)(y1) ([39], Section 5, where the
notation in definition (5.1) is slightly different and the
coefficients σij and bi do not depend on y2) as
[TABLE]
whenever the integral exists as a Bochner
integral in S−p. An interesting situation arises when
the measure μ arises as the marginals of a stochastic process
(Zt). Let {μs(dz),s≥0} be the corresponding family
of probability measures. Consider the case when
σij(⋅,⋅),bi(⋅,⋅) are uniformly bounded
and (Zt) satisfies the equation
[TABLE]
where μt(dz):=P(Zt∈dz) is the law of Zt and y∈S−p. Then applying
Ito^’s formula and taking expectations we get that Yt:=EτZty=y∘μt=:ψ(t,y) satisfies the non linear
evolution equation
[TABLE]
with ψ(t,L(y,y2)):=L(⋅,y2)∘μt(y) and y2=ψ(t,y). Equation (37) becomes a special case of (42) when
σij(y1,y2),bi(y1,y2) are independent of y1. When
we consider y=δx where x∈Rd is fixed then
Xt:=x+Zt and we get the Mckean-Vlasov equation from (41).**
Example 7
Let p>4d. We now consider the Feynman-Kac formula
for the solution of the equation
[TABLE]
where
[TABLE]
Here we assume f,V,σˉij,bˉi are given functions
in Sp. Then we define L as in equation (3) with
coefficients σij(⋅) and bi(⋅) given via the
duality between Sp and S−p as
σij(y)=⟨y,σˉij⟩,bi(y)=⟨y,bˉi⟩ where σˉij and bˉi are as
above and y∈S−p.
Denoting by (Xtx) the diffusion corresponding to Lˉ and by
(Yt(y)),y=δx the corresponding lift on Sp
satisfying equation (18), it is easy to see that the solution
u(t,x) arises from a transformation on path space
C([0,∞),S−p) viz.
[TABLE]
where c(s,y):=−⟨ys,V⟩, y∈C([0,∞),S−p). In particular, since Yt(δx)=δXtx, c(s,Y)=V(Xsx). For ease of calculations, we assume η=∞,a.s. Next, with u(t,x):=PtVf(x) where (PtV) is
the Feynman-Kac semi-group, we have
[TABLE]
We can show that the process (Y^t) satisfies an SPDE with
time dependent coefficients L^(s,y),A^i(s,y),i=1,⋯n,s≥0,y∈C([0,∞),S−p) given in
the form
[TABLE]
Here the coefficients σ^ij,b^i,c^ are
induced on [0,∞)×C([0,∞),S−p) by the
coefficients σij,bi of L,Ai appearing in L and the
transformation y→yc given by the unique solution of
the equation
[TABLE]
on the
path space C([0,∞),S−p) and satisfying σ^ij(t,y)=σij(ytc), b^i(t,y)=bi(ytc), c^(t,y)=c(t,yc).
It is then easy to see using integration by parts and the fact that
(Yt) solves (18), that, Y^t satisfies the SPDE
[TABLE]
The uniqueness of solutions of the above SPDE can be proved using
the uniqueness of the solutions of the equation ytc=yt e−0∫tc(s,yc)ds and the ‘invertibility’of the
map y→yc. The details can be seen in [41].**
7 Conclusion
Translation invariance also appears to be
a reflection of a possibly more basic,‘duality’ relation between the
finite dimensional SDE and the corresponding SPDE. Let p>0,f∈Sp,y∈S−p.
Let σij,bi,(Yt(y),η) be as in Theorem (5.3), with Y≡y.
We consider the case η=∞. Let (Zt(y)) be as in
equation (23). We observe the following duality relation between Yt(y) and Zt(y) viz.
[TABLE]
whenever the
relevant expectations are finite.
Finally we note that in the model we have introduced in this paper,
it becomes meaningful to talk about diffusions with coefficients
σ and b, in the state y, for any tempered
distribution y. The ’state y’ becomes an initial state for the
SPDE, but in the context of the SDE, allows for representation of
more complex initial states than just y=δx. The
distribution y is more intutively, thought off as an initial
distribution of the mass of the solvent particles in the diffusion
model. An interpretation of ‘translation invariance’ in the case of
non interacting particles could be that it is linked by ‘symmerty
principles’ to conservation of the mass of the particles.
Thus we may interpret the parameter y∈S′ in the process
(Xtx,y) by saying that the diffusion with parameters σˉ,bˉ and starting at x is in the state y or
that the diffusion with parameters σˉ,bˉ is in the
state (x,y). This of course, corresponds to the process
(Yt(τxy)) being in the initial state τxy. When we
consider questions such as ergodicity and existence of an invariant
measure, we replace the (initial) deterministic state x by a
random state with a distribution μ. In the context of our
results this raises the question of wether the existence of an
invariant measure and questions of ergodicity can be answered by
randomising both x and y. We refer to [4], Chapter (5),
for some results in this direction.
8 Appendix
We present the proofs of Proposition (5.4) and
Theorem (5.5).
Proof of Proposition (5.4) : Given r>0 we first construct
a pair (Y~r(t,ω,y),η~r(ω,y)) jointly
measurable in (t,ω,y) and (ω,y) respectively such
that for each (t,y),
[TABLE]
where for each y∈Sp, (Ytr(y),ηr,y) is the
solution of equation (19) constructed in Theorem (4.3) with Y≡y. In the construction below we drop the superscript r
until further notice. Recall from Section 2, that {hn,p;n∈Z+d} is the ONB in the Hilbert space Sp.
Since for each y∈Sp,
[TABLE]
where ∣n∣=0∑∞⟨Y(t,ω,y),hn,p⟩p2<∞ for all t≥0 and ω∈Ω, it suffices to show the existence of the map
η~(ω,y) and for each n≥1, a B[0,∞)⊗F∞B⊗B(Sp)/B(R) measurable map Y~n(t,ω,y)
satisfying
[TABLE]
for all (t,ω,y), and satisfying, for each t≥0,y∈Sp, η~(ω,y)=ηy(ω), and Y~n(t,ω,y)=⟨Yt(y),hn,p⟩p almost
surely on the set {t<ηy}. One can then define Y~(t,ω,y) by
[TABLE]
Recall the process (Ytk) satisfying equation (20) (which we now
denote by (Ytk(y)) to make the dependence on y explicit),
constructed for each k≥1,y∈Sp, in the proof of
Theorem (4.3), satisfying for each t≥0,y∈Sp,
[TABLE]
as k→∞; where ηy:=k→∞limηk,y, as in the proof of Theorem (4.3).
It is easy to see that there exists jointly measurable maps
(t,ω,y)→Y~k(t,ω,y) and (ω,y)→η~(ω,y) satisfying, for each t≥0,Y~k(t,ω,y)=Ytk(y) and η~(ω,y)=ηy(ω) almost surely. For the first map we define Y~k(t,ω,y):=τZk(t,ω,y)(y) where the
Rd valued process (Zk(t,ω,y)) is a jointly
measurable version which is indistinguishable for each y from the
process (Ztk(y)) defined in terms of (Yk−1(t,ω,y)) in
the proof of Theorem (4.3). Note that Y0(t,ω,y)≡y.
Thus the joint measurability of Zk(t,ω,y) follows from that
of the stochastic integrals defining Ztk and an induction
argument. Consequently, the map Y~k(t,ω,y) is for each
y, indistinguishable from the process Ytk(y)=τZtk(y)(y). To define the map η~(ω,y), we
first define
[TABLE]
It is easy to check that the map (ω,y)→σ~j(ω,y)
is jointly measurable and satisfies for each y,σ~j(ω,y)=σjy(ω), almost surely ; where we
have explicitly denoted the dependence on y of the stopping time
σj constructed in the proof of Theorem (4.3). The map
η~(ω,y) is now constructed from the map σ~j(ω,y) in the same way as ηy was constructed from
the σj’s and ηk’s in the proof of Theorem (4.3) viz.
η~j(ω,y):=σ~1(ω,y)∧⋯∧σ~j(ω,y) and η~(ω,y):=j→∞limη~j(ω,y).
Fix t≥0,y∈Sp. Since E∥Yt(y)−Yt∧ηyk(y)∥q2→0,q≤p−1,there exists a
subsequence {nk} such that Yt∧ηynk(y)→Yt(y) almost surely in Sq. In particular, for all n=(n1,⋯,nd) and for almost all ω,
[TABLE]
We now construct a set G in the product (t,ω,y)-space using
the subsequence {nk} above as follows. Let G:=n⋂Gn⋂G0 where the intersection is over all
n=(n1,⋯,nd),ni∈Z+d and where the sets
G0,Gn are defined as G0:={(t,ω,y):k→∞lim ∥Y~nk(t∧η~,ω,y)∥q<∞} and Gn:={(t,ω,y):k→∞lim ⟨Y~nk(t∧η~,ω,y),hn,q⟩q exists}. Fix n=(n1,⋯,nd). Define
[TABLE]
Then from the joint measurability of G,η~, and Y~nk(t,ω,y) we get that the map
(t,ω,y)→Yn(t,ω,y) is jointly
measurable. If (t,ω,y)∈G, then
[TABLE]
Since for fixed t≥0,y∈Sp we have Y~nk(t∧η~y,ω,y)=Yt∧ηynk(y) almost surely, it follows from the preceding definitions
that
[TABLE]
almost surely on {t<ηy}. We can now define
[TABLE]
Note that this is not the same as Y~k(t,ω,y) defined earlier in this proof, which were
approximations to Yt(y). Then Y~n(t,ω,y)=⟨Yt(y),hn,p⟩p on {t<ηy} almost surely and
since q≤p−1,
[TABLE]
for every (t,ω,y). Then, as mentioned above, we
construct the map (t,ω,y)→Y~(t,ω,y)
using Y~n(t,ω,y) as its n-th Fourier-Hermite
coefficient, n∈Z+d.
Since the maps Y~,η~ constructed above depend on
r>0, we now make the dependence explicit and patch up the maps
Y~r≡Y~(t,ω,y),η~r≡η~(ω,y) for different r>0. Let rk↑∞. We denote by Y~k(t,ω,y):=Y~rk(t,ω,y),η~k(ω,y):=η~rk(ω,y). Let
[TABLE]
and define
[TABLE]
Then for fixed y, η~(ω,y)=ηy(ω) almost surely follows from the
corresponding equality η~k(ω,y)=ηrk,y(ω), almost surely. Thus, part b) in the statement
of the theorem holds.
For k=1,⋯, define
[TABLE]
and H:=n≥1⋃k≥n⋂Hk. We define
[TABLE]
For fixed (t,y), that Y~(t,ω,y)=Yt(y) almost surely on t<ηy(ω)
follows from the fact that almost surely, Y~k(t,ω,y)=Ytrk(y) on t<ηrk,y(ω). Clearly Y~(t,ω,y) can be extended as a S^p:=Sp⋃{δ} in an obvious manner for t≥η~
to satisfy part a) of the theorem. □
Proof of Theorem (5.5): The proof consists in checking, at
each stage of the construction of measurable maps (t,ω,y)→Y~(t,ω,y) carried out in the previous
theorem, that composition with Y0(ω) at time t yields the
corresponding (approximate) solution with initial value Y0 at
time t.
Recall that for r>0,Y~r(t,ω,y),η~r(ω,y) are the measurable
versions of (Ytr(y),ηr,y) constructed in the previous
proposition. It is sufficient to show that if Y~r(t,ω):=Y~r(t,ω,Y0(ω)),η~r(ω):=η~r(ω,Y0(ω)), then almost surely,
[TABLE]
and that η~r(ω)=ηr(ω)
almost surely. Once this is done for each r>0, we take rk↑∞, define η~k(ω):=η~rk(ω),Y~k(t,ω):=Y~rk(t,ω)
and observe that by pathwise uniqueness of (19), for each t,
almost surely, (t,ω,Y0(ω))∈H,(ω,Y0(ω))∈H0, where the sets H,H0 are as in the previous proposition.
Then Y~(t,ω)=Yt(ω) on {t<η(ω)},
almost surely, follows by pathwise uniqueness.
Recall the approximations (Ytr,k(y),ηr,k,y) k≥1,
for fixed r>0, of the solutions (Ytr(y),ηr,y), of
equation (19) with initial value Y0=y∈Sp in a ball
of radius r around y. It is clear by induction and uniqueness of
the linear equation (20) satisfied by Ytr,k(ω,y), and the
independence of Y0 and (Bt) that for fixed t, (Y~r,k(t,ω,Y0(ω)),η~r,k(ω,Y0(ω))) is the kth approximant to
(Ytr,ηr), the solutions of equation (19) on [0,ηr),
with initial value Y0. Note that ηr(ω)=k↑∞limηr,k,Y0(ω)(ω)=k↑∞limη~r,k(ω,Y0(ω))=η~r(ω,Y0(ω))=:η~r, almost surely, where
the second equality follows from the preceding observation. Thus
from the above observations, we have for each t,
[TABLE]
as k→∞. It
remains to identify the limit as k→∞ of Y~r,k,(t∧ηr,Y0) with Y~r(t,ω,Y0(ω)).
From the above L2 convergence we get the subsequential
convergence
[TABLE]
almost surely. Let G be the set constructed in the
proof of Proposition (5.4), with the above subsequence. Let Yˉr,n(t,ω,y) and Y~r,n(t,ω,y) be as in the
previous proposition, where we have now made the dependence on r
explicit. Then, for fixed t and almost every ω,(t,ω,Y0(ω))∈G, and hence on t<ηr,
[TABLE]
where the last equality follows from the almost sure
subsequential convergence in Sp. Since this is true for
all n=(n1,…,nd), we have
[TABLE]
almost surely on {t<ηr(ω)}.
\hfill□