Disertacija
Disertacija
FACULTY OF SCIENCE
DEPARTMENT OF
MATHEMATICS AND INFORMATICS
Dora Seleši
- doctoral dissertation -
and d are (in the most general setting) assumed to be Colombeau generalized
random processes. Stability and regularity properties of the solutions are also
obtained. Another SDE considered in Chapter 3 is the stochastic Helmholtz
equation, which is solved by means of the Fourier transformation.
I would like to take this opportunity to express my huge thank to my
parents, Erika and György for their love, encouragement, kindness and pa-
tience. Special thanks go to my friends and colleagues for providing constant
support during my work on this dissertation. A complete list of their names
would fill pages. I wish to thank Professor Zagorka Lozanov–Crvenković and
Assistant Professor Danijela Rajter–Ćirić for many fruitful discussions about
Probability and Stochastics. To Professor Marko Nedeljkov I owe a debt of
gratitude for a careful reading of this text and providing a sequence of im-
provements. I would like to mention here his excellent lectures on PDEs;
some of the examples in this dissertation were inspired by things I learned at
the time I was conducting exercises for his lectures. A great thank goes to
Professor Michael Oberguggenberger for his hospitality and support during
my stay at the Insitute of Engineering Mathematics, Geometry and Com-
puter Science in Innsbruck, Austria. In particular, I wish to thank him for
his interesting Seminar on Stochastic Analysis and for our discussions during
my stay in Innsbruck. Last but not least, I want to express my deep ap-
preciation for the constant support, help and encouragement in many areas
given to me by my mentor, Professor Stevan Pilipović. He represents an
inexhaustible source of knowledge, who keeps the doors of his office always
open for his students. I wish to thank him for his comprehensive guidance
through the mathematical universe, showing me some of its most interesting
parts like Functional Analysis, Measure Theory, Generalized Functions (just
to mention some of them), which finally converged to this dissertation.
Contents
Introduction 1
Epilogue 137
References 138
Index 144
Chapter 1
S
inclusion Ep ⊆ E0 ⊆ E−p is continuous. Let E 0 = p∈N0 E−p and endow it
with the inductive topology. Since E is dense in each Ep , the dual of E is
(in algebraic sense) ³E 0 and the´inductive topology on E 0 is equivalent to the
T∞ 0
strong topology on p∈N0 E p .
Nuclear spaces
A countably Hilbert space E is called nuclear , if for every p ∈ N there
exists q ≥ p such that the inclusion mapping Eq → Ep is nulcear. Recall
that the
P inclusion Eq → Ep is nuclear if there exists a CONS {ek } in Eq such
that k∈N |ek |2p < ∞. All finite dimensional spaces are nuclear. Moreover,
nuclear spaces have similar properties as finite dimensional spaces, e.g. if E
is nuclear, then A ⊂ E is compact if and only if it is bounded and closed.
Infinite dimensional Banach spaces are not nuclear. The notion of a nuclear
space is not reduced only to countably Hilbert spaces; one can define also
nuclear topological vector spaces (for this we refer to [Tr]). For example,
the dual of a nuclear countably Hilbert space is a nuclear topological vector
space. Some examples of nuclear spaces we will use in this dissertation are the
Schwartz spaces S, S0 , D, D0 , but also C0 , C0∞ . Also, the Zemanian spaces
A, A0 which are yet to be introduced, are nuclear under some additional
conditions.
Gel’fand triples
Let V be a nuclear space and V 0 its dual. If there exists a Hilbert space H
such that V is dense in H, then the triple V ⊆ H ⊆ V 0 is called a Gel’fand
triple. In fact, here we identified H with its dual H 0 . In case of nuclear
countably Hilbert spaces we have the Gel’fand triple E ⊆ E0 ⊆ E 0 .
Tensor products
In order to define a topological structure on the tensor product of two
topological vector spaces X⊗Y , there are two possibilities at our disposal (see
[Tr]). One possibility is to construct a seminorm topology relying directly on
the seminorm topologies of X and Y (this will be the so called π-topology).
The other possibility is to use an embedding of X ⊗ Y in some space related
to X and Y , over which a topology already exists (this construction will
lead to the so called ε-topology). In general, the π-topology is finer than the
ε-topology.
The π-topology on X ⊗ Y is the strongest locally convex topology such
that the canonical bilinear mapping X × Y → X ⊗ Y , (x, y) 7→ x ⊗ y
10 Generalized Function Spaces
For p ∈ N define |f |−p = |A−p f |L2 (Rd ) and let S−p (Rd ) be the com-
pletition of S(Rd ) with respect to | · |−p . The | · |−p norm topology is
equivalent to the strong topology on the dual space Sp0 (Rd ) defined by
| · |−p = supψ∈S(Rd ),|ψ|p ≤1 |h·, ψi|. Thus, S−p (Rd ) is isomorphic to Sp0 (Rd ).
S
In set–theoretical sense we have S0 (Rd ) = p∈N S−p (Rd ), and moreover
the locally convex inductive limit of the spaces S−p (Rd ) is isomorphic to
12 Generalized Function Spaces
S0 (Rd ) equipped with the strong topology. Thus, S(Rd ) ⊆ L2 (Rd ) ⊆ S0 (Rd )
form a Gel’fand triple. Moreover, we have continuous inclusions
S(Rd ; H) ∼
= S(Rd ) ⊗ H, S0 (Rd ; H) ∼
= (S(Rd ) ⊗ H)0 ∼
= S0 (Rd ) ⊗ H,
where the symbol ⊗ is understood as the tensor product on the test space,
and the π–completition (or in this case the equivalent ²–completition) of
the tensor product of dual spaces in the second case. Also, since H is a
Hilbert space, we may identify it with its dual and thus consider S0 (Rd ; H)
as B(S(Rd ), H) i.e. as the space of bilinear continuous mappings S(Rd )×H →
R. We will denote by h·, ·i this dual bilinear pairing.
Since the family {ei ηk : i, k ∈ N} is an orthogonal basis of L2 (Rd )⊗H, we
get a Hermite basis characterization also in the H–valued P case:
P A function
d
f belongs
P toP S(R ; H) if and only if it is of the form f = i∈N k∈N ai,k ηk ei ,
where i∈N k∈N λpk |ai,k |2 < ∞ for each p ∈ N and the coefficients are given
by ai,k = hf, ηk ei i ∈
PC. PA function f belongsP to S0P
(Rd ; H) if and only if it
is of the form f = i∈N k∈N ai,k ηk ei , where i∈N k∈N λ−p 2
k |ai,k | < ∞ for
some p ∈ N and the coefficients are given by ai,k = hf, ηk ei i ∈ C.
1.3 Zemanian Spaces 13
The set
m
X
2
S = {f ∈ L (I) : f = an ψn , an ∈ C, m ∈ N}, (1.2)
n=1
14 Generalized Function Spaces
i.e. the linear span of the set {ψn : n ∈ N}, is dense in each Ak , k ∈ Z.
Define:
( ∞ ∞
)
\ X X 2k
2 2f
A= Ak = f ∈ L (I) : f = an ψn , ∀k, |an | λn < ∞ ,
k∈N0 n=1 n=1
( ∞ ∞
)
[ X X −2k
A0 = A−k = f= bn ψn : ∃k, fn
|bn |2 λ <∞ .
k∈N0 n=1 n=1
P f 2k
equipped with the norm kϕkp,k = ∞ 2
n=1 |an | (expp λn ) , k ∈ N0 . Thus,
\ [
expp A = expp,k A, expp A0 = expp,−k A.
k∈N0 k∈N0
Laplace operator. By its Hilbert structure, we also may identify W0k,2 (I)
with W −k,2 (I). Thus, we obtain a Gel’fand triple
For further notions and properties of Sobolev spaces we refer to the mono-
graph [Ad].
lim Mp (x)Mp−1
0 (x) = 0.
|x|→∞
We refer the reader to [GS] for the properties of K{Mp } and its strong dual
K0 {Mp }. In this paper we will consider a subclass of such spaces. Namely, as
in [GV, p.82], we will assume that {Mp , p ∈ N0 } are smooth functions such
that
(I) for every k, p ∈ N0 there exist p0 ∈ N0 and C > 0 such that
K−
d = {φ ∈ Kd : supp φ ⊂ (−∞, b) for some b ∈ R},
K0d = Kd ∩ C0 ,
K+
d = {φ ∈ Kd : supp φ ⊂ (b, ∞) for some b ∈ R}.
Let for t ∈ R,
H(t)ts−1
, s>0
Γ(s)
fs (t) = (1.4)
(n)
fs+n (t), s ≤ 0, s + n > 0, n ∈ N,
where H denotes the Heaviside function. Note that for s < 0 we have
(s)
fs = δ, where δ denotes the Dirac delta distribution . For s, r ∈ R we
have fs ∗ fr = fs+r .
For fixed x ∈ R denote
fsx− (t) = fs (x − t), t ∈ R,
fsx+ (t) = fs (t − x), t ∈ R.
18 Generalized Function Spaces
Clearly, fsx− , fsx+ ∈ C s−2 (R) and supp fsx− ⊂ (−∞, x], supp fsx+ ⊂ [x, ∞).
Note that for φ ∈ K{Mp } following holds:
µZ ¶(d+2) µZ x ¶(d+2)
x−
fd+2 (t)φ(t)dt = fd+2 (x − t)φ(t)dt = φ(x),
R −∞
µZ ¶(d+2) µZ ∞ ¶(d+2)
x+
fd+2 (t)φ(t)dt = fd+2 (t − x)φ(t)dt = φ(x).
R x
Now, let g ∈ K0 {Mp } such that supp g ⊂ (a, ∞) for some a ∈ R. There
exist C > 0 and d ∈ N such that
|hg, φi| ≤ Ckφkd , φ ∈ K{Mp }.
Then, g can be extended onto every φ ∈ Kd . Moreover, g can be extended
onto every ϕ ∈ C ∞ (R) such that supp ϕ ⊂ (−∞, b), for some b ∈ R. Thus, g
x− x−
can be extended to every fd+2 , x ∈ R. Note that the mapping x 7→ hg, fd+2 i
x−
is continuous. Thus, for G(x) = hg, fd+2 i, x ∈ R, we have
G(d+2) = g, (1.5)
where the derivative is understood in distributional sense. A similar ar-
x+
gument holds for G(x) = hg, fd+2 i, x ∈ R, where g ∈ K0 {Mp } such that
supp g ⊂ (−∞, a) for some a ∈ R.
processes (we will put some space of generalized random variables in place
of H).
holds. Here hω, φi is the dual pairing of S0 (Rd ) and S(Rd ). The triplet
(S0 (Rd ), B, µ) is called the white noise space and µ is called the white noise
measure or the Gaussian measure on S0 (Rd ).
From now on we assume that the basic probability space (Ω, F, P ) is
the space (S0 (Rd ), B, µ). Put (L)2 = L2 (S0 (Rd ), B, µ). It is a Hilbert space
equipped with the inner product (F |G)(L)2 = Eµ (F G).
QmIn a multi-dimensional
d d
case, for a given m ∈ N, m > 1, define
d 0
Qm 0Sm d=
i=1 iS (R ), where S i (R ) is a copy of S(R ), and let S m = i=1 Si (R ).
0
Equip the space Sm with the product Borel σ-algebra and with the product
measure µm = µ×· · ·×µ. The triple (S0m , B, µm ) is called the m-dimensional
d-parameter white noise L space. Put (L)2,m = L2 (S0m , µm ), and for N ∈ N,
N
N > 1, let (L)2,m,N = k=1 (L)2,m be the direct sum of N identical copies of
the m-dimensional d-parameter white noise space. Here the finite dimension
N is the state space dimension. Later, instead of the N –dimensional state
space, we will also consider the infinite–dimensional case, when the state
space is a separable Hilbert space H.
1.7 Spaces of Generalized Random Variables 21
exp(S)m,N
ρ ⊆ (S)m,N
ρ ⊆ (L)2,m,N ⊆ (S)m,N m,N
−ρ ⊆ exp(S)−ρ ,
∞ X
X ∞ X
X
f (ω) = ai,α Hα(m) (ω)ei , ai,α = hf, Hα(m) ei i ∈ R, α!|ai,α |2 < ∞,
i=1 α∈I i=1 α∈I
X X
f (ω) = aα Hα(m) (ω), aα = hf, Hα(m) i(L)2,m ∈ H, α!kaα k2H < ∞.
α∈I α∈I
24 Generalized Function Spaces
X ∞
X
kf k2ρ,p;H = α!1+ρ kaα k2H (2N)pα = kai k2ρ,p .
α∈I i=1
X ∞
X
kF k2−ρ,−q;H = α! 1−ρ
kbα k2H (2N)−qα = kbi k2−ρ,−q .
α∈I i=1
The same construction can be carried out for the exponential growth rate
spaces. Let exp S(H)m ρ be the space of functions f ∈ L
2,m
(Ω; H), f (ω) =
P∞ P (m)
i=1 α∈I ai,α Hα (ω)ei , ai,α ∈ R, such that for all p ∈ N0 ,
∞ X
X ∞
XX
α α
kf k2ρ,p,exp;H = α!1+ρ |ai,α |2 ep(2N) = α!1+ρ |ai,α |2 ep(2N) < ∞.
i=1 α∈I α∈I i=1
Both S(H)m m
ρ and exp S(H)ρ are countably Hilbert spaces and
exp S(H)m m
ρ ⊆ S(H)ρ ⊆ L
2,m
(Ω; H) ⊆ S(H)m m
−ρ ⊆ exp S(H)−ρ .
In general, S(H)m m
ρ and exp S(H)ρ are not nuclear spaces. They would
be nuclear e.g. if H were finite–dimensional (see [Tr, Prop.50.1,Cor.]).
Note that since (S)m m
ρ and exp(S)ρ are nuclear spaces, by [Tr, Prop.50.7]
we have again (as for the Schwartz spaces in the deterministic case) an iso-
morphism with tensor product spaces:
S(H)m ∼ = (S)m ⊗ H,
−ρ −ρ exp S(H)m ∼ = exp(S)m ⊗ H. −ρ −ρ
Chapter 2
Generalized Stochastic
Processes
D. Rajter and others (see [OR1], [OR2], [RO], [AHR1], [AHR2], [NR].) In
this framework a GRP is defined as a mapping X : Ω → G(T ), (G(T ) is the
simplified Colombeau algebra of generalized functions) having the property
that there exists a net of functions R² : Ω × T → C, ² ∈ (0, 1), such that
R² (ω, ·) represents X(ω) for a.e. ω ∈ Ω, and for every ² ∈ (0, 1), (ω, t) →
R² (ω, t) is a measurable (classical) stochastic process. Recall that D0 (T ) can
be embedded into G(T ), thus this is a generalization of the Walsh definition
of a GRP.
Also, in [LPP], Colombeau GRPs of second order were considered as
elements of the vector–valued Colombeau algebra G(T ; L2 (Ω)), where L2 (Ω)
is the Hilbert space of random variables with finite second moments.
We will follow a similar idea in this dissertation and replace L2 (Ω) with
the more general (S)−1 and consider the Colombeau extension of a GRP
(I) i.e. we will consider GRPs as elements of G(Rd ; (S)−1 ). Similarly as in
the deterministic theory of distributions, where S0 (Rd ) can be embedded into
G(Rd ), we will have the vector–valued analogy of it: The space of GRPs (I)
L(S(Rd ), (S)−1 ) can be embedded into G(Rd ; (S)−1 ).
(iii) Φ can be represented in the form (2.4) and there exists k1 ∈ N0 such
that ∞
X j
kfj k2−k (2N)−k1 α < ∞. (2.6)
j=1
32 Generalized Stochastic Processes
such that ai (t) ∈ L1loc (R), i ∈ N. Then due to (2.1), there is a GRP (I),
denoted by Ũ associated with U , such that
Z X
∞ ∞
X
[Ũ , ϕ](ω) = ai (t)Hαi (ω)ϕ(t)dt = hãi , ϕiHαi (ω), ω ∈ S0 (R),
R i=1 i=1
where ãi ∈ S0 (R) is the generalized function associated with the function
ai (t) ∈ L1loc (R), i ∈ N. Thus, Ũ has expansion
∞
X
Ũ = ãi ⊗ Hαi .
i=1
This means that the expansion theorems for GRPs (I) are consistent with
the expansion theorems for GRPs (O).
exp
Expansion theorems for A∗ = L(A, exp(S)−1 )
Analogously to the previous theorems concerning GRPs (I) with values in
(S)−1 , one can consider GRPs (I) taking values in other spaces of generalized
stochastic functions. The expansion theorems can be carried over, mutatis
mutandis, to these GRPs.
Consider for example the space exp(S)−1 , which will provide a larger class
of GRP (I). Let exp A∗ = L(A, exp(S)−1 ) and exp A∗k = L(Ak , exp(S)−1 ) be
GRPs (I) and GRPs (I) of R-order k, respectively. Further, let all the other
terms be defined analogously as for A∗ ; i.e. we replace (S)−1 with exp(S)−1
in Definition 2.1.1 and else where necessary.
(iii) Φ can be represented in the form (2.7) and there exists k1 ≥ 0, such
that ∞
X αj
kfj k2−k e−k1 (2N) < ∞.
j=1
For examples of GRPs (I) and definition of the Wick product in spaces
A∗ refer to [PS1].
Lemma 2.2.1 Let hfn in∈N be a sequence in A0 and hθαj ij∈N be a sequence
in S(H)−1 . If there exists k0 ∈ N0 such that for any bounded set B ∈ A,
∞
X
sup |hfj , ϕi| · kθαj k−1,−k0 ;H < ∞, (2.11)
ϕ∈B
j=1
P∞
then j=1 fj ⊗ θαj defined by
∞
X m
X
fj ⊗ θαj = lim fj ⊗ θαj
m→∞
j=1 j=1
is an element of A(H)∗ .
P Pm
Proof. Denote Υ = ∞ j=1 fj ⊗ θαj and Υm = j=1 fj ⊗ θαj , m ∈ N. Clearly,
∗
Υm ∈ A . If (2.11) holds, then Υm ∈ L(A, S(H)−1,−k0 ).
The sequence of partial sums Υm , m ∈ N, is a Cauchy sequence in A(H)∗
because, for given ε > 0 and m > n,
m
X
k[Υm , ϕ] − [Υn , ϕ]k−1,−k0 ;H = |hfj , ϕi| · kθαj k−1,−k0 ;H < ε,
j=n+1
2.2 Hilbert Space Valued GRPs (I) 35
it follows that
à m
!
X
sup kΥm k∗ ≤ sup sup{ |hfj , ϕi| · kθαj k−1,−k0 ;H : ϕ ∈ A, kϕk ≤ 1}
m∈N m∈N
j=1
∞
X
≤ sup{ |hfj , ϕi| · kθαj k−1,−k0 ;H : ϕ ∈ A, kϕk ≤ 1} < ∞,
j=1
if we choose n, m large enough. This yields that hΦm im∈N is a Cauchy se-
quence in A(H)∗k . Also, (2.13) implies that
X∞ X
∞
j
sup kΦm k∗2
−k;H ≤ sup{ |hfij , ϕi|2 (2N)−k0 α : ϕ ∈ Ak , kϕkk ≤ 1} < ∞.
m∈N
i=1 j=1
then Φ ∈ A(H)∗k .
Proof. According to the Cauchy-Schwartz inequality, it is obvious that
(2.15) implies (2.13). ¤
such that aij (t) ∈ L1loc (R), i, j ∈ N. Then there is a H–valued GRP (I),
denoted by Ũ associated with U , such that
Z X∞ X∞ ∞ X
X ∞
[Ũ , ϕ](ω) = aij (t)Hαj (ω)ei ϕ(t)dt = hãij , ϕiHαj (ω)ei , ω ∈ S0 (R),
R i=1 j=1 i=1 j=1
38 Generalized Stochastic Processes
where ãij ∈ S0 (R) is the generalized function associated with the function
aij (t) ∈ L1loc (R), i, j ∈ N. Thus, Ũ has expansion
∞ X
X ∞
Ũ = ãij ⊗ Hαj ei .
i=1 j=1
The expansion theorems for exp A(H)∗ = L(A, exp S(H)−1 ) can also be
stated as in the case of a one–dimensional state space:
Corollary 2.2.2 If Φ can be represented in the form (2.16) and there exists
k1 ≥ 0, such that
X∞ X ∞
αj
kfij k2−k e−k1 (2N) < ∞,
i=1 j=1
exp
then Φ ∈ A(H)∗k .
if we choose p > 1.
Note also that κ̃k is an element of S0 (R; H).
5
(ii) Let l > 12 and t1 , t2 , t3 , . . . ∈ R such that t1 ≤ t2 ≤ t3 ≤ · · · → ∞.
It is known that the Dirac delta distributions δtj , j ∈ N, belong to
A−l = S−l (R). Let
(to avoid confusion: the first delta is the Kronecker symbol, the second
one is the Dirac distribution). With
∞
X
∆k ⊗ Hαk (2.18)
k=1
is given a GRP (I) which is not a GRP (O). We will check condition
P 1
(2.15). Since δtj (x) = ∞n=1 ξn (tj )ξn (x), j ∈ N, and ξn = O(n
− 12
),
n ∈ N, we have
∞ X
X ∞ ∞
X
j k
δn(i,j),k kδtj k2−l (2N)−pα = kδtk k2−l (2N)−pα
i=1 j=1 k=1
X∞ X∞
k
= |ξn (tk )|2 (2n)−2l (2N)−pα
k=1 n=1
X∞ ∞
X 1
−pαk
≤C (2N) n− 6 (2n)−2l < ∞,
k=1 n=1
for some constant C > 0 and p > 1. Hence, the process given by (2.18)
meets the definition of a GRP (I).
40 Generalized Stochastic Processes
Corollary 2.2.3 If Φ can be represented in the form (2.19) and there exists
k1 ∈ N0 such that
X∞
j
kfj k2−k;H (2N)−k1 α < ∞, (2.21)
j=1
then Φ ∈ A(H)∗k .
∂
where f
∂t ij
is the distributional derivative of fij in A0 .
t ∈ I, ω ∈ S0 (R),
P∞respectively, where 0
P∞ uijk, gij ∈ A , i, j ∈ N.
k
Let uij = k=1 aij ψk , gij = k=1 bij ψk , i, j ∈ N. Then we have
∞ X
X ∞
P (R)u = P (R)uij ⊗ Hαj ei
i=1 j=1
∞ X
∞
à ∞
! ∞ X
∞ X
∞
X X X
= P (R) akij ψk ⊗ Hαj ei = akij P (λek )ψk ⊗ Hαj ei .
i=1 j=1 k=1 i=1 j=1 k=1
(2.24)
In order to solve (2.23) we may use the method of undetermined coefficients.
From (2.23) and (2.24) we have
X∞
akij P (λek )ψk = gij , i, j ∈ N.
k=1
Finally, we obtain the system
akij P (λek ) = bkij , i, j, k ∈ N.
b k
First case: If P (λek ) 6= 0 for all k ∈ N, then akij = P (λijf ) , i, j, k ∈ N, and the
k
solution of the equation is given by
∞ X ∞
̰ !
X X bkij
u= ψ ⊗ Hαj ei .
P (λek ) k
i=1 j=1 k=1
In this case the solution exists and it is unique. Note, there exists a constant
C > 0 such that P (λek ) ≥ C, for all k ∈ N. Thus,
°∞
∞ °X
°2 °∞ °2
X∞ X
bk ° 1 X∞ X∞ °X °
° ij ° −pαj ° ° j
° ψk ° (2N) ≤ 2 ° k
bij ψk ° (2N)−pα < ∞,
° e
P (λk ) ° C ° °
i=1 j=1 k=1 −r i=1 j=1 k=1 −r
First we recall the definition and some basic properties of the Wick prod-
uct in the Kondratiev spaces (see [HØUZ]).
Definition
P 2.3.1 Let F, G ∈ P (S)−1 be given by their chaos expansion
F (ω) = α∈I fα Hα (ω), G(ω) = β∈I gβ Hβ (ω), fα , gβ ∈ R. The Wick prod-
uct of F and G is the unique element in (S)−1 defined by:
à !
X X
F ♦G(ω) = fα gβ Hγ (ω).
γ∈I α+β=γ
Definition
P∞2.3.2
P Let F, G ∈ S(H)−1 bePgiven by their chaos expansion
∞ P
F (ω) = i=1 α∈I fi,α Hα (ω)ei , G(ω) = i=1 β∈I gi,β Hβ (ω)ei , fi,α , gi,β ∈
R. The Wick product of F and G is the unique element in S(H)−1 defined
2.3 The Wick Product of GRPs 45
by:
∞ X
à !
X X
F ♦G(ω) = fi,α gi,β Hγ (ω)ei
i=1 γ∈I α+β=γ
X∞
= (Fi ♦Gi (ω))ei ,
i=1
P P
where Fi (ω) = α∈I fi,α Hα (ω), Gi (ω) = β∈I gi,β Hβ (ω) ∈ (S)−1 .
This definition is legal, since S(H)1 and S(H)−1 are closed under Wick
multiplication. In the same manner we can define F ♦G for F, G ∈
exp S(H)−1 and it is an easy exercise to show (combining the methods in
[MFA] and [PS1]) that F ♦G ∈ exp S(H)−1 . Also, if F, G ∈ exp S(H)1 then
F ♦G ∈ exp S(H)1 .
Note that one can also Wick–multiply a R–valued and a H–valued gen-
eralized P
random variable: If F P ∈ (S)P −1 and G ∈ S(H)−1 have the forms
F (ω) = α∈I fα Hα (ω), G(ω) = ∞ g H (ω), then F ♦G ∈ S(H)−1
P ∞ P ³P
i=1 ´ i,β β
β∈I
Definition
P∞ 2.3.3 LetPf, g ∈ A0 be generalized functions given by expansions
f = k=1 ak ψk , g = ∞ k=1 bk ψk . Define f ¦ g to be the generalized function
0
from A given by
∞
X X
f ¦g = ai bj ψn . (2.25)
n=1 i,j∈N
i+j=n+1
46 Generalized Stochastic Processes
P P∞
Proposition 2.3.1 If f = ∞ i=1 ai ψi ∈ A−k and g = i=1 bi ψi ∈ A−l , then
f ¦ g ∈ A−(k+l+p) . Moreover,
λg ee
i+j ≤ C λi λj , i, j ∈ N.
P P∞
If f = ∞i=1 ai ψi ∈ Ak and g = i=1 bi ψi ∈ Ak , then f ¦ g ∈ A(k−p) .
Example 2.3.1 Let I = (−π, π), R = −iD. Then, the orthonormal basis of
einx
A is given by the family of trigonometric functions ψn (x) = √ 2π
, and λn = n,
n = 0, ±1, ±2, ±3, . . . etc. In this case, we leave the enumeration running
through the set of integer numbers Z. Instead doing a reordering in the space
A, we make a slight modification in³ the definition ´ of the Wick product and
P∞ P
put instead of (2.25) f ¦ g = n=1 i,j∈Z ai bj ψn . Then,
i+j=n
1
ψi ¦ ψj = ψi+j , ψi · ψj = √ ψi+j i, j ∈ Z,
2π
i.e. the Wick product and the ordinary product coincide up to a constant.
Theorem 2.3.1 The Wick product Φ¨Ψ from the previous definition is a
GRP (I), precisely Φ¨Ψ ∈ A∗k+l+p .
The Wick product for GRPs (I) taking values in exp(S)−1 can be defined
in an analogue way.
Theorem 2.3.3 The Wick product Φ¨Ψ from the previous definition is an
H–valued GRP (I), precisely Φ¨Ψ ∈ A(H)∗k+l+p .
48 Generalized Stochastic Processes
Proof. Due toPLemma 2.3.1 it follows that fi,s ¦gi,r ∈ A−k−l−p for all i, s, r ∈
N,
P∞and Pthus, αs +αr =αn+1 fi,s ¦ gi,r ∈ A−k−l−p , n ∈ N. Since [Φ¨Ψ, ϕ] =
∞ P
i=1 n=1 αs +αr =αn+1 hfi,s ¦ gi,r , ϕiHαn ei , it remains to prove
¯ ¯2
∞ X
X ∞ ¯ X ¯
¯ ¯ n
sup ¯ hfi,s ¦ gi,r , ϕi¯ (2N)−r3 α < ∞,
ϕ∈B n=1
¯ ¯
i=1 αs +αr =α n+1
¯ ¯2
∞ X
X ∞ ¯ X ¯
¯ √ ¯ n
≤ sup ¯ kfi,s k−k kgi,r k−l kϕkk+l+p M ¯ (2N)−r3 α ≤
ϕ∈B n=1
¯ ¯
i=1 αs +αr =αn+1
¯ ¯2
∞
X ∞ ¯
X X ¯
¯ r αs r αr ¯
−qαn − 12 − 22
M (2N) ¯ kfi,s k−k (2N) kgi,r k−l (2N) sup kϕkk+l+p ¯
n=1
¯ ϕ∈B ¯
i=1 αs +αr =αn+1
∞ ∞
Ã∞ ∞
!
X X X X
−qαn 2 −r1 αs 2 −r2 αr 2
≤M (2N) kfi,s k−k (2N) kgi,r k−l (2N) sup kϕkk+l+p
ϕ∈B
n=1 i=1 s=1 r=1
∞
X ∞ X
X ∞ ∞ X
X ∞
n s r
=M (2N)−qα kfi,s k2−k (2N)−r1 α kgi,r k2−l (2N)−r2 α sup kϕk2k+l+p
n=1 ϕ∈B
i=1 s=1 i=1 r=1
s +αr s r
< ∞, where we used the property (2N)α = (2N)α (2N)α . ¤
Note that the Wick product ♦ acting on S(H)−1 can be embedded into
¨ acting on A(H)∗ .
Also, following theorem holds, similarly as in finite dimensional case.
i=1
2
j=1 kfi,j k−k e
−r1 (2N)α
< ∞ and ∞ i=1
2 −r2 (2N)α
j=1 kgi,j k−l e < ∞.
2.3 The Wick Product of GRPs 49
where ½
ψj (t), n = m(i, j)
wi,n =
0, else
and m(i, j) is defined as in Example 1.7.3.
Let g ∈ A(H)∗r be of the form
∞ X
X ∞
g(t, ω) = gi,n (t) ⊗ Hεn (ω)ei , t ∈ I, ω ∈ Ω,
i=1 n=1
such that gi,n ∈ A−r , i, n ∈ N; and assume there exist q ≥ 0, such that
∞ X
X ∞ ∞ X
X ∞
kgi,n k2−r (2N)−qεn = kgi,n k2−r (2n)−q < ∞. (2.30)
i=1 n=1 i=1 n=1
which yields
∞
X −2r
|gi,n,k |2 λek < ∞. (2.31)
k=1
Then,
∞ X
∞
̰ !
X X
P (R)X(t, ω) = ai,n,k P (λek )ψk (t) ⊗ Hεn (ω)ei , t ∈ I, ω ∈ Ω,
i=1 n=1 k=1
(2.32)
and due to the definition of Wick product:
Since
∞
X ∞
X
ai,s (t) ¦ wi,r (t) = ai,s,k ψk (t) ¦ wi,r (t) = ai,s,k ψk (t) ¦ ψj (t)
k=1 k=1
∞
X
= ai,s,k ψk+j−1 (t),
k=1
X ∞
X X∞ ∞
X
ai,s,k ψk+j−1 (t)+ gi,n,k ψk (t) = ai,n,k P (λek )ψk (t), i, n ∈ N.
s,r∈N k=1 k=1 k=1
s+r=n+1,r=m(i,j)
(2.35)
2.3 The Wick Product of GRPs 51
From the system of equations (2.35) one can recursively determine the
coefficients ai,n,k , i, n, k ∈ N.
For i = 1, n = 1 we have only one possibility how to get r + s = 2
(s = 1, r = 1) and for r = 1 the corresponding j is j = 1. Thus, (2.35) gives
∞
X ∞
X ∞
X
a1,1,k ψk (t) + g1,1,k ψk (t) = a1,1,k P (λek )ψk (t), t ∈ I,
k=1 k=1 k=1
which implies
g1,1,k
a1,1,k = , k = 1, 2, . . . (2.36)
P (λek ) − 1
∞
X ∞
X ∞
X
a1,1,k ψk (t) + g1,2,k ψk (t) = a1,2,k P (λek )ψk (t), t ∈ I,
k=1 k=1 k=1
and since a1,1,k are known from the previous step, now we get
a1,1,k + g1,2,k
a1,2,k = , k = 1, 2, . . .
P (λek )
∞
X ∞
X ∞
X ∞
X ∞
X
a1,1,k ψk+1 + a1,3,k ψk + a1,2,k ψk + g1,3,k ψk = a1,3,k P (λek )ψk .
k=1 k=1 k=1 k=1 k=1
a1,2,1 + g1,3,1
a1,3,1 = ,
P (λe1 ) − 1
a1,1,k−1 + a1,2,k + g1,3,k
a1,3,k = , k = 2, 3, . . .
P (λek ) − 1
and thus,
kai,n k2−(r+p) < ∞, i, n ∈ N,
which yields that all coefficients ai,n , i, n ∈ N, belong to A−(r+p) .
Also, it holds that,
∞ X
X ∞ ∞ X
X ∞
εn 1 −2sn
kai,n k2−(r+p) e−s(2N) 2
≤D M (2n)2q Q2n ( )e .. (2.40)
i=1 n=1 i=1 n=1
K
X
Cα2 (2N)−pα < ∞. (2.42)
α∈I
(ii) Let us check first that N(A; (S)−1 ) is a subalgebra of EM (A; (S)−1 ).
Let F² , G² ∈ N(A; (S)−1 ) be of the form as in (i). Let k ∈ N0
56 Generalized Stochastic Processes
X X 1
≤ Cα Dβ ²b+nα −bβ ≤ Cα Dβ ²bβ +nα −bβ ≤ C −k M − 2 Nγ ²nγ , ² < ²0 .
α+β=γ α+β=γ
P
It is clear that γ∈I Nγ2 (2N)−rγ < ∞ for r = p + q + 2. Thus, G² ♦F² ∈
N(A; (S)−1 ). ¤
The same holds true if A is generated by the operator R = −iD (see Example
2.3.1). In this case the Wick product and the ordinary product coincide up
to a constant.
• Let N(I; (S)−1 ) denote the vector space of functions R² ∈ EM (I; (S)−1 )
with the property that for each compact subset Uk ⊂⊂ I, k ∈ N, there
exist a sequence {Cα }α∈I of positive numbers, ²0 ∈ (0, 1), p ∈ N0 ,
and for all sequences {aα }α∈I bounded from below, and for all |β| ≤ k
following hold:
• Define the multiplication in EMP (I; (S)−1 ) and N(I; (S)−1 ) in the
following
P way: For F² (x, ω) = α∈I fα,² (x) ⊗ Hα (ω), G² (x, ω) =
α∈I gα,² (x) ⊗ Hα (ω) let
à !
X X
F² ♦G² (x, ω) = fα,² (x) · gα,² (x) ⊗ Hγ (ω), (2.51)
γ∈I α+β=γ
i.e. we use the Wick product for multiplication in (S)−1 and the ordi-
nary product in C ∞ (I). Clearly, (C ∞ (I), ·) is an algebra.
• Then N(I; (S)−1 ) is an ideal of EM (I; (S)−1 ) and the quotient space
d
y(ω, t) = f (ω, t), y(ω, 0) = X(ω), (2.52)
dt
where the value at zero, y(ω, t) at t = 0, is understood in the sense of
Lojasiewicz.
To solve (2.52) we use the decomposition of K{Mp }: Every φ ∈ K{Mp }
is of the form
Z
φ = Cφ φ0 + ψ, where φ0 , ψ ∈ K{Mp }, Cφ = φ(t)dt,
R
Z
d
φ0 = 1, and ψ = θ, for some θ ∈ K{Mp }.
R dt
Now, fixing ω ∈ Ω, and using the continuity of f (ω, ·) : K{Mp } → R, we
solve this equation in the usual way:
∞
X ∞
X k0
X
|Xi (ω)|2 |Xi (ω)|2
≤ αi2 kφi k2V · = kvk2V .
i=1 i=1
kφi k2V i=1
kφi k2V
nicely to x0 ∈ Rn , if there exists α > 0 such that each En lies in an open ball
B(x0 , rn ) and m(En ) ≥ αm(B(x0 , rn )), n ∈ N, rn → 0.
Put
κ(En )(x)
δn (x − x0 ) = , n ∈ N, x ∈ Rn . (2.55)
m(En )
The proof of the following results can be found in [LP2] and [LPP].
Q
Theorem 2.5.1 Let G = ni=1 (αi , βi ) ⊂ Rn , −∞ ≤ αi < βi ≤ ∞, i =
1, 2, ..., n, and let ξ be a GRP on Ω × Lr (G), r > 1.
Cξ (ϕ, ψ)
d Z Z
d X
X
= E(fα (·, t)fβ (·, s))Md (t)Md (s)ϕ(α) (t)ψ (β) (s)dtds.
α=0 β=0 R R
Proof. a) Note that in [LP] and [LP1] more restricted problems were con-
sidered (see also [Ha] and [Ul]).
Since for every ω ∈ Ω, ξ(ω, ·) is in K0 {Mp }, it follows that for every ω ∈ Ω
there exist C(ω) > 0 and p(ω) ∈ N such that
We can assume that p(ω) ≥ C(ω). For every ϕ ∈ K{Mp } and N ∈ N, put
\
AN (ϕ) = {ω ∈ Ω : |ξ(ω, ϕ)| < N kϕkN,2 }, AN = AN (ϕ).
ϕ∈K{Mp }
Since K{M
T p } is separable, it contains a countable dense subset D and
AN = ϕ∈D AN (ϕ) ∈ F. Thus, from
∞
[
Ω= AN and AN ⊂ AN +1 , N ∈ N,
N =1
it follows that for given ε > 0 there exists an integer d such that P (Ad ) ≥
1 − ε. Denote M = Ad . It follows
then ξ1 (ω, ϕν ) → 0, ν → ∞.
Qd
Let Γd = i=0 L2 (R) and endow it with the scalar product ((ϕα ), (ψα )) =
Pd R
α=0 R ϕα ψα dt, (ϕα ), (ψα ) ∈ Γd . Clearly, Γd is a Hilbert space. Define a
mapping θ : K{Mp } → Γd by θ(ϕ) = (Md ϕ, Md ϕ0 , . . . , Md ϕ(d) ), ϕ ∈ K{Mp },
which is injective, and denote ∆ = θ(K{Mp }). Note that
d
∆ = θ(HM ). (2.62)
66 Generalized Stochastic Processes
The existence of the limit follows from (2.61) and (2.62). Thus,
¯ ψ⊥ ∈ ∆
F (ω, ψ̃) = F (ω, ψ), ψ̃ ∈ Γd , ψ̃ = ψ + ψ ⊥ , ψ ∈ ∆, ¯ ⊥.
It is endowed with the natural norm such that it is isometric to L2 (R), for
every α ∈ {0, 1, ..., d}. Let ψ ∈ Γd . Denote by [ψ]α the corresponding element
in Γd,α . The βth coordinates of [ψ]α are equal to zero for β 6= α and the αth
coordinate is equal to ψα . Since F is a GRP, it follows that Fα = F |Γd,α is a
GRP on Ω × Γd,α i.e. on Ω × L2 (R), for every α, α ≤ d.
By Theorem 2.5.1, for every α = 0, 1, . . . , d, there exists a function fα :
Ω × R → C such that fα (·, t) is measurable for every t ∈ R, fα (ω, ·) ∈ L2 (R),
ω ∈ Ω, and
Z
Fα (ω, ϕ) = fα (ω, t)ϕ(t)dt, ϕ ∈ L2 (R), ω ∈ Ω.
R
and
d
X
kF (ω, ·)k0Γd = kfα (ω, ·)kL2 (R) ≤ S(ω) ≤ d, ω ∈ Ω,
α=0
2.5 GRPs of Type (II) 67
where k · k0Γd is the dual norm. Now, the assertion follows by (2.59).
The proof of the last assertion in a) follows by repeating the previous
proof starting from relation (2.60). It follows that ξ is of the form (2.56) for
every ω ∈ Ω.
b) Obviously, Cξ (ϕ, ψ) = E(ξ(·, ϕ)ξ(·, ψ)), ϕ, ψ ∈ K{Mp } is bilinear. The
continuity follows from
Xd Z Xd Z
(α)
= E(( fα (·, t)Md (t)ϕ (t)dt)( fα (·, s)Md (s)ψ (α) (s)ds))
α=0 R α=0 R
d X
X d Z Z
= E( fα (·, t)Md (t)ϕ(α) (t)fβ (·, s)Md (s)ψ (β) (s)dtds)
α=0 β=0 R R
d Z Z
d X
X
= ( E(fα (·, t)fβ (·, s))Md (t)Md (s)ϕ(α) (t)ψ (β) (s)dtds.
α=0 β=0 R R
supp Ξ1 (ω, ·) ⊂ [−∞, −a], supp Ξ2 (ω, ·) ⊂ [−2a, 2a], supp Ξ3 (ω, ·) ⊂ [a, ∞];
|Ξi (ω, t)|
sup < ∞, i = 1, 2, 3;
t∈R Md˜i (t)
ξ(ω, ·) = Ξ1 (ω, ·)(s1 ) + Ξ2 (ω, ·)(s2 ) + Ξ3 (ω, ·)(s3 ) , ω ∈ M. (2.64)
b) In particular, if there exist D(ω) > 0, ω ∈ Ω and d ∈ N0 such that
|ξ(ω, ϕ)| ≤ D(ω)kϕkd,2 , ω ∈ Ω, ϕ ∈ K{Mp }, (2.65)
then representation (2.64) is valid on the whole Ω.
c) Let ξ be as in part a) and ε ∈ (0, 1). Let C(ω), ω ∈ Ω, be a Gaussian
random variable. The GRP ξ has the point value C(ω) at t = t0 ,
denoted by ξ(ω, t0 ) = C(ω), ω ∈ Ω, if and only if there exist M ∈ F
such that P (M ) ≥ 1 − ε, s3 , d3 ∈ N, a ∈ R and a Gaussian process
Ξ3 (ω, t), (ω, t) ∈ Ω × R with following properties:
supp Ξ3 (ω, ·) ⊂ [t0 − 2a, ∞),
ξ(ω, ·) = Ξ3 (ω, ·)(s3 ) , ω ∈ M, in a neighborhood of t0 ,
|Ξ3 (ω, t)|
sup < ∞,
t∈R Md3 (t)
and
Ξ3 (ω, t) C(ω)
lim s
= , ω ∈ M. (2.66)
t→t0 (t − t0 ) 3 s3 !
In particular, if there exist D ∈ Z 2 , D(ω) > 0, ω ∈ Ω and d ∈ N0 such
that
|ξ(ω, ϕ)| ≤ D(ω)kϕ ∗ fd3 + 1 kd,2 , ω ∈ Ω, ϕ ∈ K{Mp }, (2.67)
2
In the sequel we consider ξ˜3 . Let fs be given as in (1.4) and let κ ∈ C ∞ (R),
supp κ ⊂ [−a, ∞), κ ≥ 0, κ(t) = 1, t ≥ 0.
We have that for s3 = d3 + 2, and every x ∈ R,
t 7→ κ(t)fsx−
3
(t) = κ(t)fs3 (x − t) ∈ K−
d3
κfsx−
3
→ κfsx30 − in K−
d3 , as x → x0 .
For every x ∈ R,
ω 7→ ξ˜3 (ω, κfsx−
3
)
70 Generalized Stochastic Processes
Z x
kκfsx−
3
k2d3 ,2 = sup |Md3 (t)[κ(t)fs3 (x − t)](i) |2 dt
i≤d3 −a
Z ¯ · ¸ ¯2
¯
x+a
κ(x − u)u d3 +1 (i) ¯
¯ ¯
= sup ¯Md3 (x − u) ¯ du
i≤d3 0 ¯ Γ(d3 + 2) ¯
Z x+a ¯¯ Xi µ ¶· d3 +1
¸(n) ¯¯2
¯ 1 i u ¯
≤ C1 · Md2˜3 (x) sup ¯ ¯ du
i≤d3 0 ¯ Md3 (u) n Γ(d3 + 2) ¯
n=0
we have
v
u
u Z ¯¯ · d3 +1
¸(i) ¯¯2
|Ξ3 (ω, x)| ¯ 1 u ¯
sup ≤ C2 tsup ¯ ¯ du < ∞,
x∈R Md˜3 (x) i≤d3 R ¯ Md3 (u) Γ(d3 + 2) ¯
is equivalent to
Ξ3 (ω, x) C(ω)
lim s
= , ω ∈ M. (2.71)
x→0 x 3 s3 !
Assume that (2.70) holds. Then, for φ(x) = fsx−
3
, x ∈ R, we obtain
Ξ3 (ω, x) 1 t
lim s
= lim hξ˜3 (ω, t), (1 − )s3 −1 i
x
x→0 3 x→0 xΓ(s3 ) x
Z 1
˜ 1 s3 −1 1 C(ω)
= lim hξ3 (ω, tx), (1 − t) i = C(ω) (1 − t)s3 −1 dt = .
x→0 Γ(s3 ) 0 Γ(s3 ) s3 !
Conversely, assume (2.71) holds. Then, for φ ∈ K{Mp } we have
1 x
limhξ˜3 (ω, εx), φ(x)i = limhξ˜3 (ω, x), φ( )i
ε→0 ε→0 ε ε
1 x 1 x
= limhΞ3 (ω, x)(s3 ) , φ( )i = lim(−1)s3 hΞ3 (ω, x), s3 +1 φ(s3 ) ( )i
ε→0 ε ε ε→0 ε ε
v ¯
u
u Z x ¯ Ã 3
!(i) ¯¯2
D(ω) u ¯ 2d + ¯
¯Md (t) κ(t) (x − t)
3 2
= d3 +2 tsup ¯ dt
|x| ¯
i≤d3 −2a ¯
5
Γ(2d3 + 2 ) ¯
¯
sZ
x
D(ω)
≤ K1 d3 +2 (x − t)2d3 +3 dt ≤ K2 · D(ω),
|x| −2a
Now, concerning equation (2.52), we will assume that the GRP on the
right hand side of the equation has a primitive with a point value at t = 0
in sense of Lojasiewicz.
d
y(ω, ·) = f (ω, ·) and y(ω, 0) = X(ω).
dt
Proof. Theorem 2.5.3 c) ensures the existence of F (ω, t), ω ∈ Ω, t ∈ [−2, ∞)
and k ∈ N such that
where F (ω, ·) is continuous for every ω ∈ Ω, and F (·, t) is Gaussian for every
t ∈ [−2, ∞). From the assumption on the initial point t = 0 we get that
F (ω, t)
must be bounded in a neighborhood of zero. Define
tk
µ Z t ¶(k)
X(ω)tk
y(ω, t) = + F (ω, s)ds , ω ∈ M, t ∈ [−2, ∞). (2.72)
k! 0
Note that the solution given in (2.72) coincides with the solution given
in (2.53) up to the set M , provided (2.54) holds and thus (2.53) exists.
/ Z 2.
although |X| ∈
This is a GRP (II) which does not map K{Mp } to Z 2 , since if the support of
a test function φ equals [0, a] for some a > 0, then for x1 being in this interval,
one has 1/X(x) = 1/x ∈ [0, a) and X(x) = x ∈ [1/a, ∞), thus ξ(·, φ) is not
in Z 2 .
Example 2.5.5 We construct a GRP (II) which does not satisfy (2.58),
and thus can not be represented as a sum of continuous processes on a set of
probability measure one, only on a set of arbitrary large probability measure.
74 Generalized Stochastic Processes
Example 2.5.6 Two examples of GRPs (I) which are also GRPs (II) are
singular pointwise one-dimensional one-parameter white noise, and Brown-
ian motion. Recall from Example 1.7.1 and Example 1.7.2 , they are given
respectively, by
∞
X
W (ω, x) = ξk (x)Hεk (ω), ω ∈ S0 (R), x ∈ R
k=1
and
∞ Z
X x
B(ω, x) = ξk (t)dtHεk (ω), ω ∈ S0 (R), x ∈ R
k=1 0
∞ Z Z
X x
S(R) 3 ϕ 7→ hB(ω, x), ϕ(x)i = ξk (t)dtϕ(x)dxHεk (ω).
k=1 R 0
d
Since for fixed ω ∈ Ω we have dx B(ω, x) = W (ω, x) in S0 (R) and W (ω, x)
is continuous with respect to x ∈ R (see [HØUZ]), this is indeed compatible
with the assertion in Theorem 3.5.3, i.e. continuity ensures representation
on the whole Ω.
2.6 Hilbert Space Valued GRPs (II) 75
d
W (ω, φ̃x ) = B(ω, φ̃x ), where φ̃x (y) = f1x− (y) = H(x − y).
dx
This means that in distributional sense i.e. in S0 (R) we obtain again the
d
well-known result dx B(ω, x) = W (ω, x).
For r > 1 denote Lr (Rn ; H) = Lr (Rn ) ⊗ H and recall that its dual is
L (Rn ; H), p = r/(r
p p n r n
R − 1). The dual pairing of f ∈ L (R ; H), ϕ ∈ L (R ; H)
can be written as Rn hf (t), ϕ(t)iH dt. It can easily be checked that following
H–valued version of Theorem 2.5.1 holds.
Q
Theorem 2.6.1 Let G = ni=1 (αi , βi ) ⊂ Rn , −∞ ≤ αi < βi ≤ ∞, i =
1, 2, ..., n, and let ξ be a GRP on Ω×Lr (G; H), r > 1. There exists a function
f : Ω × G → H such that
(i) for every x ∈ G, f (·, x) is measurable and for every ω ∈ Ω, f (ω, ·) ∈
Lp (G; H), p = r/(r − 1).
(ii) Z
ξ(ω, ϕ) = hf (ω, t), ϕ(t)iH dt, ω ∈ Ω, ϕ ∈ Lr (G; H).
G
Theorem 2.6.2 a) Let ξ be a H–valued GRP (II). Then for every ε > 0
there exist d ∈ N0 , M ∈ F satisfying P (M ) ≥ 1 − ε, and functions
fα : Ω × R → H, α = 0, 1, ..., d, such that fα (·, t) is measurable for
every t ∈ R, fα (ω, ·) is in L2 (R; H) for every ω ∈ M , α = 0, 1, ..., d
and
X d Z
ξ(ω, ϕ) = hfα (ω, t), Md (t)ϕ(α) (t)iH dt, ω ∈ M, ϕ ∈ K{Mp }(H),
α=0 R
(2.75)
d
X
kfα (ω, ·)kL2 (R;H) ≤ d, ω ∈ M. (2.76)
α=0
Proof. a) Since for every ω ∈ Ω, ξ(ω, ·) is in K0 {Mp }(H), it follows that for
every ω ∈ Ω there exist C(ω) > 0 and p(ω) ∈ N such that
We can assume that p(ω) ≥ C(ω). For every ϕ ∈ K{Mp }(H) and N ∈ N,
put
\
AN (ϕ) = {ω ∈ Ω : |ξ(ω, ϕ)| < N kϕkN,2;H }, AN = AN (ϕ).
ϕ∈K{Mp }(H)
Since K{M
T p }(H) is separable, it contains a countable dense subset D and
AN = ϕ∈D AN (ϕ) ∈ F. Thus, from
∞
[
Ω= AN and AN ⊂ AN +1 , N ∈ N,
N =1
it follows that for given ε > 0 there exists an integer d such that P (Ad ) ≥
1 − ε. Denote M = Ad . It follows
then ξ1 (ω, ϕν ) → 0, ν → ∞.
Qd
Let Γd = L2 (R; H) and endow it with the scalar product
Pdi=0 R
((ϕα ), (ψα )) = α=0 R hϕα , ψα iH dt, (ϕα ), (ψα ) ∈ Γd . Clearly, Γd is a
Hilbert space. Define a mapping θ : K{Mp }(H) → Γd by θ(ϕ) =
(Md ϕ, Md ϕ0 , . . . , Md ϕ(d) ), ϕ ∈ K{Mp }(H), which is injective, and denote
∆ = θ(K{Mp }(H)). Note that
d
∆ = θ(HM ). (2.81)
The existence of the limit follows from (2.80) and (2.81). Thus,
¯ ψ⊥ ∈ ∆
F (ω, ψ̃) = F (ω, ψ), ψ̃ ∈ Γd , ψ̃ = ψ + ψ ⊥ , ψ ∈ ∆, ¯ ⊥.
It is endowed with the natural norm such that it is isometric to L2 (R; H), for
every α ∈ {0, 1, ..., d}. Let ψ ∈ Γd . Denote by [ψ]α the corresponding element
in Γd,α . The βth coordinates of [ψ]α are equal to zero for β 6= α and the αth
2.6 Hilbert Space Valued GRPs (II) 79
and
d
X
kF (ω, ·)k0Γd = kfα (ω, ·)kL2 (R;H) ≤ S(ω) ≤ d, ω ∈ Ω,
α=0
where k · k0Γd
is the dual norm. Now, the assertion follows by (2.78).
The proof of the last assertion in a) follows by repeating the previous
proof starting from relation (2.79). It follows that ξ is of the form (2.56) for
every ω ∈ Ω.
b) Obviously, Cξ (ϕ, ψ) = E(hξ(·, ϕ), ξ(·, ψ)iH ), ϕ, ψ ∈ K{Mp }(H) is bi-
linear. The continuity follows from
d X
X d Z Z
= E( hfα (·, t), Md (t)ϕ(α) (t)iH hfβ (·, s), Md (s)ψ (β) (s)iH dtds)
α=0 β=0 R R
d Z Z
d X
X
= ( E(hfα (·, t), fβ (·, s)iH )Md (t)Md (s)hϕ(α) (t), ψ (β) (s)iH dtds.
α=0 β=0 R R
80 Generalized Stochastic Processes
Chapter 3
Applications to Singular
Stochastic Partial Differential
Equations
The simplest example is when the coefficients are constants, e.g. for
aij = δij (the Kronecker symbol), bi = ci = d = 0, i, j = 1, 2 . . . , n, we obtain
the Laplace operator L = ∆. Note that the operator L in (3.2) is given in
divergence form, which will make it suitable to work with in Sobolev spaces
in terms of weak derivatives, since its coefficients are singular. If its principal
coefficients aij , bi , i, j = 1, 2, . . . , n are assumed to be differentiable, then L
can also be written in form
n
X n
X
Lu= âij Dij u + ˆ
b̂i Di u + du, (3.3)
i,j=1 i=1
P P
where dˆ = d + ni=1 Di bi and b̂i = bi + ci + nk=1 Dk aki , i = 1, 2, . . . , n.
A physical interpretation of equation (3.1) in nondivergent form (3.3)
can be given as in [Ev]: u can be interpreted as the density of some quantity
(e.g. a chemical
P concentration) at equilibrium within a region I. The prin-
cipal term ni,j=1 âij Di Dj u represents the diffusion of u within I, and the
coefficients âij describe the anisotropic, heterogeneous
Pn nature of the medium
i ij
(e.g. wood or liquid crystal). The vector F = j=1 a Dj u, i = 1, 2, . . . , n
may
Pn be jinterpreted as the flux density. The ellipticity condition will imply
j=1 F Dj u ≥ 0, meaning the flow direction is from regions of lower to
higher concentration (to get the opposite direction, one hasP to put a minus
sign in front of the coefficients âij ). The first order term ni=1 b̂i Di u rep-
resents transport within I, while the zeroth order term du describes local
creation or depletion (owing, for example, to reactions).
In the framework we consider, the coefficients of L will be stochastic
processes; thus in physical interpretation equation (3.1) can be understood
as a diffusion process in a stochastic anisotropic medium, with transport
and creation also dependent on some random factors, and with a stochastic
boundary value. Example of a stochastic anisotropic medium is a medium
consisting of two randomly mixed immiscible fluids.
In neurology, elliptic PDEs showed as a good model for brain func-
tion measurements. The elliptic PDE defined in the brain at a few points
3.1 A Linear Elliptic Dirichlet Problem 83
(the
Pn location Pn ofijthe sensors) involves the divergence form operator L =
ij
i=1 D i ( j=1 a D j u). The fact that the matrix [a ] is discontinuous along
the boundaries separating different layers of the brain, means that the PDE
must deal with singular data. Now, if there are also some random perturba-
tions happening in the brain, the PDE is no longer a PDE but a SPDE.
If the coefficients of L are deterministic functions, then our assumptions
on L will follow the approach in [GT]: We assume that L is strictly elliptic in
I and has bounded measurable coefficients; that is, there exist λ > 0, Λ > 0,
and ν ≥ 0 such that
n
X
aij (x)ξi ξj ≥ λ|ξ|2 , x ∈ I, ξ ∈ Rn , (3.4)
i,j=1
n
X
|aij (x)|2 ≤ Λ2 , x ∈ I, (3.5)
i,j=1
n
1 X¡ i 2 i 2
¢ 1
2
|b (x)| + |c (x)| + |d(x)| ≤ ν 2 , x ∈ I. (3.6)
λ i=1 λ
In fact, λ and Λ are the minimal and maximal, respectively, eigenvalues of
the matrix [aij ]. Without loss of generality we may and will assume that
λ ≤ 1 (else we divide equation (3.1) by λ to obtain this situation). We
seek for a weak solution of (3.1) in the Hilbert space of generalized random
processes L(W01,2 (I), (S)−1 ), where W01,2 (I) denotes the Sobolev space and
(S)−1 denotes the Kondratiev space. As we will see in Theorem 3.1.6, this
generalized solution exists and is unique. In Proposition 3.1.1 we investigate
some stability and regularity properties of this solution. In [LØUZ] regular-
ization was used and an error-function tending to zero was introduced to give
meaning to a distributional valued solution, as well as to give meaning to the
product of a non-smooth deterministic function (the coefficients of L) and
the distribution u(x, ·). The solution was expressed in terms of the Green
function and the Itô integral. We will have a similar stability result, but
using the Hilbert space structure. The goal in this section is to extend the
approach for even more singular coefficients and data, namely Colombeau
generalized functions. We develop the necessary Sobolev–Colombeau spaces
for generalized stochastic processes and prove existence and uniqueness of
the solution for the Dirichlet problem also in this setting.
However, if the coefficients of the operator L are also generalized ran-
dom processes, then a further problem arises: How to interpret the product
between two generalized random processes? In [Va] the product was inter-
preted as the Wick product, and the solution was found as an element of
84 Applications to Singular SPDEs
for all v ∈ W01,2 (I), and u − g ∈ W01,2 (I). The concept of the weak solution
is based on the fact that one can identify the operator L with its unique
3.1 A Linear Elliptic Dirichlet Problem 85
extension (denoted by the same symbol) L : W 1,2 (I) → W −1,2 (I) which
generates a bilinear form (u, v) 7→ B(u, v) = −hLu, vi, (u, v) ∈ W 1,2 (I) ×
W01,2 (I) where hLu, vi is defined by the left hand side of (3.8). Thus, u ∈
W 1,2 (I) is a weak solution if B(u, v) = hh+∇·f, vi = hh, vi−hf, ∇vi, for every
v ∈ W01,2 (I). Here h·, ·i denotes the standard scalar product in L2 (I), which
can be extended to the canonical dual pairing on W −1,2 (I) × W01,2 (I). Thus,
existence of a weak solution is equivalent to surjectivity, while uniqueness is
equivalent to injectivity of the mapping L. The standard norm on L2 (I) is
denoted by k · k2 .
We state the following theorem from [GT], which will be needed in the
sequel:
Theorem 3.1.1 Let the operator L given by (3.2) satisfy conditions (3.4),
(3.5) and (3.6). Moreover, assume that
Z n
X
(d(x)v(x) − bi (x)Di v(x)) dx ≤ 0, v ≥ 0, v ∈ W01,2 (I). (3.9)
I i=1
where h = (h, f 1 , f 2 , . . . f n ).
These results are based on Hilbert space methods for PDEs. The Lax-
Milgram theorem and the Fredholm alternative are used to prove existence
and uniqueness of the weak solution. For details, refer to [GT, Chapter 8].
Here we state the three main classical theorems, since they are the keystones
that cap the arch of our proofs.
In case (ii), the operator (I − T )−1 whose existence is asserted there is also
bounded.
λ x − T x = y, λ x − T ∗x = y (3.11)
W −1,2 (I) ⊗ (S)−1 . The latter isomorphisms hold since (S)1 is a nuclear space
and W −1,2 (I) is a Frèchet space. Also note that the notation ⊗ stands for the
π–completition of the tensor product space, which is in this case equivalent
to the ²–completition by the nuclearity of (S)1 .
By W0 S we denote W01,2 (I) ⊗ (S)1 , and by WS we denote W 1,2 (I) ⊗ (S)1 .
We state now Theorem 2.1.1 in this new context:
Theorem 3.1.5 Following conditions are equivalent:
(i) u ∈ WS∗ .
(ii) u can be represented in the form
X
u(x, ω) = fα (x) ⊗ Hα (ω), x ∈ I, ω ∈ Ω, fα ∈ W −1,2 (I), α ∈ I,
α∈I
(3.12)
and there exists p ∈ N0 such that for each bounded set B ⊂ W01,2 (I)
X
sup |hfα , ϕi|2 (2N)−pα < ∞. (3.13)
ϕ∈B
α∈I
(iii) u can be represented in the form (3.12) and there exists p ∈ N0 such
that X
kfα k2W −1,2 (2N)−pα < ∞. (3.14)
α∈I
P
Note that α∈I fα (x)Hα (ω), fα ∈ W ±1,2 is also a generalized random
process as defined in [Va]. The tensor product space W −1,2 (I) ⊗ (S)−1 was
used also in [Be]; we find it more convenient to deal with the structure given in
Theorem 3.1.5 and also later in the next section to deal with Wick products,
which are all special cases (Sobolev versions) of the structures introduced in
[PS1] and [Se].
Proof. Let θ ∈ (S)1 and p > 0. Put gθ (x) = hg(x, ω), θ(ω)i, hθ (x) =
hh(x, ω), θ(ω)i, fθi (x) = hf i (x, ω), θ(ω)i, i = 1, 2, . . . , n. Clearly, gθ ∈
W 1,2 (I), while hθ , fθi ∈ L2 (I), i = 1, 2, . . . , n. According to Theorem 3.1.1
there exists a unique weak solution uθ (x), x ∈ I in W 1,2 (I) of the boundary
problem
n
X
L uθ (x) = hθ (x) + Di fθi (x), x ∈ I,
i=1 (3.15)
uθ (x) ¹∂I = gθ (x).
Due to (3.10) and continuity of h and g, there exist constants C > 0, M > 0
(note: C may depend on the coefficients of L and on I, but not on θ) such
that
Due to Theorem 3.1.1, for every α ∈ I there exists a unique uα ∈ W 1,2 (I),
which solves (3.17), and there exists C > 0 (uniformly in α) such that
¡ ¢
kuα k2W 1,2 ≤ C khα k22 + kgα k2W 1,2
where hα = (hα , fα1 , fα2 , . . . fαn ). Now, according to the assumptions made on
h, f i , i = 1, 2 . . . , n, there exists p > 0 such that
à !
X X X
kuα k2W 1,2 (2N)−pα ≤ C khα k22 (2N)−pα + kgα k2W 1,2 (2N)−pα < ∞.
α∈I α∈I α∈I
Thus, u ∈ WS. Again, one can consider it also as an element of the dual
space i.e. u ∈ WS∗ .
90 Applications to Singular SPDEs
Corollary 3.1.1 Let u ∈ WS∗ be the generalized solution of (3.1). Then its
generalized expectation E(u) coincides with the weak solution of the deter-
ministic Dirichlet problem
n
X
L v(x) = h̃(x) + Di f˜i (x), x ∈ I,
i=1 (3.18)
v(x) ¹∂I = g̃(x),
Proof. The assertion follows from the proof of Theorem 3.1.6 if we choose
θ = 1. ¤
every ϕ ∈ W01,2 (I) the expression |huθ − ũθ , ϕi| is bounded by the operator
norms of kL − L̃k and khθ − h̃θ k, where h = (h, f 1 , f 2 , . . . f n ). Here h·, ·i
stands for the usual scalar product in L2 (I).
By subtracting (3.20) from (3.19) we obtain
n
X
Luθ − L̃ũθ + L̃uθ − L̃uθ = hθ − h̃θ + Di (fθi − f˜θi ),
i=1
that is
n
X
L̃(uθ − ũθ ) = hθ − h̃θ + Di (fθi − f˜θi ) − (L − L̃)uθ .
i=1
|huθ − ũθ , ϕi| = |huθ − ũθ , L̃∗ ψi| = |hL̃(uθ − ũθ ), ψi|
n
X
= |hhθ − h̃θ + Di (fθi − f˜θi ) − (L − L̃)uθ , ψi|.
i=1
n
X
|huθ − ũθ , ϕi| ≤ khθ − h̃θ k2 kψk2 + kfθi − f˜θi k2 kDi ψk2 + kL − L̃k2 kuθ k2 kψk2
i=1
≤ khθ − h̃θ k2 kψkW 1,2 + kL − L̃kW −1,2 kuθ kW 1,2 kψkW 1,2 .
(3.21)
From Theorem 3.1.1 we also know there exists C > 0 such that
Thus,
where h = (h, f 1 , f 2 , . . . f n ).
92 Applications to Singular SPDEs
Now, consider the general Dirichlet problems (3.19), (3.20) with nonzero
boundary conditions. We transform them into zero boundary conditions in
the following way: uθ is a weak solution
Pn of (3.19) if and only if uP
0,θ = uθ −gθ is
a weak solution ofPLu0,θ = h0,θ − i=1 Di f0,θ , where h0,θ = hθ − ni=1 ci Di gθ −
i
i
dgθ , f0,θ = fθi − ni=1 aij Dj gθ − bi gθ . We do the same for (3.20) and apply
the estimate obtained for |hu0 − ũ0 , ϕθi|.
We summarize the stability result in the following proposition.
Proposition 3.1.1 Let L be a strictly elliptic operator of the form (3.2) with
coefficients aij , bi , ci , d, satisfying conditions (3.4), (3.5) and (3.6). Let L̃
be another strictly elliptic operator of the form (3.2) with coefficients ãij , b̃i ,
˜ satisfying all given conditions. Let h, h̃, f i and f˜i , i = 1, 2, . . . , n be
c̃i , d,
generalized random processes from L((S)1 , L2 (I)). Let g and g̃ be generalized
random processes from WS. Let u, ũ ∈ WS∗ be the generalized solutions of
the Dirichlet problems
n
X
L u(x, ω) = h(x, ω) + Di f i (x, ω), x ∈ I, ω ∈ Ω,
i=1
u(x, ω) ¹∂I = g(x, ω)
and
n
X
L̃ ũ(x, ω) = h̃(x, ω) + Di f˜i (x, ω), x ∈ I, ω ∈ Ω,
i=1
ũ(x, ω) ¹∂I = g̃(x, ω)
respectively. There exist C > 0, p ∈ N0 such that for every θ ∈ (S)1,p :
³
|hu(x, ω) − ũ(x, ω), ϕ(x)θ(ω)i| ≤ C kh − h̃kL2 (I)⊗(S)−1,−p
+kL − L̃kW −1,2 ku − gkW 1,2 ⊗(S)−1,−p + kL − L̃kW −1,2 kg̃kW 1,2 ⊗(S)−1,−p
´
+kg − g̃kW 1,2 ⊗(S)−1,−p kLkW −1,2 kϕkW 1,2 kθk1,p
where h = (h, f 1 , f 2 , . . . f n ).
Remark. In particular, to address a similar stability question as in
[LØUZ], let us consider the net of operators L² , ² ∈ (0, 1], given by
n
X n
X
L² u(x, ·) = Di ( aij i
² (x)Dj u(x, ·) + b² (x)u(x, ·))
i=1 j=1
n
X
+ ci² (x)Di u(x, ·) + d² (x)u(x, ·)
i=1
3.1 A Linear Elliptic Dirichlet Problem 93
and a net of data h² (x, ω), f²i (x, ω), i = 1, 2, . . . , n, where aij ij
² (x) = a ∗ ρ̌² (x),
i i i i
b² (x) = b ∗ ρ̌² (x), c² (x) = c ∗ ρ̌² (x), d² (x) = d ∗ ρ̌² (x),h² (x, ω) = h(·, ω) ∗ ρ̌² (x),
f²i (x, ω) = f i (·, ω)∗ ρ̌² (x), i, jR = 1, 2, . . . , n. Here ρ is a mollifier, i.e. a positive
smooth function such that I ρ(x)dx = 1, ρ̌(x) = ρ(−x), ρ² (x) = ²−n ρ(x/²),
and ∗ denotes convolution with respect to x. Denote by u² the solution of
n
X
L² u² = h² + Di f²i , u² ¹∂I = 0.
i=1
Now, from (3.21) we get that kuθ,² − uθ kW −1,2 is bounded by the sum of the
operator norm kL − L² kW −1,2 and of khθ,² − hθ kL2 . It can be shown (see
[NP]) that L² inherits the strict ellipticity and boundedness properties from
L. Also, it is a well-known property of regularization in Sobolev spaces that
kp ∗ ρ² (x) − p(x)kW 1,2 → 0, ² → 0, for p ∈ W 1,2 (I).
Thus,
kuθ,² − uθ kW −1,2 → 0, ² → 0, (3.22)
which establishes the stability result.
This motivates us to consider the stochastic Dirichlet problem in the
Colombeau setting, which will enable us to solve problem (3.1) also with more
singular coefficients and data, i.e. with Colombeau generalized deterministic
functions. In order to do this, we first have to introduce appropriate spaces
of Colombeau algebras of generalized random processes, which will be done
in the next section.
Proposition 3.1.1 and the Remark after it also relate our approach with
those in [HØUZ] and [LØUZ]. Processes involved therein as data in SPDEs
are first smoothed out in order to get an appropriate form for applying the
S-transform and its inverse. For example, singular white noise given by
∞
X
W (x, ω) = ηj (x)Hεj (ω), x ∈ Rn , ω ∈ Ω,
j=1
involving smoothed white noise are then solved by standard methods, and
then letting ϕ → δ in S0 (Rn ) one obtains a solution of the SPDE involving
singular white noise. Due to Theorem 3.1.6 we can directly solve the Dirich-
let problem (3.1) where h, f or g is singular white noise, since singular white
noise belongs to L((S)1 , L2 (I)).
For example, consider the Dirichlet problems Lu = W (ρ² , x, ω) and Lu =
W (x, ω), u ¹∂I = 0, with generalized solutions u² , u, respectively. Now, using
the fact that W (ρ² , x, ω) = W (·, ω) ∗ ρ² (x), we get from the Remark after
Proposition 3.1.1
This stability result shows that the solution of the Dirichlet problem
where the right hand side is smoothed white noise, in fact is an adequate
approximation of the solution where the data is singular white noise.
Lu = h, u ¹∂I = g.
where K = max{kaij , bi kC 0,1 (I) , kci , dkL∞ (I) } and d0 = dist(∂I, I 0 ). Addi-
tionally, u satisfies the equation
n
X n X
X n n
X
ij ij i i
Lu= a Dij u + ( (Dj a + b + c )Di u + ( Di bi + d)u = h
i,j=1 i=1 j=1 i=1
for a.e. x ∈ I.
3.1 A Linear Elliptic Dirichlet Problem 95
The next theorem (in deterministic case) was proved in [MP], but we may
state it also in our setting. The polynomial growth rate estimated in (3.23)
will later be crucial for the Colombeau solutions.
• Let N(W 2,2 ) denote the vector space of functions r² ∈ EM (W 2,2 ) with
the property that for every a ∈ N0 , kr² (x)kW 2,2 = O(²a ) holds. Ele-
ments of N(W 2,2 ) are called negligible functions.
• Since N(W 2,2 ) is an ideal of the algebra EM (W 2,2 ), the quotient space
{aα }α∈I bounded from above (i.e. there exists a ∈ R such that aα ≤ a,
α ∈ I), with following properties:
X
Cα2 (2N)−pα < ∞. (3.25)
α∈I
i.e. we use the Wick product for multiplication in (S)−1 and the ordi-
nary product in W 2,2 (I).
98 Applications to Singular SPDEs
Proof. (i) Let F² , G² ∈ EMP (W 2,2 ; (S)−1 ) and prove that F² ♦G² ∈
2,2 2,2
EM (W ; (S)P −1 ). Since F² = α∈I fα,² (x) ⊗ Hα (ω) ∈ EM (W ; (S)−1 )
2,2
and G² = g
β∈I β,² (x) ⊗ H β (ω) ∈ E M (W ; (S)−1 ), there exist
a, b ∈ R, p, q ∈ N0 , ²1 , ²2 ∈ (0, 1), and there exist sequences {Cα }α∈I ,
{Dα }α∈I , {aα }α∈I , {bα }α∈I , such that aα ≤ a, bα ≤ b for all α ∈ I,
kf
Pα,² kW 2,2 ≤ Cα ²−aα forP² < ²1 , kgα,² kW 2,2 ≤ Dα ²−bα for ² < ²2 , and
2 −pα
α∈I Cα (2N) < ∞, α∈I Dα2 (2N)−qα < ∞. We will prove that for
²0 = min{²1 , ²2 } there exist r > 0, a sequence P {Mγ }γ∈I and a sequence
{mγ }γ∈I bounded from above such that k α+β=γ fα,² (x)gα,² (x)kW 2,2 ≤
P
Mγ ²−mγ for ² < ²0 and γ∈I Mγ2 (2N)−rα < ∞.
P
For a fixed multiindex γ ∈ I put Mγ = α+β=γ Cα Dβ and mγ =
max{aα + bβ : α, β ∈ I, α + β = γ} (note, for γ fixed, there are only
finite many α and β which give the sum γ). Now using the fact that
W 2,2 (I) is an algebra for n ≤ 3, we get
X X
k fα,² (x)gα,² (x)kW 2,2 ≤ kfα,² (x)kW 2,2 kgα,² (x)kW 2,2 ≤
α+β=γ α+β=γ
X X
Cα Dβ ²−(aα +bβ ) ≤ ²−mγ Cα Dβ = ²−mγ Mγ .
α+β=γ α+β=γ
P 2 −qα
DP
α∈Iα (2N) < ∞. Let ²0 = min{²1 , ²2 }, r = p + q + 2 and
Mγ = α+β=γ α Dβ , γ ∈ I. Now we may
C proceed as in (i) to get
P mγ
P
α+β=γ kfα,² kW 2,2 kgα,² kW 2,2 ≤ Mγ ² and γ∈I Mγ2 (2N)−rα < ∞.
In order to prove that N(W 2,2 ; (S)−1 ) is an ideal of EM (W 2,2 ; (S)−1 ),
we must check that for all G² ∈ N(W 2,2 ; (S)−1 ) and all F² ∈
EM (W 2,2 ; (S)−1 ), we have G² ♦F² = F² ♦G² ∈ N(W 2,2 ; (S)−1 ). Let
{nγ }γ∈I be an arbitrary sequence bounded from below (i.e. nγ ≥ n,
γ ∈ I). Since F² ∈ EM (W 2,2 ; (S)−1 ), there exist p ∈ N0 , ²1 ∈ (0, 1),
{Dβ }β∈I , b ∈ R and a sequenceP{bβ }β∈I such that bβ ≤ b, β ∈ I,
kfβ,² kW 2,2 ≤ Dβ ²−bβ , ² < ²1 , and β Dβ2 (2N)−pβ < ∞.
For a fixed multiindex α ∈ I, let aα = b + nα . Clearly, the se-
quence {aα }α∈I is bounded from below by a = b + n, thus since
G² ∈ N(W 2,2 ; (S)−1 ), there exist q ∈ N0 , P
²2 ∈ (0, 1), {Cα }α∈I , such
−aα 2 −qα
that kgα,² kW 2,2 ≤ Cα ² ,P² < ²2 , and α Cα (2N) < ∞. Let
²0 = min{²1 , ²2 } and Nγ = α+β=γ Cα Dβ , γ ∈ I.
P P
Now for G² ♦F² = γ∈I α+β=γ gα,² fα,² ⊗ Hγ we have that
X X
kgα,² kW 2,2 kfβ,² kW 2,2 ≤ Cα Dβ ²aα −bβ
α+β=γ α+β=γ
X X
≤ Cα Dβ ²b+nα −bβ ≤ Cα Dβ ²bβ +nα −bβ ≤ Nγ ²nγ , ² < ²0 .
α+β=γ α+β=γ
P
It is clear that γ∈I Nγ2 (2N)−rγ < ∞ for r = p + q + 2. Thus, G² ♦F² ∈
N(W 2,2 ; (S)−1 ). ¤
• Since N(W 2,2 ; (S)−1 ) is an ideal of the algebra EM (W 2,2 ; (S)−1 ), the
quotient space
(iii) R² (x, ω), x ∈ I, ω ∈ Ω, belongs to N(W 2,2 ; (S)−1 ) if and only if there
exist ²0 ∈ (0, 1), C > 0 and q ≥ 0 such that for every a ∈ R,
X Z
kDβ R² (x, ·)k2−1,−q dx ≤ C²a , ² < ²0 . (3.31)
0≤|β|≤2 I
(iv) R² (x, ω), x ∈ I, ω ∈ Ω, belongs to N(W 2,2 ; (S)−1 ) if and only if there
exist ²0 ∈ (0, 1), C > 0 and q ≥ 0 such that for every a ∈ R,
XZ X X
= |Dβ rα,² (x)|2 dx(2N)−pα = krα,² k2W 2,2 (2N)−pα
α∈I I 0≤|β|≤2 α∈I
X X
≤ Cα2 ²−2bα (2N)−pα ≤ ²−2b Cα2 (2N)−pα = C²−a .
α∈I α∈I
3.1 A Linear Elliptic Dirichlet Problem 101
à !4 à !2
X X
krα,² kW 2,2 (2N)−qα ≤ K krα,² k2W 2,2 (2N)−qα
α∈I α∈I
2
Z X X
≤ K2 |Dβ rα,² (x)|2 (2N)−qα dx
I α∈I 0≤|β|≤2
2
Z X X
≤ K 2 |I| |Dβ rα,² (x)|2 (2N)−qα dx
I α∈I 0≤|β|≤2
Z Ã !2
X X
2 β 2 −qα
≤ 2K |I| |D rα,² (x)| (2N) dx
I 0≤|β|≤2 α∈I
Z X
= K̃ kDβ R² (x, ·)k2−1,−q dx ≤ K̃C²−a .
I 0≤|β|≤2
P p
4 a
Thus, α∈I krα,² kW 2,2 (2N)
−qα
≤ K̃C²− 4 , which implies that for each α ∈ I
a p
krα,² kW 2,2 ≤ Cα ²− 4 , where Cα = (2N)qα 4
K 2 |I|C,
and X
Cα2 (2N)−pα < ∞, for p > 1 + 2q.
α∈I
(ii) Since
kR² (x, ω)k2W 2,2 (I)⊗(S)−1,−q = hR² (x, ω), R² (x, ω)i
X X
= hrα,² , rα,² iW 2,2 hHα , Hα i−1,−q = krα,² k2W 2,2 (2N)−qα ≤ C²−a ,
α∈I α∈I
it follows similarly as in (i) that krα,² k2W 2,2 ≤ C(2N)qα ²−a for each α ∈ I.
(iii) and (iv) follow similarly as (i) and (ii). ¤
102 Applications to Singular SPDEs
2. there exist M > 0 and b > 0 such that for all ² ∈ (0, 1) and α ∈ I:
kaij i i
² kW 2,2 , kb² kW 2,2 , kc² kW 2,2 , kd² kW 2,2 ≤ Λ² ;
Λ² ≤ M ²b .
104 Applications to Singular SPDEs
Using the estimate derived in (3.23) we obtain that there exist C > 0,
s > 0 (from the proof in [MP] one can see that C and s may depend on n, λ,
|I|, ∂I, and the coefficients of L, but they are uniform in ² and α) such that:
µ ¶s ³ ´
Λ²
kuα,² kW 2,2 ≤ C sup |gα,² (x)| · |I| + kgα,² kW 2,2 + khα,² kW 0,2 . (3.39)
λ² x∈∂I
Now, since gα,² and hα,² , α ∈ I, are all bounded by Λ² , which has polyno-
mial growth with respect to ², we obtain
and
L2,² u² (x, ω) = h2,² (x, ω) in W 2,2 (I) ⊗ (S)−1
(3.41)
u² (x, ω) ¹∂I = g2,² (x, ω) ¹∂I , ² ∈ (0, 1),
Thus,
∇ · (A1,² · ∇(u2,² − u1,² ) + b1,² · (u2,² − u1,² )) + c1,² · ∇(u2,² − u1,² ) + d1,² (u2,² − u1,² )
= h2,² − h1,² − (∇ · ((A2,² − A1,² ) · ∇u1,²
+ (b2,² − b1,² ) · u1,² ) + (c2,² − c1,² ) · ∇u1,² + (d2,² − d1,² )u1,² )
= H²
(u2,² − u1,² ) ¹∂I = (g2,² − g1,² ) ¹∂I , ² ∈ (0, 1),
where H² ∈ N(W 2,2 ; (S)−1 ). Using the estimate as in (3.39) we finally obtain
that there exist p > 1 and C > 0 such that for all a ∈ R
X X
kuα,2,² (x) − uα,1,² (x)k2W 2,2 (2N)−pα ≤ C 2 ²2a (2N)−pα < ∞.
α∈I α∈I
106 Applications to Singular SPDEs
P P
α∈I fα (x)Hα (ω), fα ∈ W01,2 (I) such that α∈I kfα k2W 1,2 (2N)pα < ∞, the
dual action is: X
hF, f i = α!hFα , fα iL2 (I) .
α∈I
qP
Clearly, the norm corresponding to h·, ·i is kF k = α∈I α!kFα k2L2 (I) .
Definition
P 3.2.2 Let F be an essentially bounded GRP given by F (x, ω) =
f
α∈I α (x) ⊗ Hα (ω), fα ∈ L∞ (I). Let G ∈ WS∗ be given by its chaos
108 Applications to Singular SPDEs
P
expansion G(ω) = β∈I gβ (x) ⊗ Hβ (ω), gβ ∈ W −1,2 (I). The Wick product of
F and G is the unique element in WS∗ defined by:
à !
X X
F ♦G(x, ω) = fα (x)gβ (x) ⊗ Hγ (ω).
γ∈I α+β=γ
The next lemma shows that the Wick product is well defined, and that
for fixed F the mapping G 7→ F ♦G is continuous.
Proof.
Since F is essentially bounded and G is an element of WS∗ , there exist
q > 0 and r > 0 such that
X X
kfα k2L∞ (2N)−qα < ∞ and kgβ k2W −1,2 (2N)−rα < ∞.
α β
By definition of the
P Wick −kγ
product and using the Cauchy-Schwartz inequality,
we have for C = γ (2N) and p = q + r + k, for arbitrary k > 1
° °2
X°
° X °
°
° fα (x)g β (x) ° (2N)−pγ
° °
γ α+β=γ W −1,2
à !à !
X X X
≤ (2N)−kγ (2N)−(q+r)γ kfα k2L∞ (I) kgβ k2W −1,2
γ α+β=γ α+β=γ
à !à !
X X X
≤ (2N)−kγ kfα k2L∞ (I) (2N)−qα kgβ k2W −1,2 (2N)−rβ .
γ α β ¤
GRPs, and will thus interpret the action of L onto a GRP u ∈ WS∗ as:
3.2 A Dirichlet Problem with Stochastic Coefficients 109
n
X n
X
L♦u(x, ω) = Di ( aij (x, ω)♦Dj u(x, ω) + bi (x, ω)♦u(x, ω))
i=1 j=1
n
(3.42)
X
+ ci (x, ω)♦Di u(x, ω) + d(x, ω)♦u(x, ω).
i=1
à !
X X
Di f ♦g(x, ω) = Di fα (x)gβ (x) ⊗ Hγ (ω)
γ α+β=γ
à !
X X
= Di (fα (x)gβ (x)) ⊗ Hγ (ω)
γ α+β=γ
à !
X X
= Di fα (x)gβ (x) + fα (x)Di gβ (x) ⊗ Hγ (ω)
γ α+β=γ
In further notation we will suppress the indexes, and simply denote h·, ·i
and k · k. From the context it is clear whether the index is L∞ (I) ⊗ (S)−1,−p
or L2 (I) ⊗ (S)±1,±p or L2 (I) ⊗ (S)±1,±2(p+1) .
Due to the previous definition we are able to develop something that
might be called a weak Wick calculus: The Wick-adjoint operator inherits
its properties from the classical Wick multiplication – most important, the
chain rule holds in weak sense.
Later we will see that all properties of the operator L (ellipticity and
other assumptions on its coefficients) can be carried over to the formal Wick–
adjoint of L.
n
X
hd♦u, vi − hbi ♦u, Di vi ≤ 0. (3.46)
i=1
Conditions (3.45) and (3.46) will be adequate to prove existence and
uniqueness of the generalized weak solution of (3.47). Later on we will imply
a stronger condition to (3.45) and a weaker condition to (3.46) in order to
prove stability and regularity properties.
We turn now back to our stochastic Dirichlet problem
n
X
L♦u(x, ω) = h(x, ω) + Di f i (x, ω), x ∈ I, ω ∈ Ω,
i=1 (3.47)
u(x, ω) ¹∂I = g(x, ω)
where the action L♦u is defined as in (3.42). Applying partial integration
we obtain
X n Xn n
X
ij i
hL♦u, vi = hDi ( a ♦Dj u + b ♦u), vi + hci ♦Di u, vi + hd♦u, vi
i=1 j=1 i=1
112 Applications to Singular SPDEs
n X
X n n
X n
X
ij i
=− ha ♦Dj u, Di vi − hb ♦u, Di vi + hci ♦Di u, vi + hd♦u, vi.
i=1 j=1 i=1 i=1
for u ∈ WS∗ , v ∈ W0 S.
As in [GT], we associate a bilinear form B : WS∗ × W0 S → R associated
with L, defined by
B(u, v) = −hL♦u, vi
n X
X n n
X n
X
ij i
= ha ♦Dj u, Di vi + hb ♦u, Di vi − hci ♦Di u, vi − hd♦u, vi.
i=1 j=1 i=1 i=1
∗
We
Pncall u i∈ WS a generalized
Pn weak solution of (3.47) if hL♦u, vi =
i
hh + i=1 Di f , vi = hh, vi − i=1 hf , Di vi for all v ∈ W0 S, i.e. if
n
X
B(u, v) = −hh, vi + hf i , Di vi
i=1
for all v ∈ W0 S.
n X
X n Xn Xn
≤ kaij kkDj ukkDi vk+ kbi kkukkDi vk+ kci kkDi ukkvk+kdkkukkvk
i=1 j=1 i=1 i=1
à n !à n
!
X X
≤C kDi uk kDi vk = CkukWS∗ kvkW0 S ,
i=0 i=0
Proof. (i) Let B(u, v) ≤ 0. Then, using (3.46) and boundedness of the
coefficients ci we obtain
Xn X n n
X
ij
ha ♦Dj u, Di vi ≤ hci ♦Di u, vi
i=1 j=1 i=1
n
(3.48)
X
≤C kDi uk kvk
i=1
and note that v ≥ 0 (since vα (x) = max{uα (x) − supx∈∂I u+α (x), 0} ≥ 0 for
each α ∈ I), and Di v = Di u, i = 1, 2, . . . , n. Now from (3.48) and the
ellipticity condition (3.44) we retain that
n
X n X
X n
2
λ kDi vk ≤ haij ♦Dj v, Di vi ≤ 0,
i=1 i=1 j=1
implies
n
X
kDi vα kL2 (I) ≤ C̃kvα kL2 (I) , α ∈ I, (3.49)
i=1
for an appropriate C̃ > 0. But now, using Sobolev type inequalities one can
prove that (3.49) implies that supp Dvα is a set of strictly positive measure.
This is a contradiction with Duα = 0. Thus, supx∈I uα ≤ supx∈∂I u+ α for all
α ∈ I. This proves the assertion. ¤
and at the same time replace (3.45) by a stronger condition: Let there exist
p, Λ, ν > 0 such that for all x ∈ I,
n
X
kaij (x, ·)k2−1,−p ≤ Λ2 and
i,j=1
n
1 X¡ i 2 i 2
¢ 1
2
kb (x, ·)k −1,−p + kc (x, ·)k −1,−p + kd(x, ·)k−1,−p ≤ ν 2 . (3.52)
λ i=1 λ
of the L♦ operator:
P P
where h̃ = h − ni=1 ci ♦Di g − d♦g and f˜i = f i − nj=1 aij ♦Dj g − bi ♦g,
i = 1, 2, . . . , n. Clearly, ũ ¹∂I = 0. Thus, any stochastic Dirichlet problem of
the form (3.47) can be reduced to the zero boundary condition. Moreover, if
h, f i ∈ L2 (I) ⊗ (S)−1 and g ∈ W 1,2 (I) ⊗ (S)−1 , then h̃, f˜i ∈ L2 (I) ⊗ (S)−1
and ũ ∈ W01,2 (I) ⊗ (S)−1 .
The following two lemmas state that the bilinear form B associated with
the operator L can be made coercive by adding a sufficiently large multiple
of the P identity operator to it. (Recall, Du denotes the total differential
Du = ni=1 Di u.)
Lemma 3.2.6 Let L satisfy conditions (3.44) and (3.45). There exist con-
stants K1 , K2 > 0 such that
Lσ ♦u = L♦u − σu
is compact.
Proof. From [GT] we know that the embedding W 1,2 (I) → L2 (I) is com-
pact. From [HØUZ] we have that the embeddings (S)1,p → (L)2 are of
Hilbert-Schmidt type for all p > 0. Thus, the embedding of their pro-
jective limit (S)1 → (L)2 is compact. It can be extended to a com-
pact embedding (S)1 → (S)−1 . From these we get that the embedding
W 1,2 (I) ⊗ (S)1 → L2 (I) ⊗ (S)−1 is also compact. ¤
Now we are ready to state the main theorem about the existence and
uniqueness of generalized weak solutions of the stochastic Dirichlet problem.
Theorem 3.2.2 Let the operator L satisfy conditions (3.44), (3.45) and
(3.46). Then for h, f i ∈ L2 (I) ⊗ (S)−1 , i = 1, 2, . . . , n and for g ∈
W 1,2 (I) ⊗ (S)−1 the stochastic Dirichlet problem (3.47) has a unique gen-
eralized weak solution u ∈ WS∗ .
118 Applications to Singular SPDEs
Proof. Choose σ > 0 such that the bilinear form Bσ is coercive (this can
be done due to Lemma 3.2.7). From Lemma 3.2.5 we know that Bσ is also
continuous. According to the Lax–Milgram theorem, for Lσ there exists the
inverse operator
(Lσ ♦)−1 : WS∗ → WS
and it is a continuous, injective mapping.
The mapping F : W0 S → R, v 7→ hF, vi = hf i , Di vi − hh, vi is linear and
continuous, i.e. F ∈ WS∗ .
The equation L♦u = F , u ∈ W0 S, F ∈ WS∗ is equivalent to Lσ ♦u+σIu =
F , i.e. to
u − (−σ(Lσ ♦)−1 I)u = (Lσ ♦)−1 F.
Due to Lemma 3.2.8 the operator T = −σ(Lσ ♦)−1 I is compact. Thus, by
the Fredholm alternative,
Case (i) is impossible, since the homogeneous equation has a unique trivial
solution according to Theorem 3.2.1. Thus, case (ii) must hold true. ¤
where h = (h, f 1 , f 2 , . . . , f n ).
for u ∈ W0 S, where aji~ , bi~ , ci~ and d~ are the Wick–adjoint multiplication
operators (recall Definition 3.2.3) of the coefficients aji , bi , ci and d.
Proof.
n
X Xn n
X
~ ji~ i~
hu, L vi = hu, Di ( a Dj v − c v)i − hu, bi~ Di vi + hu, d~ vi
i=1 j=1 i=1
n
X Xn X n
=− hDi u, aji~ Dj v − ci~ vi − hu, bi~ Di vi + hu, d~ vi
i=1 j=1 i=1
Xn n
X n
X n
X
ji i
=− h a ♦Di u, Dj vi + hc ♦Di u, vi − hbi ♦u, Di vi + hd♦u, vi
i=1 j=1 i=1 i=1
= hL♦u, vi. ¤
120 Applications to Singular SPDEs
Theorem 3.2.3 Let the operator L satisfy conditions (3.44) and (3.45).
There exists a countable, discrete set Σ ⊂ R with following properties:
P
(i) If σ P∈/ Σ, the Dirichlet problems Lσ ♦u = h + ni=1 Di f i and L~ σu =
h + ni=1 Di f i , u ¹∂I = g, are uniquely solvable in WS for arbitrary
h, f i ∈ L2 (I) ⊗ (S)−1 and g ∈ WS. Moreover, there exists C > 0
(depending on L, I and σ) such that kuk ≤ C(khk + kgk).
n
X
L̃♦ũ = h̃ + Di f˜i , ũ ¹∂I = g̃, (3.63)
i=1
respectively. There exist C > 0, p ∈ N0 such that for every v ∈ W01,2 (I) ⊗
(S)1,p following estimate holds:
³
|hu−ũ, vi| ≤ C kh−h̃kL2 (I)⊗(S)−1,−p +kL−L̃kW −1,2 ⊗(S)−1,−p ku−gkW 1,2 ⊗(S)−1,−p
where h = (h, f 1 , f 2 , . . . f n ).
˜
Note that kL − L̃k = max1≤i,j≤n {kaij − ãij k, kbi − b̃i k, kci − c̃i k, kd − dk}.
Since v is arbitrary, we finally obtain
³ ´
ku − ũk ≤ K kh − h̃k + kL − L̃kkuk . ¤
Now, from Theorem 3.2.4 we get that ku² − ukWS∗ is bounded by the sum of
the operator norm kL − L² kWS∗ and of kh² − hkL2 (I)⊗(S)−1 . Thus,
ku² − ukWS∗ → 0, ² → 0.
processes i.e. elements of the algebra G(W 2,2 ; (S)−1 ). But in order to do this,
first we have to establish some necessary polynomial growth rate estimates
for the generalized weak solution u of (3.47).
In order to get these polynomial growth rates we need stricter assump-
tions on L: Assume that beside the ellipticity condition (3.44), conditions
(3.51) and (3.52) hold. By careful investigation of the proofs in [GT], we
see that one may carry over the regularity properties also in our setting to
the Hilbert space WS∗ . One only needs to consider (S)−1 –valued Sobolev
functions u(x, ω) ∈ W 2,2 (I) ⊗ (S)−1 , but this is no problem since the chain
rule and all other necessary tools hold (one can define differential quotients
and carry out the calculations as in [GT]). We state now our analogies of
[GT, Theorem 8.8.] and [GT, Theorem 8.12.] without proof.
L♦u = h, u ¹∂I = g.
for n ≤ 3 and β < 12 . Thus, kukL2 (I)⊗(S)−1,−p ≤ supx∈I ku(x, ·)k−1,−p |I|. Com-
bining the results of Theorem 3.2.6 and Theorem 3.2.7 we get the existence
of s > 0 and C > 0 such that
³ Λ ´s ³
kukW (I)⊗(S)−1,−p ≤ C
2,2 sup ku(x, ·)k−1,−p |I| + kgkW 2,2 (I)⊗(S)−1,−p
λ x∈I
´
+ khkL2 (I)⊗(S)−1,−p + khkLq (I)⊗(S)−1,−p .
(3.66)
Proposition 3.2.1 If u² ∈ N(W 2,2 ; (S)−1 ), then L² ♦u² ∈ N(W 0,2 ; (S)−1 ).
Theorem 3.2.8 Let h and g belong to G(W 2,2 ; (S)−1 ). Let the operator L =
[L² ] be given by (3.68) with coefficients [aij i i
² ], [b² ], [c² ], [d² ] in G(W
2,2
; (S)−1 ).
Assume that following conditions hold:
1. L is strictly elliptic, i.e. (3.69) holds,
2. there exist M > 0 and b > 0 such that for all ² ∈ (0, 1):
kaij i i
² k, kb² k, kc² k, kd² k, kh² k, kg² k ≤ Λ² ≤ M ²
b
where k · k denotes the norm in W 2,2 ⊗ (S)−1,−p for some p > 0 fixed.
Proof. For ² ∈ (0, 1) fixed, there exists a unique generalized weak solution
u² ∈ W 2,2 (I) ⊗ (S)−1 of problem (3.71). This follows from Theorem 3.2.2
and Theorem 3.2.6.
Using the estimate derived in (3.66) we obtain that there exist C > 0,
s > 0 (independent of ²) such that:
µ ¶s ³
Λ²
ku² kW 2,2 ⊗(S)−1,−p ≤ C sup kg² (x, ·)k−1,−p |I| + kg² kW 2,2 ⊗(S)−1,−p
λ² x∈∂I
´
+ kh² kL2 ⊗(S)−1,−p .
(3.72)
Now, since g² and h² , are all bounded by Λ² , which has polynomial growth
with respect to ², we obtain
and
n ³X
X n ´
(aij
1,² ♦Dj (u2,² − u1,² ) + bi1,² ♦(u2,² − u1,² )
i=1 j=1
Xn
+ ci1,² ♦Di (u2,² − u1,² ) + d1,² ♦(u2,² − u1,² )
i=1
³X
n n
¡X ¢
= h2,² − h1,² − (aij ij i i
2,² − a1,² )♦Dj u1,² + (b2,² − b1,² )♦u1,²
i=1 j=1
n
X ´
+ (ci2,² − ci1,² )♦Di u1,² + (d2,² − d1,² )♦u1,²
i=1
= H²
(u2,² − u1,² ) ¹∂I = (g2,² − g1,² ), ² ∈ (0, 1),
where H² ∈ N(W 2,2 ; (S)−1 ). Using the estimate (3.66) we finally obtain
that there exist p > 0 and C > 0 such that for all a ∈ R
ku2,² − u1,² kW 2,2 ⊗(S)−1,−p ≤ C²a .
Thus, by Proposition 3.1.4 (iv), (u2,² − u1,² ) ∈ N(W 2,2 ; (S)−1 ). ¤
d
∂ 1
• f (y) = yj fˆ(y), j = 1, 2, . . . n,
∂xj i
d 1 d∂
• x j f (y) = f (y), j = 1, 2, . . . n,
i ∂yj
\
• f (x − h)(y) = e−ihy fˆ(y), h ∈ Rn , (translation formula)
y
• f\
(λx)(y) = |λ|−n fˆ( ), λ ∈ C. (dilatation formula)
λ
In the Hilbert space valued case, the Fourier transformation on S0 (Rn ; H)
is defined (see [Tr]) by
F ⊗ Id : S0 (Rn ) ⊗ H → S0 (Rn ) ⊗ H
Definition
P∞ 3.3.2
P∞ Let Φ be a H–valued GRP (I) given by the expansion
n
Φ = f ij ⊗ H j
α i e , fij ∈ S −k (R ), i, j = 1, 2, . . . such that
P∞ P∞i=1 j=1 2 −pαj
i=1 j=1 kfij k−k (2N) < ∞ for some p ≥ 0. The Fourier transfor-
mation of Φ, denoted by F(Φ) is the unique H–valued GRP (I) defined by
the expansion
X ∞ X
∞
F(Φ) = fcij ⊗ Hαj ei . (3.79)
i=1 j=1
A simple application
Let A ∈ (S)−1 , δ ∈ S0 (R) be the Dirac delta distirbution, and consider
the SDE
d2
u(x, ω) − u(x, ω) = δ 0 (x) ⊗ A(ω), x ∈ R, ω ∈ Ω. (3.80)
dx2
1
u(x, ω) = sgn(x)e−|x| ⊗ A(ω).
2
P∞
Proof. Let the chaos expansion of A be A(ω) P∞ = j=1 aj Hαj (ω), ω ∈ Ω. We
seek for the solution in form of u(x, ω) = j=1 uj (x) ⊗ Hαj (ω). Now (3.80)
P ³ d2 ´ P
obtains the form ∞ j=1 dx 2 u j (x) − u j (x) ⊗Hαj (ω) = δ 0 (x)⊗ ∞ j=1 aj Hαj (ω),
from which we get the system of ODEs
d2
2
uj (x) − uj (x) = aj δ 0 (x), j ∈ N.
dx
Applying the Fourier transformation we get
1
−y 2 ûj (y) − ûj (y) = aj iy √ ,
2π
132 Applications to Singular SPDEs
i.e.
−iyaj
ûj (y) = √ .
2π(1 + y 2 )
Now we apply the inverse Fourier transformation to get
1
uj (x) = aj sgn(x)e−|x| , j ∈ N,
2
where sgn ∈ S0 (R) is understood in distributional sense. ¤
We can check that this solution is bona fide, also using the results from
Example 2.2.3. The function
½ aj −x
e , x>0
uj (x) = 2
aj x
−2e , x<0
and ½ aj −x
d2 e , x>0
uj (x) = 2
aj x + aj δ 0 (x) = uj (x) + aj δ 0 (x).
dx2 −2e , x<0
defines a GRP (I). This process was introduced in [Se], where it was shown
5
that it is an element of L(Sk (R), (S)−1 ) for k > 12 .
Consider now the one dimensional stochastic Helmholtz equation
d2
Φ(x, ω) + k 2 Φ(x, ω) = ∆(x, ω), x ∈ R, ω ∈ Ω, (3.82)
dx2
where k is a constant.
The (deterministic) Helmholtz equation is closely related to the wave
2 d2
equation: Considering the wave equation dtd 2 u(x, t) − c2 dx 2 u(x, t) = 0 and
3.3 Applications of the Fourier Transformation to GRPs (I) 133
applying a separation of variables u(x, t) = e−jvt X(x) one gets that X(x)
d2 2 v
satisfies ( dx 2 + k )X(x) = 0, k = c . Thus, the solution of the Helmholtz
we note that fbj ¶∈ exp S0 (R), since it has expansion fbj (y) =
First µ
P∞ − 1 t2 √
e 2 j 4π n
n=1 k2 −y 2
i hn (tj ) ξn (y), where hn are the Hermite polynomials, ξn
−t2j
√ P
are the Hermite functions, and (ke 2 −y2π)2 ∞ 2 2n −2k
n=1 |hn (tj )| (e ) < ∞ for k large
enough. Thus, we can apply the inverse Fourier transformation in exp S0 (R)
to obtain
1
fj (x) = sin(k(x − tj )) sgn(x − tj ),
2k
where sgn is understood again in distributional sense as sgn(x) = 2H(x) − 1,
and H is the Heaviside function.
Thus, the solution of the Helmholtz equation (3.82) is given by
∞
1 X
Φ(x, ω) = sin(k(x − tj )) sgn(x − tj ) ⊗ Hαj (ω), x ∈ R, ω ∈ Ω. (3.83)
2k j=1
It remains to prove that (3.83) is a well-defined GRP (I) i.e. the series
converges in L(exp S(R), (S)−1 ).
Denote by an,j = P hsin(k(x − tj )), ξn i the nth coefficient in the expan-
sion sin(k(x − tj )) = ∞ n=1 an,j ξn (x). By elementary calculus we get a1,j =
1 k2 √ 1 − k2
2 √ k2
−√ 4 πe
− 2 2π sin(ktj ), a2,j = √ 4 πe k 2π cos(ktj ), a3,j = 12 e− 2 (2k 2 −
√ k2
1) 4 π sin(ktj ), ..., an,j = Pn−1 (k)e− 2 f (ktj ), where Pn−1 is some polynomial
of order n − 1 and f is either the sine or the cosine function. By boundedness
of f we get
∞
X ∞
X
¡ ¢
2n −2l −k2
k sin(k(x − tj ))k2−l,exp = |an,j | 2
e =e P2n−2 (k)e−4nl < ∞
n=1 n=1
for l > 0 large enough. Moreover, k sin(k(x − tj ))k2−l,exp does not depend on
j ∈ N. Denote A = k sin(k(x − tj ))k2−l,exp . Thus,
∞
X ∞
X
j j
k sin(k(x − tj )) sgn(x − tj )k2−l,exp (2N)−pα =A (2N)−pα < ∞,
j=1 j=1
for p > 1. ¤
3.3 Applications of the Fourier Transformation to GRPs (I) 135
By considering only the zeroth terms in the expansion of the GRPs in-
volved in the Helmholtz equation we obtain following consequence of Propo-
sition 3.3.2.
Corollary 3.3.1 Let Φ be the solution of (3.82). Its generalized expectation
E(Φ) is the solution of the deterministic Helmholtz equation
d2
u(x) + k 2 u(x) = E(∆(x, ·)), x ∈ R,
dx2
where E(∆) is the generalized expectation of ∆.
Note that a similar consideration can be carried out for the Hilbert space
valued case, if we consider the H–valued GRP (I) defined in Example 2.2.1
(ii) as the right hand side of equation (3.82).
1 p
f3 (x) = √ sin( k1 (x − t3 )) sgn(x − t3 )
2 k1
k3 ³ p p p ´
− p 3 sin( k1 (x − t1 )) + k1 (t1 − x) cos( k1 (x − t1 )) sgn(x − t1 )
4 k1
k2 sgn(x − t2 ) ³p p p p ´
− √ k1 sin( k2 (x − t2 )) − k2 sin( k1 (x − t2 )) .
2(k1 − k2 ) k1 k2
White noise theory offers one of the most compelling instances of infinite
dimensional analysis, leading to a clear understanding of many empirical
phenomena. As a contribution to this noble theory, in Chapter 2 of the dis-
sertation, fundamental theorems were obtained characterizing the structure
of generalized stochastic processes. To begin the harvest of consequences, a
few examples have found place in Chapter 3 to illustrate the applications to
singular SPDEs. Further applications to SPDEs serving more concrete de-
mands, and the modeling of probabilistical properties of their solutions (e.g.
such as the martingale or the Markov property) are deferred but certainly
not denied. They remain as enticing possibilities for future investigations.
Bibliography
[AHR1] Albeverio S., Haba Z., Russo F., Trivial solutions for a non–linear
two–space dimensional wave equation perturbed by space–time white
noise, Stochastics Stochastics Rep. 56, No.1–2, pp. 127–160, 1996.
[AHR2] Albeverio S., Haba Z., Russo F., A two–space dimensional heat
equation perturbed by (Gaussian) white noise, Probab. Theory Re-
lat. Fields 121, pp. 319–366, 2001.
[DHPV] Delcroix A., Hasler M., Pilipović S., Valmorin V., Generalized func-
tion algebras as sequence space algebras, arXiv:math.FA/0206039
v1, 2002.
[Ev] Evans L.C., Partial Differential Equations, Grad. Stud. Math., Vol.
19, Amer. Math. Soc., Providence, 1998.
[FW] Floret K., Wloka J., Einführung in die Theorie der lokalkonvexen
Räume, Springer Verlag, 1968.
BIBLIOGRAPHY 139
[GS] Gel’fand I.M., Shilov G.E., Generalized functions, Vol. I, II, Aca-
demic Press, New York and London, 1964.
[HKPS] Hida T., Kuo H.H., Pothoff J., Streit L., White Noise - An Infinite
Dimensional Calculus, Kluwer Academic Publishers, 1993.
[HØUZ] Holden H., Øksendal B., Ubøe J., Zhang T., Stochastic Partial
Differential Equations - A Modeling, White Noise Functional Ap-
proach, Birkhäuser, 1996.
[In] Inaba H., Tapley B.T., Generalized Processes: A Theory and the
White Gaussian Process, SIAM J. Control, 13, pp 719–735, 1975.
[It] Itô K., Stationary Random Distributions, Mem. Coll. Sci. Kyoto
Univ., Ser. A, 28, pp 209–223, 1954.
[Ku] Kuo H.H., White Noise Distribution Theory, CRC Press, 1996.
[LØUZ] Lindstrøm T., Øksendal B., Ubøe J., Zhang T., Stability Properties
of Stochastic Partial Differential Equations, Stochastic Analysis and
Applications 13, pp 177–204, 1995.
140 BIBLIOGRAPHY
[NPR] Nedeljkov M., Pilipović S., Rajter-Ćirić D., Heat equation with sin-
gular potential and singular data, Proc. Roy. Soc. Edinburgh Sect.
A 135, No.4, pp 863–886, 2005.
[NR] Nedeljkov M., Rajter D., Nonlinear stochastic wave equation with
Colombeau generalized stochastic processes, Math. Models Methods
Appl. Sci. 12, No. 5, pp 665–688, 2002.
[Nu] Nualart D., Rozovski B., Weighted Stochastic Sobolev Spaces and
Bilinear SPDEs Driven by Space-Time White Noise, Journal of
Functional Analysis 149, pp 200–225, 1997.
BIBLIOGRAPHY 141
[Ru] Rudin W., Real and Complex Analysis, Tata McGraw-Hill, New
Delhi, 1977.
[Va] Våge G., Hilbert Space Methods Applied to Elliptic Stochastic Par-
tial Differential Equations, Stochastic analysis and related top-
ics, V(Silivri, 1994), pp 281–294, Progr.Probab., 38, Birkhäuser,
Boston, 1996.
B G
bilinear form, 10, 85, 112, 116 Gaussian
bounded, 85, 112 measure, 10, 20
coercive, 85, 117, 118 process, 21, 60, 63, 68–70
Brownian motion, 21, 22, 25, 26, 28, random variable, 60, 68, 70, 72
29, 57, 74 Gel’fand triple, 9, 12, 16, 20, 23
generalized
C expectation, 22, 86, 90, 106, 118,
chaos expansion, 21, 29, 44, 88, 96, 135
107, 108, 123, 131 function, 7, 10–20, 27–29, 31, 32,
Colombeau algebra, 7, 19–20, 29, 30, 38–40, 42, 44, 45, 57, 83, 95,
53–59, 83, 93, 95–106, 125– 96
128 random process, 27–136
random variable, 20–27, 29, 40,
D 44, 45, 86
Dirichlet problem, 81–128 solution, 87, 90, 92, 94, 104, 112,
distribution, 10–12, 16, 29, 30, 41, 74, 118
83, 96
Dirac delta, 17, 39, 41, 42, 73, 93, H
95, 129, 132, 133 Helmholtz equation, 132–136
point value, 18, 67
tempered, 11, 12, 14, 16, 74, 129
O
E operator
expectation, 62 adjoint, 11, 91, 110
convolution, 59, 93, 96
correlation, 59, 62–64
F differential, 13, 42, 49, 57, 102,
Fourier transformation, 129–136 109, 125
function elliptic, 82, 83, 90, 92, 95, 102
Fourier-Hermite, 21 Hilbert-Schmidt, 12
Heaviside, 17, 134 Laplace, 11, 16, 82, 86, 133
Hermite, 10–12, 26, 38, 74, 93, Wick–adjoint, 110, 111, 119
129 Wick-multiplication, 109
INDEX 145
P
product
tensor, 9–10, 12, 15, 18, 25, 84,
86, 87, 107, 123
Wick, 33, 55, 57, 58, 83, 87, 97,
135
S
space
countably Hilbert, 8, 9, 25
Hida, 7, 22
Hilbert, 8, 9, 11–13, 15, 18, 20,
23, 28, 29, 33, 44, 61, 75, 83–
86, 106, 110, 119, 128, 130,
131, 133
Kondratiev, 7, 21–24, 29, 33, 44,
83, 84, 86
nuclear, 8–12, 14, 15, 18, 23, 25,
28, 37, 40, 45, 47, 53, 55, 57,
59, 87, 98, 107
Schwartz, 7, 9, 11–12, 19, 25, 57,
130
Sobolev, 15–16, 19, 59, 65, 77,
82–84, 86, 87, 90, 93, 115,
123, 124
white noise, 20, 23
Zemanian, 7, 9, 13–15, 29, 33, 45,
53, 59, 86
W
white noise
coordinate process, 75
singular, 22, 25, 26, 29, 38, 39,
49, 57, 74, 93, 94
smoothed, 93, 94
Wick product, 44–49, 106–128
146
Curriculum Vitae
Personal Data
Name and surname: Dora Seleši (in Hungarian: Szöllössy Dóra)
Born on April 22nd 1978 in Zrenjanin, Serbia.
Education
2004 – Present Doctoral studies
2000 – 2004 Master studies (Mathematics, Major: analysis and probability)
1996 – 2000 Undergraduate studies (Mathematics, Major: graduate mathematician)
Interests
Language Skills
Hungarian, Serbian, English, German (Grosses Deutsches Sprachdiplom, Goethe Institut, 2005.)
Membership
MENSA Yugoslavia, 1998 – 2005.
Hobbys
Swimming, foreign languages, science fiction.
Work Experience
2004 – Present Assistant, Department of mathematics and informatics, Faculty of
Science, Novi Sad.
Assistant, Chair of Pharmacology, Faculty of Medicine, Novi Sad.
2000 – 2004 Junior Assistant, Department of mathematics and informatics, Faculty of
Science, Novi Sad.
147
Research
2006 – Present Engaged at the project Methods of functional analysis, ODEs and PDEs
with singularities, No. 144016 financed by the Ministry of science and
environment of Republic Serbia
2002 – 2006 Engaged at the project Functional analysis, ODEs and PDEs with
singularities, No. 1835 financed by the Ministry of science, technology and
development of Republic Serbia
Teaching
2000 – Present Conducting exercises for following courses:
Probability (for students of Mathematics and Informatics)
Statistics (for students of Mathematics and Informatics)
Statistical Modeling (for students of Mathematics)
Functional Analysis (for students of Informatics)
Applied Analysis (for students of Informatics)
Metric and normed spaces (for students of Mathematics)
Partial Differential Equations (for students of Mathematics)
Financial Mathematics II (for students of Mathematics)
Introduction to Calculus (for students of Mathematics)
Analysis I (for students of Mathematics)
Mathematics I (for students of Physics)
Mathematics II (for students of Physics)
Mathematics (for students of Pharmacology, also one group in English
language)
Scientific Impact
Attended several national and international conferences, seminars, summer and winter schools,
apart from those listed below.
Presentations at Conferences
D. Seleši, Stochastic Models of Financial Markets, XV Conference on Applied
Mathematics, Zlatibor, 2002.
D. Seleši, Expansion Theorems for Generalized Random Processes, Wick Products and
Applications to Stochastic Differential Equations, Generalized Functions 2004: Topics
in PDE, Harmonic Analysis and Mathematical Physics, Novi Sad, 2004.
Publication List
D. Seleši, Expansion of Generalized Stochastic Processes with Applications to
Equations, Master Thesis, University of Novi Sad, 2004.
S. Pilipović, D. Seleši, Structure Theorems for Generalized Random Processes, Acta
Mathematica Hungarica, 2007, to appear.
S. Pilipović, D. Seleši, Expansion Theorems for Generalized Random Processes, Wick
Products and Applications to Stochastic Differential Equations, Infinite Dimensional
Analysis, Quantum Probability and Related Topics, 2007, to appear in Vol.10.
D. Seleši, Hilbert Space Valued Generalized Random Processes, Novi Sad Journal of
Mathematics, submitted,
148
149
Biografija
Rođena sam 22. aprila 1978. godine u Zrenjaninu, gde sam završila
osnovnu školu i gimnaziju sa prosekom 5,00 i Vukovom diplomom,
kao i nižu muzičku školu. Od 1995. do 2001. bila sam stipendista
Republičke fondacije za razvoj naučnog i umetničkog podmlatka.
Studije matematike sam započela 1996. na Prirodno-matematičkom
fakultetu u Novom Sadu na smeru Diplomirani matematičar.
Diplomirala sam juna 2000. sa prosečnom ocenom 9,92 i nagradom
Mileva Marić-Einstein za najbolji seminarski rad.
Od oktobra 2000. godine sam zaposlena kao asistent-pripravnik, a od
juna 2004. kao asistent na Departmanu za matematiku i informatiku, PMF-a u Novom Sadu.
Tokom proteklih godina držala sam vežbe studentima matematike, informatike i fizike iz
predmeta: Verovatnoća, Statistika, Statističko modeliranje, Funkcionalna analiza, Metrički i
normirani prostori, Parcijalne diferencijalne jednačine, Finansijska matematika II, Uvod u
analizu, Analiza I, Matematika I, Matematika II. Od oktobra 2004. godine zaposlena sam kao
asistent i na katedri za farmaciju Medicinskog fakulteta, gde izvodim vežbe iz predmeta
Matematika, kao i vežbe iz predmeta Mathematics za studijsku grupu na engleskom jeziku.
Postdiplomske studije sam upisala novembra 2000. godine iz oblasti Analiza i verovatnoća. U
toku studija do 2002. položila sam sve ispite sa prosečnom ocenom 10,00. Magistarsku tezu pod
naslovom Ekspanzija uopštenih stohastičkih procesa sa primenama u jednačinama, odbranila
sam 16. februara 2004. na Prirodno-matematičkom fakultetu u Novom Sadu i time stekla zvanje
magistra matematičkih nauka. Od februara do jula 2005. godine boravila sam u Austriji kao
stipendista ÖAD na Institutu za tehničku matematiku, geometriju i informatiku (Institut für
Technische Mathematik, Geometrie und Bauinformatik), univerziteta u Innsbrucku.
U novembru 2005. položila sam sertifikacioni ispit Goethe Instituta Grosses Deutsches
Sprachdiplom za poznavanje nemačkog jezika na nivou maternjeg jezika sa najvišom ocenom.
Tečno govorim i engleski kao strani jezik, a maternji jezik mi je mađarski.
Oblasti mog naučnog interesovanja su stohastički procesi, teorija uopštenih funkcija,
diferencijalne jednačine, finansijska matematika i statistika. Imala sam izlaganja na
konferenciji PRIM XV (Zlatibor, 2002.) i na međunarodnoj konferenciji GF04 (Novi Sad,
2004.). Pored toga učestvovala sam još na više konferencija, seminara i letnjih/zimskih škola.
Autor ili koautor sam tri naučna rada, kao i nekoliko skripti za studente. Od 2002. angažovana
sam na projektu Metode funkcionalne analize, ODJ i PDJ sa singularitetima, koji finansira
Ministarstvo nauke i zaštite životne sredine Republike Srbije.
Oblasti interesovanja su mi pored nauke i naučna fantastika, učenje stranih jezika i plivanje. Od
1998. do 2005. godine bila sam član MENSA Jugoslavije.
Accession number:
ANO
Identification number:
INO
Document type: Monograph type
DT
Type of record: Printed text
TR
Contents code: PhD dissertation
CC
Author: Dora Seleši, MSc
AU
Mentor: Academic Stevan Pilipović, PhD
MN
UC:
Holding data: Library of the Department of Mathematics and Informatics,
Novi Sad
HD
Note:
N
Abstract: Subject of the dissertation are various classes of generalized
stochastic processes and their applications to solving singular stochastic
partial differential equations. Basically, the dissertation can be divided into
two parts. The first part (Chapter 2) is devoted to structural characteri-
zations of generalized random processes in terms of chaos expansions and
integral representations. The second part of the dissertation (Chapter 3)
involves applications of the obtained results to solving a stochastic Dirichlet
problem, where multiplication is modeled by the Wick product, and the
coefficients of the elliptic differential operator are Colombeau generalized
random processes.
AB
Redni broj:
RBR
Identifikacioni broj:
IBR
Tip dokumentacije: Monografska dokumentacija
TD
Tip zapisa: Tekstualni štampani materijal
TZ
Vrsta rada: Doktorska disertacija
VR
Autor: mr Dora Seleši
AU
Mentor: Akademik dr Stevan Pilipović
MN
UDK:
Čuva se: u biblioteci Departmana za matematiku i informatiku, Novi Sad
ČU
Važna napomena:
VN
Izvod: Doktorska disertacija je posvećena raznim klasama uopštenih
stohastičkih procesa i njihovim primenama na rešavanje singularnih sto-
hastičkih parcijalnih diferencijalnih jednačina. U osnovi, disertacija se
može podeliti na dva dela. Prvi deo disertacije (Glava 2) je posvećen
strukturnoj karakterizaciji uopštenih stohastičkih procesa u vidu haos
ekspanzije i integralne reprezentacije. Drugi deo disertacije (Glava 3) čini
primena dobijenih rezultata na rešavanje stohastičkog Dirihleovog problema
u kojem se množenje modelira Vikovim proizvodom, a koeficijenti eliptičnog
diferencijalnog operatora su Kolomboovi uopšteni stohastički procesi.
IZ
Članovi komisije:
Predsednik: dr Marko Nedeljkov, redovni profesor, Prirodno-matematički
fakultet, Univerzitet u Novom Sadu
Član: Akademik dr Stevan Pilipović, redovni profesor, Prirodno-matematički
fakultet, Univerzitet u Novom Sadu
Član: dr Michael Oberguggenberger, vanredni profesor, Fakultet za
grad̄evinu i arhitekturu, Univerzitet u Insbruku (Austrija)
Član: dr Zagorka Lozanov-Crvenković, redovni profesor, Prirodno-
matematički fakultet, Univerzitet u Novom Sadu
Član: dr Danijela Rajter-Ćirić, docent, Prirodno-matematički fakultet,
Univerzitet u Novom Sadu
KO