# Canonical Sample Spaces

for Random Dynamical Systems

###### Abstract.

This is an overview about natural sample spaces for differential equations driven by various noises. Appropriate sample spaces are needed in order to facilitate a random dynamical systems approach for stochastic differential equations. The noise could be white or colored, Gaussian or non-Gaussian, Markov or non-Markov, and semimartingale or non-semimartingale. Typical noises are defined in terms of Brownian motion, Lévy motion and fractional Brownian motion. In each of these cases, a canonical sample space with an appropriate metric (or topology that gives convergence concept) is introduced. Basic properties of canonical sample spaces, such as separability and completeness, are then discussed.

Moreover, a flow defined by shifts, is introduced on these canonical sample spaces. This flow has an invariant measure which is the probability distribution for Brownian motion, or Lévy motion or fractional Brownian motion. Thus canonical sample spaces are much richer in mathematical structures than the usual sample spaces in probability theory, as they have metric or topological structures, together with a shift flow (or driving flow) defined on it. This facilitates dynamical systems approaches for studying stochastic differential equations.

###### Key words and phrases:

Random dynamical systems, Brownian motion, fractional brownian motion, Wiener measure, Lévy motion, colored noise.This work was partially supported by NSF grants 0620539 and 0731201, and a Graduate Research Fellowship from Illinois Institute of Technology

###### 2000 Mathematics Subject Classification:

Primary: 37L55, 35R60; Secondary: 58B99, 35L20J. Duan]duan@iit.edu X. Kan]xkan@iit.edu Björn Schmalfuß]schmalfuss@upb.de

## 1. Random dynamical systems

Stochastic differential equations (SDEs) or stochastic partial differential equations (SPDEs) arise as mathematical models for complex systems under various random influences in engineering and science. Here we only consider random dynamical systems defined by SDEs. Such SDEs define random dynamical systems (RDS) with appropriate sample spaces, much as ordinary differential equations define deterministic dynamical systems.

Theory of random dynamical systems allows to discuss the qualitative behavior of stochastic systems that are not only driven by a white noise, Markov processes and semimartingales, but also driven by non-Markov processes or by non-semimartingales (e.g., fractional Brownian motion). To analyze these more general noise cases, appropriate sample spaces and ergodic theory play an important role. In this article, we discuss canonical or natural sample spaces for SDEs with various noises.

We recall the definition of a random dynamical system (RDS) in the state space , with the underlying probability space , and with time varying in or , as in Arnold [1]. The state space is equipped with the Euclidean norm (or length) and the usual scalar product .

Note that a deterministic dynamical system on the state space is a mapping , such that the flow property is satisfied:

for all and .

For a random dynamical system, we need an extra ingredient, namely, a model for the noise. Moreover, the flow property has a twist (thus called cocycle property) due to the effect of noise.

###### Definition 1.1.

(Random dynamical system)

A random dynamical system (RDS), denoted by
, consists of two ingredients:

(i) Model for the noise: A driving flow on a
probability space , i.e., a flow on the sample space , such that
is invariant, namely for all , and is measurable from to .

(ii) Model for the evolution: A cocycle over
, i.e. a measurable mapping , such that the family
of
random mappings satisfies the cocycle property:

(1) |

Here the driving flow describes stationary dynamics of noise in an appropriately chosen sample space (see below). The mathematical model for noises in engineering and science is usually a stationary generalized stochastic process [22].

When is continuous for all , we say that is a continuous RDS. Since the continuity in space is quite common for RDS generated by stochastic differential equations, we usually do not specifically mention this spatial continuity. We often call a continuous-time or discrete-time RDS when it is continuous or discrete in time .

It follows from [1] that , , is a homeomorphism of and

## 2. Dynamical systems driven by white noises

### 2.1. Brownian Motion

The physical phenomenon Brownian motion^{1}^{1}1The
phenomenon was first observed by Jan Ingenhouz in 1785, but was
subsequently rediscovered by Brown in 1828, according to sources
used by Eric Weisstein’s World of Physics, which can be
found on the Internet at
http://scienceworld.wolfram.com/physics/BrownianMotion.html is
due to the incessant hitting of pollen by the much smaller
molecules of the liquid. The hits occur a large number of times in
any small time interval, independently of each other and the
effect of a particular hit is small compared to the total effect
[43]. The physical theory of this motion, set up by Albert
Einstein in 1905, suggests
that the motion is random, and has the following properties:

i) the motion is continuous;

ii) it has independent increments;

iii) the increments are stationary and Gaussian random variables.

Figure 1 shows a sample path of the Brownian motion.

Intuitively speaking, property i) says that the sample path of the
Brownian motion is continuous. Property ii) means that the
displacements of a pollen particle over disjoint time intervals are
independent random variables. Property iii) is natural considering
the Central Limit Theorem.

We now describe the Brownian motion in the mathematical language,
i.e., introduce the first definition of the Brownian
motion[22, 2, 6].

Definition A: A stochastic process
defined on a probability space
is called a Brownian motion or a
Wiener process if the following conditions
hold:

1) a.s.;

2) the sample paths are a.s. continuous;

3) has stationary independent increments;

4) the increments has the normal
distribution with mean and variance , i.e.
for any .

Since the stochastic process is a mapping from a
probability space to a metric space and is governed by its law, i.e.
the probability measure on the metric space induced by
, we want to find it in order to have a better
understanding of the Brownian motion. In fact, we can find the
finite dimensional distribution of the Brownian motion using
condition 3) and 4) in Definition A, and this gives another
definition of the Brownian motion as we shall see. We first write
down the second definition[6] and then prove it is equivalent to Definition A.

Definition B: A stochastic process
defined on a probability space
is called a Brownian motion or a
Wiener process if the following conditions
hold:

1’) a.s.;

2’) the sample paths are a.s. continuous;

3’) for any finite sequence of times
and Borel sets

(2) |

where

(3) |

defined for any and is called the
transition density.

We can see the conditions 1’),2’) in Definition B are completely the
same as 1),2) in Definition A, respectively, the only thing we need
to do in the proof of the equivalence of the two is to show 3),4)
3’).

Proof 3),4) 3’)

Firstly, we show that Brownian motion has Markov property,
by proving that the conditional distribution of given
is the same as that given , in terms of
moment generating function [10]. In fact,

Secondly, we compute the joint distribution of the Brownian motion.

Finally, we have the finite dimensional distribution.

which is obviously the same as
3’).

Conversely, we need to show 3’) 3),4)

However, this part of work is completely done in
[6](see page 153-155).

### 2.2. Wiener Measure

In Section 1, we treat the Brownian motion as a stochastic
process, i.e., a collection of time-parameterized random
variables. Since a stochastic process is governed by its
law, i.e., the probability measure on the space
it maps to, one may ask the question why the Brownian motion can
be determined only by its finite dimensional distributions
although we have proved the so-defined Bownian motion (Definition
B) coincides with the definition describing the phenomenon
(Definition A). Motivated by this question, we shall treat the
Brownian motion as a “random variable”, which we call
random function [11], and this point of view
introduces the Wiener measure, the probability measure
defined on an appropriate space which gives us the Brownian
motion. We will sketch the ideas showing the existence and
uniqueness of this special probability measure, starting from
three different spaces, i.e, , the spaces of arbitrary functions, continuous
functions, and continuous functions passing zero at time , defined on ,
respectively.

First approach This way of showing the existence and
uniqueness of the Wiener measure is the one most often used to
construct a Markov Process and is rather technical. The main idea
is the following: First define a set function on the algebra
generated by the cylinder sets in according
to the finite-dimensional distribution of Brownian motion. Note
that it is a set function rather than a probability measure since
it is just finitely additive (not countably additive); then by the
celebrated Kolmogorov extension theorem [22], this
set function is uniquely extended to a probability measure on the
algebra generated by the algebra mentioned above.

More precisely, we now give the definitions and theorems.

denotes the set of all real-valued
functions
on .

Cylinder set [22] A subset of
of the form

(4) |

where is a Borel subset of is called a
cylinder set.

Fixing but varying over the entire
Borel subsets of , the class of such all cylinder
sets forms a algebra
.

Set function A set function
is defined on the measurable space
, given
by

By varying the choice of the finite time points , we get a class of such set functions , and this class is independent of the cylindrical expression of the functions in , i.e. if the cylinder set of (4) has another expression, say

then

Algebra and -algebra

denotes the algebra of subsets of
consisting all cylinder sets

denotes the smallest
algebra containing

is a finitely additive measure on
such that the
restriction of to
coincides with

###### Theorem 2.1 (Kolmogorov extension theorem).

The set function on is uniquely extendable to a probability measure on .

It has been proved [22, 32] that there exist a unique probability measure on , under which the coordinate mapping process

satisfies condition 3) and 4) of Definition A in section 2.1, so it does not introduce a “Brownian motion” as we defined. Fortunately, by another famous theorem of Kolmogorov[41], the continuity problem has been solved.

###### Theorem 2.2 (Kolmogorov continuity theorem).

Suppose that the process satisfies the following condition: for all there exist positive constants such that

then there exists a continuous modification of .

Modification [41] Suppose that are stochastic processes on , then we say that is a modification of if

Note that if is a modification of , then they have
the same finite-dimensional distributions.

Now there is one problem needs to be solved: Why is
a.s.?

Definition [11] Wiener
measure, denoted by , is a probability measure on a
measurable space
having the following two properties:

i) each is normally distributed under
with mean and variance , i.e,

For this is interpreted to mean that
.

ii) the stochastic process has independent
increments under , i.e., ,
are independent under .

Remark It may not seem so obvious that this approach in
fact proves the existence of the Wiener measure. Also, we cannot
obtain the Wiener measure simply by assigning measure one to
; see [32]. One might hope to
construct the measure directly on . Indeed, this is
the main idea of
the second approach we are going to present.

Second approach There are some advantages if we start
from since we can make it Polish (complete
and separable) by assigning an appropriate metric. It has
been proved that the Borel algebra generated by the open
sets in is equal to the algebra generated
by all the cylinder sets [32]. Thus there is a
totally
different way of constructing Wiener measure.

Metric on

(5) |

This metric introduces the topology of uniform convergence on
compact intervals. For convenience, we denote by
and its Borel -algebra by .

The main idea is to construct a sequence of probability measure
on such that
. Since (2.1) gives the f.d.d of
the Wiener measure , we must let the f.d.d of
weakly converges to those of . Although weak convergence
in need not follow the weak convergence of the f.d.d alone in
general, it does if we add the condition that is tight
[32, 11]. In fact, if a metric space is
separable and complete, then each probability measure on the
measurable space is tight (see Theorem 1.3 in [11]).
It can be shown that equipped with the metric defined by
(5), the canonical space is separable and complete
[11].
For detail, see [32][11].

Remark What is still not natural is that
a.s. So next we will talk about
constructing Wiener measure on . It appears that
the above two different approaches may be adapted to this case.

We can define the Brownian motion for as follows: Taking two independent Brownian motions and , we define

(6) |

We can work on the space for two-sided Brownian motion. The Wiener measure defined above should be similar defined in this space [1].

### 2.3. Canonical sample space

We consider a SDE

(7) |

The canonical sample space is , space of continuous functions that are zero at time zero, equipped with the compact open topology, the Borel field , and Wiener (probability) measure .

## 3. Dynamical systems driven by colored noises

Colored noise, or noise with non-zero correlation (’memory’) in time, are common in the physical, biological and engineering sciences [20]. A good candidate for modeling colored noise is the fractional Brownian motion.

### 3.1. Fractional Brownian motion

A fractional Brownian motion (fBM) , , with the Hurst parameter, is still a Gaussian process. But it is characterized by the stationarity of its increments and a memory property. The increments of the fractional Brownian motion are not independent, except in the standard Brownian motion case (). Thus it is not a Markov process except when . Specifically, a.s., mean , covariance , and variance . It also exhibits power scaling and path regularity properties with Hölder parameter , which are very distinct from Brownian motion. The standard Brownian motion is a special fBM with . Figure 2 is a sample path of the fractional Brownian motion with .

### 3.2. Canonical sample space

The stochastic calculus involving fBM is currently being developed; see e.g. [40, 54] and references therein. This will lead to more advances in the study of SDEs driven by colored fBM noise:

(8) |

Since the fBM is not Markov, the solution process is not Markov either. Thus the usual techniques from Markov processes will not be applicable to the study of SDEs driven by fBms. However, the random dynamical systems approach, as described in §1 above, looks promising [38, 19]. The theory of RDS, developed by Arnold and coworkers [1], describes the qualitative behavior of systems of stochastic differential equations in terms of stability, Lyapunov exponents, invariant manifolds, and attractors.

See Theorem 2.3 in [38, 9], the canonical sample space is , the set of continuous functions that is zero at zero, but the probability measure is generated by under the compact-open topology as defined in Section 2. The Borel field is .

We can introduce a flow on this canonical sample space defined by the shifts

In this case the measure is invariant, i.e.

for all .

## 4. Dynamical systems driven by non-Gaussian noises

In the last two sections, we considered dynamical systems driven by Gaussian noises (white or colored), in terms of Brownian motion or fractional Brownian motion. In this section, we discuss differential equation driven by non-Gaussian Lévy noises.

### 4.1. Lévy Motions

Gaussian processes like Brownian motion have been widely used to model fluctuations in engineering and science. For a particle in Brownian motion, its sample paths are continuous in time almost surely (i.e., no jumps), its mean square displacement increases linearly in time (i.e., normal diffusion), and its probability density function decays exponentially in space (i.e., light tail or exponential relaxation) [41]. But some complex phenomena involve non-Gaussian fluctuations, with peculiar properties such as anomalous diffusion (mean square displacement is a nonlinear power law of time) [12] and heavy tail (non-exponential relaxation) [57]. For instance, it has been argued that diffusion in a case of geophysical turbulence [52] is anomalous. Loosely speaking, the diffusion process consists of a series of “pauses”, when the particle is trapped by a coherent structure, and “flights” or “jumps” or other extreme events, when the particle moves in a jet flow. Moreover, anomalous electrical transport properties have been observed in some amorphous materials such as insulators, semiconductors and polymers, where transient current is asymptotically a power law function of time [49, 21]. Finally, some paleoclimatic data [15] indicates heavy tail distributions and some DNA data [52] shows long range power law decay for spatial correlation.

Lévy motions are thought to be appropriate models for non-Gaussian processes with jumps [48]. Let us recall that a Lévy motion , or , is a non-Gaussian process with independent and stationary increments, i.e., increments are stationary (therefore has no statistical dependence on ) and independent for any non overlapping time lags . Moreover, its sample paths are only continuous in probability, namely, as for any positive . With a suitable modification [4], these paths may be taken as càdlàg, i.e., paths are continuous on the right and have limits on the left. This continuity is weaker than the usual continuity in time.

This generalizes the Brownian motion or , as satisfies all these three conditions. But Additionally, (i) Almost every sample path of the Brownian motion is continuous in time in the usual sense and (ii) Brownian motion’s increments are Gaussian distributed.

Dynamical systems driven by non-Gaussian Lévy noises have attracted much attention recently [4, 30, 50]. Under certain conditions, the SDEs driven by Lévy motion generate stochastic flows [35, 4], and also generate random dynamical systems (or cocycles) in the sense of Arnold [1]. Recently, exit time estimates have been investigated by Imkeller and Pavlyukevich [27, 28] , and Yang and Duan [56] for SDEs driven by Lévy motion. This shows some qualitatively different dynamical behaviors between SDEs driven by Gaussian and non-Gaussian noises.

Lévy motions are named in honor of the
French probabilist Paul Lévy, who first studied them in 1930s.
From a mathematical point of view, there
are so many reasons why they are so important [4], such as:

There are many important examples, such as Brownian
motion, the Poisson process, stable processes, and subordinators.

They are generalizations of random walks to continuous
time.

They are the simplest class of processes whose paths
consist of continuous motion interspersed with jump
discontinuities
of random size appearing at random times.

Definition [4] A stochastic process
defined on a probability space is a Lévy motion if:

(L1) a.s.;

(L2) has independent and stationary increments;

(L3) is stochastically continuous, i.e. for
all and for all

With a suitable modification [4], these paths may be taken as càdlàg, i.e., paths are continuous on the right and have limits on the left. This continuity is weaker than the usual continuity in time.

Figure 3 is a sample path for a Lévy motion.

One way to understand the structure of the Lévy motions is to
employ Fourier analysis. It may be shown that each is
infinitely divisible. The infinitely divisible random
variables are characterized completely through their
characteristic functions by a beautiful formula, established by
Paul Lévy and A. Ya. Khintchine in the 1930s. The
related definitions and theorems are as follows.

Definition [4] The
characteristic function of a stochastic process
taking values in is the mapping defined by

where is the distribution of .

Definition [4] is
infinitely divisible if for each , there
exists a probability measure on with
characteristic function such that
, for each .

###### Theorem 4.1 (The Lévy-Khintchine Formula [4]).

If is a Lévy motion, then

where

(9) |

for some , a non-negative definite symmetric
matrix and a Borel measure , called Lévy
jump measure, on for which
.
Here is the indicator function of the set .

Conversely, given a mapping of the form (9), we can always
construct a Lévy motion for which and
call it a Lévy motion with the characteristics or generating
triple .

Remark: The so-defined Borel measure in (9)
is called the Lévy jump measure, which should not be
confused with the
probability measure induced by the Lévy motion, to be defined below.

The different terms which appear in the Lévy-Khintchine formula
have a probabilistic significance emphasized in [46].
Every Lévy motion is obtained as a sum of independent processes
with three types of characteristics , and
. Thus, the Lévy measure accounts for the jumps of
and the knowledge of permits to give a probabilistic
construction of ; see [46].

### 4.2. Canonical sample space

Consider a SDE driven by non-Gaussian Levy noise

(10) |

The canonical space has to be enlarged to include all the cadlag functions, i.e. functions that are right-continuous and have left limits, defined on and taking values in . This space is denoted as .

We adopt the same point of view as in Section 2
that a stochastic process is also a
random variable, i.e.

###### Remark 4.2.

However, it can be made complete and separable when endowed with the Skorohod metric [11]. With this special metric, we call a Skorohod space. The Skorohod metric on is defined as

where , with

and

where denotes the set of strictly increasing, continuous
functions from to itself.

The Borel field under this topology is denoted as
.
For studying the weak convergence
and tightness in , the same approach adopted in can
be applied except that the fact the natural projections are not
continuous need to be noticed [11].

Definition The probability measure, , in that makes every element in a sample Lévy path is called the Lévy probability measure. Note that this measure is not to be confused with the Lévy jump measure mentioned above.

## References

- [1] L. Arnold. Random Dynamical Systems. Springer, New York, 1998.
- [2] L. Arnold, Stochastic Differential Equations, John Wiley & Sons, New York, 1974.
- [3] L. Arnold and M. Scheutzow, Perfect cocycles through stochastic differential equations. Prob. Theory Rel. Fields 101 (1995), 65-88.
- [4] D. Applebaum, Lévy Processes and Stochastic Calculus. Cambridge University Press, Cambridge, UK, 2004.
- [5] R. B. Ash. Probability and Measure Theory. Second Edition, Academic Press, New York, 2000.
- [6] Z. Brzezniak, and T. Zastawniak. Basic Stochastic Processes. Springer, 1998.
- [7] M. Capinski and E. Kopp. Measure, Integral and Probability. Springer-Verlag, New York, 1999.
- [8] P. L. Chow, Stochastic Partial Differential Equations. Chapman & Hall/CRC, New York, 2007.
- [9] L. Decreusefond and D. Nualart, Flow properties of differential equations driven by fractional Brownian motion. In Stochastic differential equations: theory and applications, 249–262, Interdiscip. Math. Sci., 2, World Sci. Publ., Hackensack, NJ, 2007.
- [10] F. Klebaner. Introduction to stochastic Calculus with Application. Imperial College Press, Second Edition, 2005.
- [11] P. Billingsley. Convergence of Probability Measure. Wiley Series in probability and Statistics, Chicago, Second Edition, 1999.
- [12] J. P. Bouchaud and A. Georges, Anomalous diffusion in disordered media: Statistic mechanics, models and physical applications. Phys. Repts. 195 (1990), 127-293.
- [13] P. Boxler, Stochastische Zentrumsmannigfaltigkeiten. Ph.D. thesis, Institut fiir Dynamische Systeme, Universitat Bremen, 1988.
- [14] G. Da Prato and J. Zabczyk, Stochastic Equations in Infinite Dimensions, Cambridge University Press, 1992.
- [15] P. D. Ditlevsen, Observation of stable noise induced millennial climate changes from an ice record. Geophys. Res. Lett. 26 (1999), 1441-1444.
- [16] A. Du and J. Duan, Invariant manifold reduction for stochastic dynamical systems. Dynamical Systems and Applications 16(2007), 681-696.
- [17] R. M. Dudley, Real Analysis and Probability. Cambridge University Press, 2002.
- [18] T. Fujiwara, H. Kunita, Stochastic differenntial equations of jump type and Leévy flows in diffeomorphisms group,J. Math. Kyoto Univ. 25, 71-106, 1985.
- [19] M. J. Garrido-Atienza, B. Maslowski and B. Schmalfuss, Random attractors for ordinary differential equations driven by a fractional Brownian motion with Hurst parameter gretaer than 1/2. Preprint, 2008.
- [20] P. Hanggi and P. Jung, Colored Noise in Dynamical Systems. Advances in Chem. Phys., 89(1995),239-326.
- [21] M. P. Herrchen, Stochastic Modeling of Dispersive Diffusion by Non-Gaussian Noise. Doctorial Thesis, Swiss Federal Inst. of Tech., Zurich, 2001.
- [22] T. Hida. Brownian Motion. Springer, New York, 1974.
- [23] Z. Huang and J. Yan, Introduction to Infinite Dimensional Stochastic Analysis. Science Press/Kluwer Academic Pub., Beijing/New York, 1997.
- [24] P. Imkeller and C. Lederer, On the cohomology of flows of stochastic and random differential equations. Probab. Theory Related Fields 120 (2001), 209–235.
- [25] N. Ikeda and S. Watanabe. Stochastic Differential Equations and Diffusion Processes. North-Holland, New York, 1981.
- [26] P. Imkeller and A. Monahan (eds.), Conceptual Stochastic Climate Models. Special Issue:Stochastics and Dynamics, 2(2002),no.3.
- [27] P. Imkeller and I. Pavlyukevich, First exit time of SDEs driven by stable Lévy processes. Stoch. Proc. Appl. 116 (2006), 611-642.
- [28] P. Imkeller, I. Pavlyukevich and T. Wetzel, First exit times for Lévy-driven diffusions with exponentially light jumps. arXiv:0711.0982.
- [29] J. Jacod and A. N. Shiryaev, Limit Theorems for Stochastic Processes. Springer, New York, 1987.
- [30] A. Janicki and A. Weron, Simulation and Chaotic Behavior of Stable Stochastic Processes, Marcel Dekker, Inc., 1994.
- [31] O. Kallenberg, Foundations of Modern Probability. Second Edition, Applied probability Trust, 2002.
- [32] I. Karatzs, and Steven E. Shreve. Brownian Motion and Stochastic Calculus. Springer, New York, Second Edition, 1991.
- [33] P. E. Kloeden and E. Platen, Numerical solution of stochastic differential equations, Springer-Verlag, 1992; second corrected printing 1995.
- [34] H. Kunita. Stochastic flows and stochastic differential equations. Cambridge University Press, 1990.
- [35] H. Kunita, Stochastic differential equations based on Lévy processes and stochastic flows of diffeomorphisms. Real and stochastic analysis, 305–373, Trends Math., Birkhäuser Boston, Boston, MA, 2004.
- [36] H. Kunita, Stochastic differenntial equations with jumps and stochastic flows of diffeomorphisms ,(Itô’s stochastic calculus and probability theory.)Springer,Tokyo, 197-211,1996.
- [37] X. Liu, J. Duan, J. Liu and P. E. Kloeden, Synchronization of dissipative dynamical systems driven by non-Gaussian Lévy noises. Preprint, 2009. arXiv:0901.2446v1 [math.DS]
- [38] B. Maslowski and B. Schmalfuss, Random dynamical systems and stationary solutions of differential equations driven by the fractional Brownian motion. Stoch. Anal. Appl. 22 (2004), 1577- 1607.
- [39] Y. S. Mishura, Stochastic Calculus for Fractional Brownian Motion and Related Processes. Springer, New York, 2008.
- [40] D. Nualart, Stochastic calculus with respect to the fractional Brownian motion and applications. Contemporary Mathematics 336, 3-39, 2003.
- [41] B. Oksendal. Stochastic Differenntial Equations. Sixth Ed., Springer-Verlag, New York, 2003.
- [42] S. Peszat and J. Zabczyk, Stochastic Partial Differential Equations with Lévy Processes, Cambridge University Press, Cambridge, UK, 2007.
- [43] J. von Plato, Creating Modern Probability. Cambridge University Press, UK, 1994.
- [44] C. Prevot and M. Rockner, A Concise Course on Stochastic Partial Differential Equations, Lecture Notes in Mathematics, Vol. 1905. Springer, New York, 2007.
- [45] P. E. Protter Stochastic Integration and Differential Equations. Springer, New York, Second Edition, 2005.
- [46] D. Revuz, and M. Yor Continuous Martingales and Brownian Motion. Springer, New York, Third Edition, 2005.
- [47] B. L. Rozovskii, Stochastic Evolution Equations. Kluwer Academic Publishers, Boston, 1990.
- [ST] G. Samorodnitsky and M. Grigoriu. Stable Non-Gaussian Random Processes. Chapman & Hall, New York, 1994.
- [48] K.-I. Sato, Lévy Processes and Infinitely Divisible Distributions, Cambridge University Press, Cambridge, 1999.
- [49] H. Scher, M. F. Shlesinger and J. T. Bendler, Time-scale invariance in transport and relaxzation. Phys. Today p.26-34, Jan. 1991.
- [50] D. Schertzer, M. Larcheveque, J. Duan, V. Yanovsky and S. Lovejoy, Fractional Fokker–Planck Equation for Nonlinear Stochastic Differential Equations Driven by Non-Gaussian Levy Stable Noises. J. Math. Phys., 42 (2001), 200-212.
- [51] M. Scheutzow, On the perfection of crude cocycles. Random $ Computational Dynamics 4 (1996), 235-255.
- [52] M. F. Shlesinger, G. M. Zaslavsky and U. Frisch, Lévy Flights and Related Topics in Physics (Lecture Notes in Physics, 450. Springer-Verlag, Berlin, 1995).
- [53] R. Situ. Theory of Stochastic Differential Equations with Jumps and Applications Springer, 2005.
- [54] C. A. Tudor and F. Viens, Statistical aspects of the fractional stochastic calculus. Annals of Statistics, Vol. 35 (3) (2007), 1183-1212.
- [55] E. Waymire and J. Duan (Eds.). Probability and Partial Differential Equations in Modern Applied Mathematics. Springer-Verlag, 2005.
- [56] Z. Yang, and J. Duan An intermediate regime for exit phenomena driven by non-Gaussian Lévy noises. Stochastics and Dynamics 8 (2008), 583-591.
- [57] F. Yonezawa, Introduction to focused session on ‘anomalous relaxation’. J. Non-Cryst. Solids, 198-200 (1996), 503-506.