Modeling and simulation with operator scaling

Self-similar processes are useful in modeling diverse phenomena that exhibit scaling properties. Operator scaling allows a different scale factor in each coordinate. This paper develops practical methods for modeling and simulating stochastic processes with operator scaling. A simulation method for operator stable Levy processes is developed, based on a series representation, along with a Gaussian approximation of the small jumps. Several examples are given to illustrate practical applications. A classification of operator stable Levy processes in two dimensions is provided according to their exponents and symmetry groups. We conclude with some remarks and extensions to general operator self-similar processes.Comment: 29 pages, 13 figure

Inversions of Levy Measures and the Relation Between Long and Short Time Behavior of Levy Processes

The inversion of a Levy measure was first introduced (under a different name) in Sato 2007. We generalize the definition and give some properties. We then use inversions to derive a relationship between weak convergence of a Levy process to an infinite variance stable distribution when time approaches zero and weak convergence of a different Levy process as time approaches infinity. This allows us to get self contained conditions for a Levy process to converge to an infinite variance stable distribution as time approaches zero. We formulate our results both for general Levy processes and for the important class of tempered stable Levy processes. For this latter class, we give detailed results in terms of their Rosinski measures

Convolution-type derivatives, hitting-times of subordinators and time-changed C0C_0-semigroups

In this paper we will take under consideration subordinators and their inverse processes (hitting-times). We will present in general the governing equations of such processes by means of convolution-type integro-differential operators similar to the fractional derivatives. Furthermore we will discuss the concept of time-changed $C_0$-semigroup in case the time-change is performed by means of the hitting-time of a subordinator. We will show that such time-change give rise to bounded linear operators not preserving the semigroup property and we will present their governing equations by using again integro-differential operators. Such operators are non-local and therefore we will investigate the presence of long-range dependence.Comment: Final version, Potential analysis, 201

Numerical approximations for the tempered fractional Laplacian: Error analysis and applications

In this paper, we propose an accurate finite difference method to discretize the $d$-dimensional (for $d\ge 1$) tempered integral fractional Laplacian and apply it to study the tempered effects on the solution of problems arising in various applications. Compared to other existing methods, our method has higher accuracy and simpler implementation. Our numerical method has an accuracy of $O(h^\epsilon)$, for $u \in C^{0, \alpha+\epsilon} (\bar{\Omega})$ if $\alpha < 1$ (or $u \in C^{1, \alpha-1+\epsilon} (\bar{\Omega})$ if $\alpha \ge 1$) with $\epsilon > 0$, suggesting the minimum consistency conditions. The accuracy can be improved to $O(h^2)$, for $u \in C^{2, \alpha+\epsilon} (\bar{\Omega})$ if $\alpha < 1$ (or $u \in C^{3, \alpha - 1 + \epsilon} (\bar{\Omega})$ if $\alpha \ge 1$). Numerical experiments confirm our analytical results and provide insights in solving the tempered fractional Poisson problem. It suggests that to achieve the second order of accuracy, our method only requires the solution $u \in C^{1,1}(\bar{\Omega})$ for any $0<\alpha<2$. Moreover, if the solution of tempered fractional Poisson problems satisfies $u \in C^{p, s}(\bar{\Omega})$ for $p = 0, 1$ and $0<s \le 1$, our method has the accuracy of $O(h^{p+s})$. Since our method yields a (multilevel) Toeplitz stiffness matrix, one can design fast algorithms via the fast Fourier transform for efficient simulations. Finally, we apply it together with fast algorithms to study the tempered effects on the solutions of various tempered fractional PDEs, including the Allen-Cahn equation and Gray-Scott equations.Comment: 21 pages, 11 figures, 3 table