
The Laplacian $\Delta$ is the infinitesimal generator of isotropic Brownian
motion, being the limit process of normal diffusion, while the fractional
Laplacian $\Delta^{\beta/2}$ serves as the infinitesimal generator of the limit
process of isotropic L\'{e}vy process. Taking limit, in some sense, means that
the operators can approximate the physical process well after sufficient long
time. We introduce the nonlocal operators (being effective from the starting
time), which describe the general processes undergoing normal diffusion. For
anomalous diffusion, we extend to the anisotropic fractional Laplacian
$\Delta_m^{\beta/2}$ and the tempered one $\Delta_m^{\beta/2,\lambda}$ in
$\mathbb{R}^n$. Their definitions are proved to be equivalent to an alternative
one in Fourier space. Based on these new nonlocal diffusion operators, we
further derive the deterministic governing equations of some interesting
statistical observables of the very general jump processes with multiple
internal states. Finally, we consider the associated initial and boundary value
problems and prove their wellposedness of the Galerkin weak formulation in
$\mathbb{R}^n$. To obtain the coercivity, we claim that the probability density
function $m(Y)$ should be nondegenerate.

In this paper, we study the numerical schemes for the twodimensional
FokkerPlanck equation governing the probability density function of the
tempered fractional Brownian motion. The main challenges of the numerical
schemes come from the singularity in the time direction. When $0<H<0.5$, a
change of variables $\partial \left(t^{2H}\right)=2Ht^{2H1}\partial t$ avoids
the singularity of numerical computation at $t=0$, which naturally results in
nonuniform time discretization and greatly improves the computational
efficiency. For $0.5<H<1$, the time span dependent numerical scheme and
nonuniform time discretization are introduced to ensure the effectiveness of
the calculation and the computational efficiency. By numerically solving the
corresponding FokkerPlanck equation, we obtain the mean squared displacement
of stochastic processes, which conforms to the characteristics of the tempered
fractional Brownian motion.

Topological defects resulted from boundary constraints in confined liquid
crystals have attracted extensive research interests. In this paper, we use
numerical simulation to study the phase transition dynamics in the context of
stochastic resonance in a bistable liquid crystal device containing defects.
This device is made of nematic liquid crystals confined in a shallow square
well, and is described by the planar LebwohlLasher model. The stochastic phase
transition processes of the system in the presence of a weak oscillating
potential is simulated using an overdamped Langevin dynamics. Our simulation
results reveal that, depending on system size, the phase transition may follow
two distinct pathways: in small systems the preexisting defect structures at
the corners hold until the last stage and there is no newly formed defect point
in the bulk during the phase transition; In large systems new defect points
appear spontaneously in the bulk and eventually merge with the preexisting
defects at the corners. For both transition pathways stochastic resonance can
be observed, but show dramatic difference in their responses to the boundary
anchoring strength. In small systems we observe a "stickyboundary" effect for
a certain range of anchoring strength in which the phase transition gets stuck
and stochastic resonance becomes deactivated. Our work demonstrates the
dynamical interplay among defects, noises, and boundary conditions in confined
liquid crystals.

Applying proper orthogonal decomposition to a usual finite element (FE)
formulation for space fractional partial differential equation, we get a
reduced FE model, which greatly reduces the complexity of computation. Then,
the stability analysis and error estimate for the reduced model are presented.
Finally, we verify the effectiveness of the algorithm by numerical experiments.

Tempered fractional Laplacian is the generator of the tempered isotropic
L\'evy process [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale
Model. Simul., 16(1), 125149, 2018]. This paper provides the finite difference
discretization for the two dimensional tempered fractional Laplacian
$(\Delta+\lambda)^{\frac{\beta}{2}}$. Then we use it to solve the tempered
fractional Poisson equation with Dirichlet boundary conditions and derive the
error estimates. Numerical experiments verify the convergence rates and
effectiveness of the schemes.

Based on the theory of continuous time random walks (CTRW), we build the
models of characterizing the transitions among anomalous diffusions with
different diffusion exponents, often observed in natural world. In the CTRW
framework, we take the waiting time probability density function (PDF) as an
infinite series in three parameter MittagLeffler functions. According to the
models, the mean squared displacement of the process is analytically obtained
and numerically verified, in particular, the trend of its transition is shown;
furthermore the stochastic representation of the process is presented and the
positiveness of the PDF of the position of the particles is strictly proved.
Finally, the fractional moments of the model are calculated, and the analytical
solutions of the model with external harmonic potential are obtained and some
applications are proposed.

The generalized Langevin equation describes anomalous dynamics. Noise is not
only the origin of uncertainty but also plays a positive role in helping to
detect signal with information, termed stochastic resonance (SR). This paper
analyzes the anomalous resonant behaviors of the generalized Langevin system
with a multiplicative dichotomous noise and an internal tempered MittagLeffler
noise. For the system with fluctuating harmonic potential, we obtain the exact
expressions of several SR, such as, the first moment, the amplitude and the
autocorrelation function for the output signal as well as the signalnoise
ratio. We analyze the influence of the tempering parameter and memory exponent
on the bona fide SR and the general SR. Moreover, it is detected that the
critical memory exponent changes regularly with the increase of tempering
parameter. Almost all the theoretical results are validated by numerical
simulations.

Levy walk is a fundamental model with applications ranging from quantum
physics to paths of animal foraging. Taking animal foraging as an example, a
natural idea that comes to one's mind is to introduce the multiple internal
states for dealing with the dependence of the PDF of waiting time on the energy
of the animal and richness of the food at a particular location, etc; the
framework can also be used to model the moving trajectories of smart animals
without returning to the directions or locations which they come from
immediately. After building the Levy walk model with multiple internal states
and deriving the governing equation of the distribution of the positions of the
particles, some applications are discussed with specific transition matrices.
The type of diffusion for nonimmediatelyrepeating L\'{e}vy walk is uncovered,
and the distribution and average of first passage time are numerically
simulated.

The fractional Laplacian $\Delta^{\beta/2}$ is the generator of
$\beta$stable L\'evy process, which is the scaling limit of the L\'evy fight.
Due to the divergence of the second moment of the jump length of the L\'evy
fight it is not appropriate as a physical model in many practical applications.
However, using a parameter $\lambda$ to exponentially temper the isotropic
power law measure of the jump length leads to the tempered L\'evy fight, which
has finite second moment. For short time the tempered L\'evy fight exhibits the
dynamics of L\'evy fight while after sufficiently long time it turns to normal
diffusion. The generator of tempered $\beta$stable L\'evy process is the
tempered fractional Laplacian $(\Delta+\lambda)^{\beta/2}$ [W.H. Deng, B.Y. Li,
W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., in press, 2017]. In the
current work, we present new computational methods for the tempered fractional
Laplacian equation, including the cases with the homogeneous and nonhomogeneous
generalized Dirichlet type boundary conditions. We prove the wellposedness of
the Galerkin weak formulation and provide convergence analysis of the single
scaling Bspline and multiscale Riesz bases finite element methods. We propose
a technique for efficiently generating the entries of the dense stiffness
matrix and for solving the resulting algebraic equation by preconditioning. We
also present several numerical experiments to verify the theoretical results.

Aging, the process of growing old or maturing, is one of the most widely seen
natural phenomena in the world. For the stochastic processes, sometimes the
influence of aging can not be ignored. For example, in this paper, by analyzing
the functional distribution of the trajectories of aging particles performing
anomalous diffusion, we reveal that for the fraction of the occupation time
$T_+/t$ of strong aging particles, $\langle (T^+(t)^2)\rangle=\frac{1}{2}t^2$
with coefficient $\frac{1}{2}$, having no relation with the aging time $t_a$
and $\alpha$ and being completely different from the case of weak (none) aging.
In fact, we first build the models governing the corresponding functional
distributions, i.e., the aging forward and backward FeynmanKac equations; the
above result is one of the applications of the models. Another application of
the models is to solve the asymptotic behaviors of the distribution of the
first passage time, $g(t_a,t)$. The striking discovery is that for weakly aging
systems, $g(t_a,t)\sim t_a^{\frac{\alpha}{2}}t^{1\frac{\alpha}{2}}$, while
for strongly aging systems, $g(t_a,t)$ behaves as $ t_a^{\alpha1}t^{\alpha}$.

This paper focuses on the adaptive discontinuous Galerkin (DG) methods for
the tempered fractional (convection) diffusion equations. The DG schemes with
interior penalty for the diffusion term and numerical flux for the convection
term are used to solve the equations, and the detailed stability and
convergence analyses are provided. Based on the derived posteriori error
estimates, the local error indicator is designed. The theoretical results and
the effectiveness of the adaptive DG methods are respectively verified and
displayed by the extensive numerical experiments. The strategy of designing
adaptive schemes presented in this paper works for the general PDEs with
fractional operators.

In this paper we derive new uniform convergence estimates for the Vcycle MGM
applied to symmetric positive definite Toeplitz block tridiagonal matrices, by
also discussing few connections with previous results. More concretely, the
contributions of this paper are as follows: (1) It tackles the Toeplitz systems
directly for the elliptic PDEs. (2) Simple (traditional) restriction operator
and prolongation operator are employed in order to handle general Toeplitz
systems at each level of the recursion. Such a technique is then applied to
systems of algebraic equations generated by the difference scheme of the
twodimensional fractional FeynmanKac equation, which describes the joint
probability density function of nonBrownian motion. In particular, we consider
the two coarsening strategies, i.e., doubling the mesh size (geometric MGM) and
Galerkin approach (algebraic MGM), which lead to the distinct coarsening
stiffness matrices in the general case: however, several numerical experiments
show that the two algorithms produce almost the same error behaviour.

This paper provides a theoretical framework of deriving the forward and
backward FeynmanKac equations for the distribution of functionals of the path
of a particle undergoing both diffusion and chemical reaction. Very general
forms of the equations are obtained. Once given the diffusion type and reaction
rate, a specific forward or backward FeynmanKac equation can be obtained. The
listed in the paper include the ones for normal/anomalous diffusions and
reactions with linear/nonlinear rates. Using the derived equations, we also
study the occupation time in halfspace, the first passage time to a fixed
boundary, and the occupation time in halfspace with absorbing or reflecting
boundary conditions.

This paper further discusses the tempered fractional Brownian motion, its
ergodicity, and the derivation of the corresponding FokkerPlanck equation.
Then we introduce the generalized Langevin equation with the tempered
fractional Gaussian noise for a free particle, called tempered fractional
Langevin equation (tfLe). While the tempered fractional Brownian motion
displays localization diffusion for the long time limit and for the short time
its mean squared displacement has the asymptotic form $t^{2H}$, we show that
the asymptotic form of the mean squared displacement of the tfLe transits from
$t^2$ (ballistic diffusion for short time) to $t^{22H}$, and then to $t^2$
(again ballistic diffusion for long time). On the other hand, the overdamped
tfLe has the transition of the diffusion type from $t^{22H}$ to $t^2$
(ballistic diffusion). The tfLe with harmonic potential is also considered.

We provide and analyze the high order algorithms for the model describing the
functional distributions of particles performing anomalous motion with
powerlaw jump length and tempered powerlaw waiting time. The model is derived
in [Wu, Deng, and Barkai, Phys. Rev. E., 84 (2016), 032151], being called the
timetempered fractional FeynmanKac equation. The key step of designing the
algorithms is to discretize the time tempered fractional substantial
derivative, being defined as $${^S\!}D_t^{\gamma,\widetilde{\lambda}}
G(x,p,t)\!=\!D_t^{\gamma,\widetilde{\lambda}} G(x,p,t)\!\!\lambda^\gamma
G(x,p,t) ~{\rm with}~\widetilde{\lambda}=\lambda+ pU(x),\, p=\rho+J\eta,\,
J=\sqrt{1},$$ where $$D_t^{\gamma,\widetilde{\lambda}} G(x,p,t)
=\frac{1}{\Gamma(1\gamma)} \left[\frac{\partial}{\partial
t}+\widetilde{\lambda} \right]
\int_{0}^t{\left(tz\right)^{\gamma}}e^{\widetilde{\lambda}\cdot(tz)}{G(x,p,z)}dz,$$
and $\lambda \ge 0$, $0<\gamma<1$, $\rho>0$, and $\eta$ is a real number. The
designed schemes are unconditionally stable and have the global truncation
error $\mathcal{O}(\tau^2+h^2)$, being theoretically proved and numerically
verified in {\em complex} space. Moreover, some simulations for the
distributions of the first passage time are performed, and the second order
convergence is also obtained for solving the `physical' equation (without
artificial source term).

For the particles undergoing the anomalous diffusion with different waiting
time distributions for different internal states, we derive the FokkerPlanck
and FeymannKac equations, respectively, describing positions of the particles
and functional distributions of the trajectories of particles; in particular,
the equations governing the functional distribution of internal states are also
obtained. The dynamics of the stochastic processes are analyzed and the
applications, calculating the distribution of the first passage time and the
distribution of the fraction of the occupation time, of the equations are
given.

For characterizing the Brownian motion in a bounded domain: $\Omega$, it is
wellknown that the boundary conditions of the classical diffusion equation
just rely on the given information of the solution along the boundary of a
domain; on the contrary, for the L\'evy flights or tempered L\'evy flights in a
bounded domain, it involves the information of a solution in the complementary
set of $\Omega$, i.e., $\mathbb{R}^n\backslash \Omega$, with the potential
reason that paths of the corresponding stochastic process are discontinuous.
Guided by probability intuitions and the stochastic perspectives of anomalous
diffusion, we show the reasonable ways, ensuring the clear physical meaning and
wellposedness of the partial differential equations (PDEs), of specifying
`boundary' conditions for space fractional PDEs modeling the anomalous
diffusion. Some properties of the operators are discussed, and the
wellposednesses of the PDEs with generalized boundary conditions are proved.

The mean first exit (passage) time characterizes the average time of a
stochastic process never leaving a fixed region in the state space, while the
escape probability describes the likelihood of a transition from one region to
another for a stochastic system driven by discontinuous (with jumps) L\'evy
motion. This paper discusses the two deterministic quantities, mean first exit
time and escape probability, for the anomalous processes having the tempered
L\'{e}vy stable waiting times with the tempering index $\lambda>0$ and the
stability index $0<\alpha \le 1$; as for the distribution of jump lengths or
the type of the noises driving the system, two cases are considered, i.e.,
Gaussian white noise and nonGaussian (tempered) $\beta$stable ($0<\beta<2$)
L\'{e}vy noise. Firstly, we derive the nonlocal elliptic partial differential
equations (PDEs) governing the mean first exit time and escape probability.
Based on the derived PDEs, it is observed that the mean first exit time depends
strongly on the domain size and the values of $\alpha$, $\beta$ and $\lambda$;
when $\lambda$ is close to zero, the mean first exit time tends to $\infty$. In
particular, we also find an interesting result that the escape probability of a
particle with (tempered) powerlaw jumping length distribution has no relation
with the distribution of waiting times for the model considered in this paper.
For the solutions of the derived PDEs, the boundary layer phenomena are
observed, which inspires the motivation for developing the boundary layer
theory for nonlocal PDEs.

This paper focuses on providing the high order algorithms for the spacetime
tempered fractional diffusionwave equation. The designed schemes are
unconditionally stable and have the global truncation error
$\mathcal{O}(\tau^2+h^2)$, being theoretically proved and numerically verified.

For describing the probability distribution of the positions and times of
particles performing anomalous motion, fractional PDEs are derived from the
continuous time random walk models with waiting time distribution having
divergent first order moment and/or jump length distribution which has
divergent second order moment. It can be noted that the fractional PDEs are
essentially dealing with the multiscale issues. Generally the regularity of the
solutions for fractional PDEs is weak at the areas close to boundary and
initial time. This paper focuses on developing the applications of wavelet
bases to numerically solving fractional PDEs and digging out the potential
benefits of wavelet methods comparing with other numerical methods, especially
in the aspects of realizing preconditioning, adaptivity, and keeping the
Toeplitz structure. More specifically, the contributions of this paper are as
follows: 1. the techniques of efficiently generating stiffness matrix with
computational cost $\mathcal{O}(2^J)$ are provided for first, second, and any
order bases; 2. theoretically and numerically discuss the effective
preconditioner for timeindependent equation and multigrid method for
timedependent equation, respectively; 3. the wavelet adaptivity is detailedly
discussed and numerically applied to solving the timedependent (independent)
equations. In fact, having reliable, simple, and local regularity indicators is
the striking benefit of the wavelet in adaptively solving fractional PDEs (it
seems hard to give a local posteriori error estimate for the adaptive finite
element method because of the global property of the operator).

Recently, nonlocal models attract the wide interests of scientist. They
mainly come from two applied scientific fields: peridyanmics and anomalous
diffusion. Even though the matrices of the algebraic equation corresponding the
nonlocal models are usually Toeplitz (denote a0 as the principal diagonal
element, a1 as the trailing diagonal element, etc). There are still some
differences for the models in these two fields. For the model of anomalous
diffusion, a0/a1 is uniformly bounded; most of the time, a0/a1 of the model for
peridyanmics is unbounded as the stepsize h tends to zero. Based on the uniform
boundedness of a0/a1, the convergence of the twogrid method is well
established [Chan, Chang, and Sun, SIAM J. Sci. Comput., 19 (1998), pp.
516529; Pang and Sun, J. Comput. Phys., 231 (2012), pp. 693703; Chen, Wang,
Cheng, and Deng, BIT, 54 (2014), pp. 623647]. This paper provides the
detailed proof of the convergence of the twogrid method for the nonlocal model
of peridynamics. Some special cases of the full multigrid and the Vcycle
multigrid are also discussed. The numerical experiments are performed to verify
the convergence.

This paper focuses on providing the computation methods for the backward time
tempered fractional FeynmanKac equation, being one of the models recently
proposed in [Wu, Deng, and Barkai, Phys. Rev. E, 84 (2016) 032151]. The
discretization for the tempered fractional substantial derivative is derived,
and the corresponding finite difference and finite element schemes are designed
with well established stability and convergence. The performed numerical
experiments show the effectiveness of the presented schemes.

Because of the finiteness of the life span and boundedness of the physical
space, the more reasonable or physical choice is the tempered powerlaw instead
of pure powerlaw for the CTRW model in characterizing the waiting time and
jump length of the motion of particles. This paper focuses on providing the
variational formulation and efficient implementation for solving the
corresponding deterministic/macroscopic models, including the space tempered
fractional equation and time tempered fractional equation. The convergence,
numerical stability, and a series of variational equalities are theoretically
proved. And the theoretical results are confirmed by numerical experiments.

Fractional operators are widely used in mathematical models describing
abnormal and nonlocal phenomena. Although there are extensive numerical methods
for solving the corresponding model problems, theoretical analysis such as the
regularity result, or the relationship between the leftside and rightside
fractional operators are seldom mentioned. In stead of considering the
fractional derivative spaces, this paper starts from discussing the image
spaces of RiemannLiouville fractional integrals of $L_p(\Omega)$ functions,
since the fractional derivative operators that often used are all
pseudodifferential. Then high regularity situationthe image spaces of
RiemannLiouville fractional integral operators on $W^{m,p}(\Omega)$ space are
considered. Equivalent characterizations of the defined spaces, as well as of
the intersection of the leftside and rightside spaces are given. The behavior
of the functions in the defined spaces at both the nearby boundary point/ponits
and the points in the domain are demonstrated in a clear way. Besides, tempered
fractional operators show to be reciprocal to the corresponding
RiemannLiouville fractional operators, which is expected to make some efforts
on theoretical support for relevant numerical methods. Last, we also provide
some instructions on how to take advantage of the introduced spaces when
numerically solving fractional equations.

The traditional wave equation models wave propagation in an ideal conducting
medium. For characterizing the wave propagation in inhomogeneous media with
frequency dependent powerlaw attenuation, the spacetime fractional wave
equation appears; further incorporating the additive white Gaussian noise
coming from many natural sources leads to the stochastic spacetime fractional
wave equation. This paper discusses the Galerkin finite element approximations
for the stochastic spacetime fractional wave equation forced by an additive
spacetime white noise. We firstly discretize the spacetime additive noise,
which introduces a modeling error and results in a regularized stochastic
spacetime fractional wave equation; then the regularity of the regularized
equation is analyzed. For the discretization in space, the finite element
approximation is used and the definition of the discrete fractional Laplacian
is introduced. We derive the meansquared $L^2$norm priori estimates for the
modeling error and for the approximation error to the solution of the
regularized problem; and the numerical experiments are performed to confirm the
estimates. For the timestepping, we calculate the analytically obtained
MittagLeffler type function.