Source author record

Mark M. Meerschaert

Mark M. Meerschaert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2016arXiv

Domain and range symmetries of operator fractional Brownian fields

An operator fractional Brownian field (OFBF) is a Gaussian, stationary increment R^n-valued random field on R^m that satisfies the operator self-similarity property {X(c^E t)}_{t in R^m} L= {c^H X(t)}_{t in R^m}, c > 0, for two matrix exponents (E,H). In this paper, we characterize the domain and range symmetries of OFBF, respectively, as maximal groups with respect to equivalence classes generated by orbits and, based on a new anisotropic polar-harmonizable representation of OFBF, as intersections of centralizers. We also describe the sets of possible pairs of domain and range symmetry groups in dimensions (m,1) and (2,2).

preprint2016arXiv

Exponents of operator self-similar random fields

If X(c^E t) and c^H X(t) have the same finite-dimensional distributions for some linear operators E and H, we say that the random vector field X(t) is operator self-similar. The exponents E and H are not unique in general, due to symmetry. This paper characterizes the possible set of range exponents H for a given domain exponent, and conversely, the set of domain exponents E for a given range exponent.

preprint2016arXiv

Reflected Spectrally Negative Stable Processes and their Governing Equations

This paper explicitly computes the transition densities of a spectrally negative stable process with index greater than one, reflected at its infimum. First we derive the forward equation using the theory of sun-dual semigroups. The resulting forward equation is a boundary value problem on the positive half-line that involves a negative Riemann-Liouville fractional derivative in space, and a fractional reflecting boundary condition at the origin. Then we apply numerical methods to explicitly compute the transition density of this space-inhomogeneous Markov process, for any starting point, to any desired degree of accuracy. Finally, we discuss an application to fractional Cauchy problems, which involve a positive Caputo fractional derivative in time.

preprint2016arXiv

Space-time fractional Dirichlet problems

This paper establishes explicit solutions for fractional diffusion problems on bounded domains. It also gives stochastic solutions, in terms of Markov processes time-changed by an inverse stable subordinator whose index equals the order of the fractional time derivative. Some applications are given, to demonstrate how to specify a well-posed Dirichlet problem for space-time fractional diffusions in one or several variables. This solves an open problem in numerical analysis.

preprint2014arXiv

Extreme value statistics for truncated Pareto-type distributions

Recently attention has been drawn to practical problems with the use of unbounded Pareto distributions, for instance when there are natural upper bounds that truncate the probability tail. Aban, Meerschaert and Panorska (2006) derived the maximum likelihood estimator for the Pareto tail index of a truncated Pareto distribution with a right truncation point $T$. The Hill (1975) estimator is then obtained by letting $T \to \infty$. The problem of extreme value estimation under right truncation was also introduced in Nuyts (2010) who proposed a similar estimator for the tail index and considered trimming of the number of extreme order statistics. Given that in practice one does not always know whether the distribution is truncated or not, we discuss estimators for the Pareto index and extreme quantiles both under truncated and non-truncated Pareto-type distributions. We also propose a truncated Pareto QQ-plot in order to help deciding between a truncated and a non-truncated case. In this way we extend the classical extreme value methodology adding the truncated Pareto-type model with truncation point $T \to \infty$ as the sample size $n \to \infty$. Finally we present some practical examples, asymptotics and simulation results.

preprint2014arXiv

Semi-Markov approach to continuous time random walk limit processes

Continuous time random walks (CTRWs) are versatile models for anomalous diffusion processes that have found widespread application in the quantitative sciences. Their scaling limits are typically non-Markovian, and the computation of their finite-dimensional distributions is an important open problem. This paper develops a general semi-Markov theory for CTRW limit processes in $\mathbb{R}^d$ with infinitely many particle jumps (renewals) in finite time intervals. The particle jumps and waiting times can be coupled and vary with space and time. By augmenting the state space to include the scaling limits of renewal times, a CTRW limit process can be embedded in a Markov process. Explicit analytic expressions for the transition kernels of these Markov processes are then derived, which allow the computation of all finite dimensional distributions for CTRW limits. Two examples illustrate the proposed method.

preprint2011arXiv

Fractal Dimensions for Continuous Time Random Walk Limits

In a continuous time random walk (CTRW), each random jump follows a random waiting time. CTRW scaling limits are time-changed processes that model anomalous diffusion. The outer process describes particle jumps, and the non-Markovian inner process (or time change) accounts for waiting times between jumps. This paper studies fractal properties of the sample functions of a time-changed process, and establishes some general results on the Hausdorff and packing dimensions of its range and graph. Then those results are applied to CTRW scaling limits.

preprint2011arXiv

Tempered stable laws as random walk limits

Stable laws can be tempered by modifying the Lévy measure to cool the probability of large jumps. Tempered stable laws retain their signature power law behavior at infinity, and infinite divisibility. This paper develops random walk models that converge to a tempered stable law under a triangular array scheme. Since tempered stable laws and processes are useful in statistical physics, these random walk models can provide a basic physical model for the underlying physical phenomena.

preprint2011arXiv

The fractional Poisson process and the inverse stable subordinator

The fractional Poisson process is a renewal process with Mittag-Leffler waiting times. Its distributions solve a time-fractional analogue of the Kolmogorov forward equation for a Poisson process. This paper shows that a traditional Poisson process, with the time variable replaced by an independent inverse stable subordinator, is also a fractional Poisson process. This result unifies the two main approaches in the stochastic theory of time-fractional diffusion equations. The equivalence extends to a broad class of renewal processes that include models for tempered fractional diffusion, and distributed-order (e.g., ultraslow) fractional diffusion. The paper also establishes an interesting connection between the fractional Poisson process and Brownian time.

preprint2011arXiv

Transient anomalous sub-diffusion on bounded domains

This paper develops strong solutions and stochastic solutions for the tempered fractional diffusion equation on bounded domains. First the eigenvalue problem for tempered fractional derivatives is solved. Then a separation of variables, and eigenfunction expansions in time and space, are used to write strong solutions. Finally, stochastic solutions are written in terms of an inverse subordinator.

preprint2009arXiv

Distributed-order fractional Cauchy problems on bounded domains

In a fractional Cauchy problem, the usual first order time derivative is replaced by a fractional derivative. The fractional derivative models time delays in a diffusion process. The order of the fractional derivative can be distributed over the unit interval, to model a mixture of delay sources. In this paper, we provide explicit strong solutions and stochastic analogues for distributed-order fractional Cauchy problems on bounded domains with Dirichlet boundary conditions. Stochastic solutions are constructed using a non-Markovian time change of a killed Markov process generated by a uniformly elliptic second order space derivative operator.

preprint2006arXiv

Operator Scaling Stable Random Fields

A scalar valued random field is called operator-scaling if it satisfies a self-similarity property for some matrix E with positive real parts of the eigenvalues. We present a moving average and a harmonizable representation of stable operator scaling random fields by utilizing so called E-homogeneous functions. These fields also have stationary increments and are stochastically continuous. In the Gaussian case critical Hölder-exponents and the Hausdorff-dimension of the sample paths are also obtained.