Source author record

Federico Polito

Federico Polito 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

36works
13topics
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

36 published item(s)

preprint2023arXiv

Semi-Markovian discrete-time telegraph process with generalized Sibuya waiting times

In a recent work we introduced a semi-Markovian discrete-time generalization of the telegraph process. We referred this random walk to as squirrel random walk (SRW). The SRW is a discrete-time random walk on the one-dimensional infinite lattice where the step direction is reversed at arrival times of a discrete-time renewal process and remains unchanged at uneventful time instants. We first recall general notions of the SRW. The main subject of the paper is the study of the SRW where the step direction switches at the arrival times of a generalization of the Sibuya discrete-time renewal process (GSP) which only recently appeared in the literature. The waiting time density of the GSP, the `generalized Sibuya distribution' (GSD) is such that the moments are finite up to a certain order $r\leq m-1$ ($m \geq 1$) and diverging for orders $r \geq m$ capturing all behaviors from broad to narrow and containing the standard Sibuya distribution as a special case ($m=1$). We also derive some new representations for the generating functions related to the GSD. We show that the generalized Sibuya SRW exhibits several regimes of anomalous diffusion depending on the lowest order $m$ of diverging GSD moment. The generalized Sibuya SRW opens various new directions in anomalous physics.

preprint2022arXiv

Squirrels can remember little: A random walk with jump reversals induced by a discrete-time renewal process

We consider a class of discrete-time random walks with directed unit steps on the integer line. The direction of the steps is reversed at the time instants of events in a discrete-time renewal process and is maintained at uneventful time instants. This model represents a discrete-time semi-Markovian generalization of the telegraph process. We derive exact formulae for the propagator using generating functions. We prove that for geometrically distributed waiting times in the diffusive limit, this walk converges to the classical telegraph process. We consider the large-time asymptotics of the expected position: For waiting time densities with finite mean the walker remains in the average localized close to the departure site whereas escapes for fat-tailed waiting-time densities (i.e. densities with infinite mean) by a sublinear power-law. We explore anomalous diffusion features by accounting for the `aging effect' as a hallmark of non-Markovianity where the discrete-time version of the `aging renewal process' comes into play. By deriving pertinent distributions of this process we obtain explicit formulae for the variance when the waiting-times are Sibuya-distributed. In this case and generally for fat-tailed waiting time PDFs a $t^2$-ballistic superdiffusive scaling emerges in the large time limit. In contrast if the waiting time PDF between the step reversals is light-tailed (`narrow' with finite mean and variance) the walk exhibits normal diffusion and for `broad' waiting time PDFs (with finite mean and infinite variance) superdiffusive large time scaling. We also consider time-changed versions where the walk is subordinated to a continuous-time point process such as the time-fractional Poisson process. This defines a new class of biased continuous-time random walks exhibiting several regimes of anomalous diffusion.

preprint2021arXiv

A fractional generalization of the Dirichlet distribution and related distributions

This paper is devoted to a fractional generalization of the Dirichlet distribution. The form of the multivariate distribution is derived assuming that the $n$ partitions of the interval $[0,W_n]$ are independent and identically distributed random variables following the generalized Mittag-Leffler distribution. The expected value and variance of the one-dimensional marginal are derived as well as the form of its probability density function. A related generalized Dirichlet distribution is studied that provides a reasonable approximation for some values of the parameters. The relation between this distribution and other generalizations of the Dirichlet distribution is discussed. Monte Carlo simulations of the one-dimensional marginals for both distributions are presented.

preprint2021arXiv

Prabhakar discrete-time generalization of the time-fractional Poisson process and related random walks

In recent years a huge interdisciplinary field has emerged which is devoted to the complex dynamics of anomalous transport with long-time memory and non-markovian features. It was found that the framework of fractional calculus and its generalizations are able to capture these phenomena. Many of the classical models are based on continuous-time renewal processes and use the Montroll Weiss continuous time random walk (CTRW) approach. On the other hand their discrete time counterparts are rarely considered in the literature despite their importance in various applications. The goal of the present paper is to give a brief sketch of our recently introduced discrete-time Prabhakar generalization of the fractional Poisson process and the related discrete-time random walk (DTRW) model. We show that this counting process is connected with the continuous time Prabhakar renewal process by a (well scaled) continuous-time limit. We deduce the state probabilities and discrete time generalized fractional Kolmogorov-Feller equations governing the Prabhakar DTRW and discuss effects such as long time memory (nonmarkovianity) as a hallmark of the complexity of the process.

preprint2020arXiv

A practical guide to Prabhakar fractional calculus

The Mittag-Leffler function is universally acclaimed as the Queen function of fractional calculus. The aim of this work is to survey the key results and applications emerging from the three-parameter generalization of this function, known as the Prabhakar function. Specifically, after reviewing key historical events that led to the discovery and modern development of this peculiar function, we discuss how the latter allows one to introduce an enhanced scheme for fractional calculus. Then, we summarize the progress in the application of this new general framework to physics and renewal processes. We also provide a collection of results on the numerical evaluation of the Prabhakar function.

preprint2020arXiv

Flexible models for overdispersed and underdispersed count data

Within the framework of probability models for overdispersed count data, we propose the generalized fractional Poisson distribution (gfPd), which is a natural generalization of the fractional Poisson distribution (fPd), and the standard Poisson distribution. We derive some properties of gfPd and more specifically we study moments, limiting behavior and other features of fPd. The skewness suggests that fPd can be left-skewed, right-skewed or symmetric; this makes the model flexible and appealing in practice. We apply the model to real big count data and estimate the model parameters using maximum likelihood. Then, we turn to the very general class of weighted Poisson distributions (WPD's) to allow both overdispersion and underdispersion. Similarly to Kemp's generalized hypergeometric probability distribution, which is based on hypergeometric functions, we analyze a class of WPD's related to a generalization of Mittag--Leffler functions. The proposed class of distributions includes the well-known COM-Poisson and the hyper-Poisson models. We characterize conditions on the parameters allowing for overdispersion and underdispersion, and analyze two special cases of interest which have not yet appeared in the literature.

preprint2020arXiv

On discrete-time semi-Markov processes

In the last years, many authors studied a class of continuous time semi-Markov processes obtained by time-changing Markov processes by hitting times of independent subordinators. Such processes are governed by integro-differential convolution equations of generalized fractional type. The aim of this paper is to develop the discrete-time version of such a theory. We show that a class of discrete-time semi-Markov chains can be seen as time-changed Markov chains and we obtain governing convolution type equations. Such processes converge weakly to those in continuous time under suitable scaling limits.

preprint2020arXiv

On the continuous-time limit of the Barabási-Albert random graph

We prove that the Barabási-Albert model converges weakly to a set of generalized Yule models via an appropriate scaling. To pursue this aim we superimpose to its graph structure a suitable set of processes that we call the planted model and we introduce an ad-hoc sampling procedure. The use of the obtained limit process represents an alternative and advantageous way of looking at some of the asymptotic properties of the Barabási-Albert random graph.

preprint2020arXiv

Studies on generalized Yule models

We present a generalization of the Yule model for macroevolution in which, for the appearance of genera, we consider point processes with the order statistics property, while for the growth of species we use nonlinear time-fractional pure birth processes or a critical birth-death process. Further, in specific cases we derive the explicit form of the distribution of the number of species of a genus chosen uniformly at random for each time. Besides, we introduce a time-changed mixed Poisson process with the same marginal distribution as that of the time-fractional Poisson process.

preprint2017arXiv

A note on Hadamard fractional differential equations with varying coefficients and their applications in probability

In this paper we show several connections between special functions arising from generalized COM-Poisson-type statistical distributions and integro-differential equations with varying coefficients involving Hadamard-type operators. New analytical results are obtained, showing the particular role of Hadamard-type derivatives in connection with a recently introduced generalization of the Le Roy function. We are also able to prove a general connection between fractional hyper-Bessel-type equations involving Hadamard operators and Le Roy functions.

preprint2017arXiv

On a class of Time-fractional Continuous-state Branching Processes

We propose a class of non-Markov population models with continuous or discrete state space via a limiting procedure involving sequences of rescaled and randomly time-changed Galton--Watson processes. The class includes as specific cases the classical continuous-state branching processes and Markov branching processes. Several results such as the expressions of moments and the branching inequality governing the evolution of the process are presented and commented. The generalized Feller branching diffusion and the fractional Yule process are analyzed in detail as special cases of the general model.

preprint2016arXiv

A Generalization of the Space-Fractional Poisson Process and its Connection to some Lévy Processes

This paper introduces a generalization of the so-called space-fractional Poisson process by extending the difference operator acting on state space present in the associated difference-differential equations to a much more general form. It turns out that this generalization can be put in relation to a specific subordination of a homogeneous Poisson process by means of a subordinator for which it is possible to express the characterizing Lévy measure explicitly. Moreover, the law of this subordinator solves a one-sided first-order differential equation in which a particular convolution-type integral operator appears, called Prabhakar derivative. In the last section of the paper, a similar model is introduced in which the Prabhakar derivative also acts in time. In this case, too, the probability generating function of the corresponding process and the probability distribution are determined.

preprint2016arXiv

Random Graphs Associated to some Discrete and Continuous Time Preferential Attachment Models

We give a common description of Simon, Barabási--Albert, II-PA and Price growth models, by introducing suitable random graph processes with preferential attachment mechanisms. Through the II-PA model, we prove the conditions for which the asymptotic degree distribution of the Barabási--Albert model coincides with the asymptotic in-degree distribution of the Simon model. Furthermore, we show that when the number of vertices in the Simon model (with parameter $α$) goes to infinity, a portion of them behave as a Yule model with parameters $(λ,β) = (1-α,1)$, and through this relation we explain why asymptotic properties of a random vertex in Simon model, coincide with the asymptotic properties of a random genus in Yule model. As a by-product of our analysis, we prove the explicit expression of the in-degree distribution for the II-PA model, given without proof in \cite{Newman2005}. References to traditional and recent applications of the these models are also discussed.

preprint2015arXiv

Fractional diffusions with time-varying coefficients

This paper is concerned with the fractionalized diffusion equations governing the law of the fractional Brownian motion $B_H(t)$. We obtain solutions of these equations which are probability laws extending that of $B_H(t)$. Our analysis is based on McBride fractional operators generalizing the hyper-Bessel operators $L$ and converting their fractional power $L^α$ into Erdélyi--Kober fractional integrals. We study also probabilistic properties of the r.v.'s whose distributions satisfy space-time fractional equations involving Caputo and Riesz fractional derivatives. Some results emerging from the analysis of fractional equations with time-varying coefficients have the form of distributions of time-changed r.v.'s.

preprint2014arXiv

Discussion on the paper "On Simulation and Properties of the Stable Law" by L. Devroye and L. James

We congratulate the authors for the interesting paper. The reading has been really pleasant and instructive. We discuss briefly only some of the interesting results given in Devroye and James "On simulation and properties of the stable law", 2014 with particular attention to evolution problems. The contribution of the results collected in the paper is useful in a more wide class of applications in many areas of applied mathematics.

preprint2014arXiv

Fractional Klein-Gordon equation for linear dispersive phenomena: analytical methods and applications

In this paper we discuss some exact results related to the fractional Klein--Gordon equation involving fractional powers of the D'Alembert operator. By means of a space-time transformation, we reduce the fractional Klein--Gordon equation to a fractional hyper-Bessel-type equation. We find an exact analytic solution by using the McBride theory of fractional powers of hyper-Bessel operators. A discussion of these results within the framework of linear dispersive wave equations is provided. We also present exact solutions of the fractional Klein-Gordon equation in the higher dimensional cases. Finally, we suggest a method of finding travelling wave solutions of the nonlinear fractional Klein-Gordon equation with power law nonlinearities.

preprint2014arXiv

Fractional Klein-Gordon equations and related stochastic processes

This paper presents finite-velocity random motions driven by fractional Klein-Gordon equations of order $α\in (0,1]$. A key tool in the analysis is played by the McBride's theory which converts fractional hyper-Bessel operators into Erdelyi-Kober integral operators. Special attention is payed to the fractional telegraph process whose space-dependent distribution solves a non-homogeneous fractional Klein-Gordon equation. The distribution of the fractional telegraph process for $α= 1$ coincides with that of the classical telegraph process and its driving equation converts into the homogeneous Klein-Gordon equation. Fractional planar random motions at finite velocity are also investigated, the corresponding distributions obtained as well as the explicit form of the governing equations. Fractionality is reflected into the underlying random motion because in each time interval a binomial number of deviations $B(n,α)$ (with uniformly-distributed orientation) are considered. The parameter $n$ of $B(n,α)$ is itself a random variable with fractional Poisson distribution, so that fractionality acts as a subsampling of the changes of directions. Finally the behaviour of each coordinate of the planar motion is examined and the corresponding densities obtained. Extensions to $N$-dimensional fractional random flights are envisaged as well as the fractional counterpart of the Euler-Poisson-Darboux equation to which our theory applies.

preprint2014arXiv

Hilfer-Prabhakar Derivatives and Some Applications

We present a generalization of Hilfer derivatives in which Riemann--Liouville integrals are replaced by more general Prabhakar integrals. We analyze and discuss its properties. Further, we show some applications of these generalized Hilfer-Prabhakar derivatives in classical equations of mathematical physics, like the heat and the free electron laser equations, and in difference-differential equations governing the dynamics of generalized renewal stochastic processes.

preprint2014arXiv

On Firing Rate Estimation for Dependent Interspike Intervals

If interspike intervals are dependent the instantaneous firing rate does not catch important features of spike trains. In this case the conditional instantaneous rate plays the role of the instantaneous firing rate for the case of samples of independent interspike intervals. If the conditional distribution of the interspikes intervals is unknown, it becomes difficult to evaluate the conditional firing rate. We propose a non parametric estimator for the conditional instantaneous firing rate for Markov, stationary and ergodic ISIs. An algorithm to check the reliability of the proposed estimator is introduced and its consistency properties are proved. The method is applied to data obtained from a stochastic two compartment model and to in vitro experimental data.

preprint2014arXiv

State-dependent Fractional Point Processes

The aim of this paper is the analysis of the fractional Poisson process where the state probabilities $p_k^{ν_k}(t)$, $t\ge 0$, are governed by time-fractional equations of order $0<ν_k\leq 1$ depending on the number $k$ of events occurred up to time $t$. We are able to obtain explicitely the Laplace transform of $p_k^{ν_k}(t)$ and various representations of state probabilities. We show that the Poisson process with intermediate waiting times depending on $ν_k$ differs from that constructed from the fractional state equations (in the case $ν_k = ν$, for all $k$, they coincide with the time-fractional Poisson process). We also introduce a different form of fractional state-dependent Poisson process as a weighted sum of homogeneous Poisson processes. Finally we consider the fractional birth process governed by equations with state-dependent fractionality.

preprint2014arXiv

The role of detachment of in-links in scale-free networks

Real-world networks may exhibit detachment phenomenon determined by the cancelling of previously existing connections. We discuss a tractable extension of Yule model to account for this feature. Analytical results are derived and discussed both asymptotically and for a finite number of links. Comparison with the original model is performed in the supercritical case. The first-order asymptotic tail behavior of the two models is similar but differences arise in the second-order term. We explicitly refer to World Wide Web modeling and we show the agreement of the proposed model on very recent data. However, other possible network applications are also mentioned.

preprint2013arXiv

Coupled systems of fractional equations related to sound propagation: analysis and discussion

In this note we analyse the propagation of a small density perturbation in a one-dimensional compressible fluid by means of fractional calculus modelling, replacing thus the ordinary time derivative with the Caputo fractional derivative in the constitutive equations. By doing so, we embrace a vast phenomenology, including subdiffusive, superdiffusive and also memoryless processes like classical diffusions. From a mathematical point of view, we study systems of coupled fractional equations, leading to fractional diffusion equations or to equations with sequential fractional derivatives. In this framework we also propose a method to solve partial differential equations with sequential fractional derivatives by analysing the corresponding coupled system of equations.

preprint2013arXiv

Fractional Non-Linear, Linear and Sublinear Death Processes

This paper is devoted to the study of a fractional version of non-linear $\mathpzc{M}^ν(t)$, $t>0$, linear $M^ν(t)$, $t>0$ and sublinear $\mathfrak{M}^ν(t)$, $t>0$ death processes. Fractionality is introduced by replacing the usual integer-order derivative in the difference-differential equations governing the state probabilities, with the fractional derivative understood in the sense of Dzhrbashyan--Caputo. We derive explicitly the state probabilities of the three death processes and examine the related probability generating functions and mean values. A useful subordination relation is also proved, allowing us to express the death processes as compositions of their classical counterparts with the random time process $T_{2 ν} (t)$, $t>0$. This random time has one-dimensional distribution which is the folded solution to a Cauchy problem of the fractional diffusion equation.

preprint2013arXiv

On Some Operators Involving Hadamard Derivatives

In this paper we introduce a novel Mittag--Leffler-type function and study its properties in relation to some integro-differential operators involving Hadamard fractional derivatives or Hyper-Bessel-type operators. We discuss then the utility of these results to solve some integro-differential equations involving these operators by means of operational methods. We show the advantage of our approach through some examples. Among these, an application to a modified Lamb--Bateman integral equation is presented.

preprint2013arXiv

On the Integral of Fractional Poisson Processes

In this paper we consider the Riemann--Liouville fractional integral $\mathcal{N}^{α,ν}(t)= \frac{1}{Γ(α)} \int_0^t (t-s)^{α-1}N^ν(s) \, \mathrm ds $, where $N^ν(t)$, $t \ge 0$, is a fractional Poisson process of order $ν\in (0,1]$, and $α> 0$. We give the explicit bivariate distribution $\Pr \{N^ν(s)=k, N^ν(t)=r \}$, for $t \ge s$, $r \ge k$, the mean $\mathbb{E}\, \mathcal{N}^{α,ν}(t)$ and the variance $\mathbb{V}\text{ar}\, \mathcal{N}^{α,ν}(t)$. We study the process $\mathcal{N}^{α,1}(t)$ for which we are able to produce explicit results for the conditional and absolute variances and means. Much more involved results on $\mathcal{N}^{1,1}(t)$ are presented in the last section where also distributional properties of the integrated Poisson process (including the representation as random sums) is derived. The integral of powers of the Poisson process is examined and its connections with generalised harmonic numbers is discussed.

preprint2013arXiv

Parameter estimation for fractional birth and fractional death processes

The fractional birth and the fractional death processes are more desirable in practice than their classical counterparts as they naturally provide greater flexibility in modeling growing and decreasing systems. In this paper, we propose formal parameter estimation procedures for the fractional Yule, the fractional linear death, and the fractional sublinear death processes. The methods use all available data possible, are computationally simple and asymptotically unbiased. The procedures exploited the natural structure of the random inter-birth and inter-death times that are known to be independent but are not identically distributed. We also showed how these methods can be applied to certain models with more general birth and death rates. The computational tests showed favorable results for our proposed methods even with relatively small sample sizes. The proposed methods are also illustrated using the branching times of the plethodontid salamanders data of \cite{hal79}.

preprint2013arXiv

Renewal processes based on generalized Mittag--Leffler waiting times

The fractional Poisson process has recently attracted experts from several fields of study. Its natural generalization of the ordinary Poisson process made the model more appealing for real-world applications. In this paper, we generalized the standard and fractional Poisson processes through the waiting time distribution, and showed their relations to an integral operator with a generalized Mittag--Leffler function in the kernel. The waiting times of the proposed renewal processes have the generalized Mittag--Leffler and stretched-squashed Mittag--Leffler distributions. Note that the generalizations naturally provide greater flexibility in modeling real-life renewal processes. Algorithms to simulate sample paths and to estimate the model parameters are derived. Note also that these procedures are necessary to make these models more usable in practice. State probabilities and other qualitative or quantitative features of the models are also discussed.

preprint2013arXiv

Simulation and estimation for the fractional Yule process

In this paper, we propose some representations of a generalized linear birth process called fractional Yule process (fYp). We also derive the probability distributions of the random birth and sojourn times. The inter-birth time distribution and the representations then yield algorithms on how to simulate sample paths of the fYp. We also attempt to estimate the model parameters in order for the fYp to be usable in practice. The estimation procedure is then tested using simulated data as well. We also illustrate some major characteristics of fYp which will be helpful for real applications.

preprint2013arXiv

Superprocesses as models for information dissemination in the Future Internet

Future Internet will be composed by a tremendous number of potentially interconnected people and devices, offering a variety of services, applications and communication opportunities. In particular, short-range wireless communications, which are available on almost all portable devices, will enable the formation of the largest cloud of interconnected, smart computing devices mankind has ever dreamed about: the Proximate Internet. In this paper, we consider superprocesses, more specifically super Brownian motion, as a suitable mathematical model to analyse a basic problem of information dissemination arising in the context of Proximate Internet. The proposed model provides a promising analytical framework to both study theoretical properties related to the information dissemination process and to devise efficient and reliable simulation schemes for very large systems.

preprint2013arXiv

Transient behavior of fractional queues and related processes

We propose a generalization of the classical M/M/1 queue process. The resulting model is derived by applying fractional derivative operators to a system of difference-differential equations. This generalization includes both non-Markovian and Markovian properties, which naturally provide greater flexibility in modeling real queue systems than its classical counterpart. Algorithms to simulate M/M/1 queue process and the related linear birth-death process are provided. Closed-form expressions of the point and interval estimators of the parameters of these fractional stochastic models are also presented. These methods are necessary to make these models usable in practice. The proposed fractional M/M/1 queue model and the statistical methods are illustrated using S&P data.

preprint2011arXiv

A note on fractional linear pure birth and pure death processes in epidemic models

In this note we highlight the role of fractional linear birth and linear death processes recently studied in \citet{sakhno} and \citet{pol}, in relation to epidemic models with empirical power law distribution of the events. Taking inspiration from a formal analogy between the equation of self consistency of the epidemic type aftershock sequences (ETAS) model, and the fractional differential equation describing the mean value of fractional linear growth processes, we show some interesting applications of fractional modelling to study \textit{ab initio} epidemic processes without the assumption of any empirical distribution. We also show that, in the frame of fractional modelling, subcritical regimes can be linked to linear fractional death processes and supercritical regimes to linear fractional birth processes. Moreover we discuss a simple toy model to underline the possible application of these stochastic growth models to more general epidemic phenomena such as tumoral growth.

preprint2011arXiv

Compositions, Random Sums and Continued Random Fractions of Poisson and Fractional Poisson Processes

In this paper we consider the relation between random sums and compositions of different processes. In particular, for independent Poisson processes $N_α(t)$, $N_β(t)$, $t>0$, we show that $N_α(N_β(t)) \overset{\text{d}}{=} \sum_{j=1}^{N_β(t)} X_j$, where the $X_j$s are Poisson random variables. We present a series of similar cases, the most general of which is the one in which the outer process is Poisson and the inner one is a nonlinear fractional birth process. We highlight generalisations of these results where the external process is infinitely divisible. A section of the paper concerns compositions of the form $N_α(τ_k^ν)$, $ν\in (0,1]$, where $τ_k^ν$ is the inverse of the fractional Poisson process, and we show how these compositions can be represented as random sums. Furthermore we study compositions of the form $Θ(N(t))$, $t>0$, which can be represented as random products. The last section is devoted to studying continued fractions of Cauchy random variables with a Poisson number of levels. We evaluate the exact distribution and derive the scale parameter in terms of ratios of Fibonacci numbers.

preprint2011arXiv

Fractional pure birth processes

We consider a fractional version of the classical nonlinear birth process of which the Yule--Furry model is a particular case. Fractionality is obtained by replacing the first order time derivative in the difference-differential equations which govern the probability law of the process with the Dzherbashyan--Caputo fractional derivative. We derive the probability distribution of the number $\mathcal{N}_ν(t)$ of individuals at an arbitrary time $t$. We also present an interesting representation for the number of individuals at time $t$, in the form of the subordination relation $\mathcal{N}_ν(t)=\mathcal{N}(T_{2ν}(t))$, where $\mathcal{N}(t)$ is the classical generalized birth process and $T_{2ν}(t)$ is a random time whose distribution is related to the fractional diffusion equation. The fractional linear birth process is examined in detail in Section 3 and various forms of its distribution are given and discussed.

preprint2011arXiv

On a fractional linear birth--death process

In this paper, we introduce and examine a fractional linear birth--death process $N_ν(t)$, $t>0$, whose fractionality is obtained by replacing the time derivative with a fractional derivative in the system of difference-differential equations governing the state probabilities $p_k^ν(t)$, $t>0$, $k\geq0$. We present a subordination relationship connecting $N_ν(t)$, $t>0$, with the classical birth--death process $N(t)$, $t>0$, by means of the time process $T_{2ν}(t)$, $t>0$, whose distribution is related to a time-fractional diffusion equation. We obtain explicit formulas for the extinction probability $p_0^ν(t)$ and the state probabilities $p_k^ν(t)$, $t>0$, $k\geq1$, in the three relevant cases $λ>μ$, $λ<μ$, $λ=μ$ (where $λ$ and $μ$ are, respectively, the birth and death rates) and discuss their behaviour in specific situations. We highlight the connection of the fractional linear birth--death process with the fractional pure birth process. Finally, the mean values $\mathbb{E}N_ν(t)$ and $\operatorname {\mathbb{V}ar}N_ν(t)$ are derived and analyzed.

preprint2011arXiv

Randomly Stopped Nonlinear Fractional Birth Processes

We present and analyse the nonlinear classical pure birth process $\mathpzc{N} (t)$, $t>0$, and the fractional pure birth process $\mathpzc{N}^ν(t)$, $t>0$, subordinated to various random times, namely the first-passage time $T_t$ of the standard Brownian motion $B(t)$, $t>0$, the $α$-stable subordinator $\mathpzc{S}^α(t)$, $α\in (0,1)$, and others. For all of them we derive the state probability distribution $\hat{p}_k (t)$, $k \geq 1$ and, in some cases, we also present the corresponding governing differential equation. We also highlight interesting interpretations for both the subordinated classical birth process $\hat{\mathpzc{N}} (t)$, $t>0$, and its fractional counterpart $\hat{\mathpzc{N}}^ν(t)$, $t>0$ in terms of classical birth processes with random rates evaluated on a stretched or squashed time scale. Various types of compositions of the fractional pure birth process $\mathpzc{N}^ν(t)$ have been examined in the last part of the paper. In particular, the processes $\mathpzc{N}^ν(T_t)$, $\mathpzc{N}^ν(\mathpzc{S}^α(t))$, $\mathpzc{N}^ν(T_{2ν}(t))$, have been analysed, where $T_{2ν}(t)$, $t>0$, is a process related to fractional diffusion equations. Also the related process $\mathpzc{N}(\mathpzc{S}^α({T_{2ν}(t)}))$ is investigated and compared with $\mathpzc{N}(T_{2ν}(\mathpzc{S}^α(t))) = \mathpzc{N}^ν(\mathpzc{S}^α(t))$. As a byproduct of our analysis, some formulae relating Mittag--Leffler functions are obtained.

preprint2011arXiv

The Space-Fractional Poisson Process

In this paper we introduce the space-fractional Poisson process whose state probabilities $p_k^α(t)$, $t>0$, $α\in (0,1]$, are governed by the equations $(\mathrm d/\mathrm dt)p_k(t) = -λ^α(1-B)p_k^α(t)$, where $(1-B)^α$ is the fractional difference operator found in the study of time series analysis. We explicitly obtain the distributions $p_k^α(t)$, the probability generating functions $G_α(u,t)$, which are also expressed as distributions of the minimum of i.i.d.\ uniform random variables. The comparison with the time-fractional Poisson process is investigated and finally, we arrive at the more general space-time fractional Poisson process of which we give the explicit distribution.