Fractional Poisson process

In probability theory, a fractional Poisson process is a stochastic process to model the long-memory dynamics of a stream of counts. The time interval between each pair of consecutive counts follows the non-exponential power-law distribution with parameter \nu , which has physical dimension [\nu ]=\sec ^{-\mu }, where  0<\mu \leq 1. In other words, fractional Poisson process is non-Markov counting stochastic process which exhibits non-exponential distribution of interarrival times. The fractional Poisson process is a continuous-time process which can be thought of as natural generalization of the well-known Poisson process. Fractional Poisson probability distribution is a new member of discrete probability distributions.

The fractional Poisson process, Fractional compound Poisson process and fractional Poisson probability distribution function have been invented, developed and encouraged for applications by Nick Laskin (2003) who coined the terms fractional Poisson process, Fractional compound Poisson process and fractional Poisson probability distribution function.[1]

Fundamentals

The fractional Poisson probability distribution captures the long-memory effect which results in the non-exponential waiting time probability distribution function empirically observed in complex classical and quantum systems. Thus, fractional Poisson process and fractional Poisson probability distribution function can be considered as natural generalization of the famous Poisson process and the Poisson probability distribution.

The idea behind the fractional Poisson process was to design counting process with non-exponential waiting time probability distribution. Mathematically the idea was realized by substitution the first-order time derivative in the Kolmogorov–Feller equation for the Poisson probability distribution function with the time derivative of fractional order.[2][3]

The main outcomes are new stochastic non-Markov process – fractional Poisson process and new probability distribution function – fractional Poisson probability distribution function.

Fractional Poisson probability distribution function

The probability distribution function of fractional Poisson process has been found for the first time by Nick Laskin (see, Ref.[1])


P_\mu (n,t)=\frac{(\nu t^\mu )^n}{n!}\sum\limits_{k=0}^\infty \frac{(k+n)!}{
k!}\frac{(-\nu t^\mu )^k}{\Gamma (\mu (k+n)+1)},\qquad 0<\mu \leq 1,

where parameter \nu has physical dimension [\nu ]=\sec ^{-\mu } and {\Gamma (\mu (k+n)+1)} is the Gamma function.

The P_\mu (n,t) gives us the probability that in the time interval [0,t] we observe n events governed by fractional Poisson stream.

The probability distribution of the fractional Poisson process P_{\mu }(n,t) can be represented in terms of the Mittag-Leffler function E_{\mu }(z) in the following compact way (see, Ref.[1]),


P_{\mu }(n,t)=(\frac{(-z)^{n}}{n!}\frac{d^{n}}{dz^{n}}E_{\mu }(z))|_{z=-\nu
t^{\mu }},


P_{\mu }(n=0,t)=E_{\mu }(-\nu t^{\mu }).

It follows from the above equations that when \mu =1 the P_\mu (n,t) is transformed into the well known probability distribution function of the Poisson process, P(n,t)=P_1(n,t), 
P(n,t)=\frac{(\overline{\nu}t )^n}{n!}\exp(-\overline{\nu} t),


P(n=0,t)=\exp(-\overline{\nu} t),

where \overline{\nu} is the rate of arrivals with physical dimension [\overline{\nu}]=\sec ^{-1 }.

Thus, P_\mu (n,t) can be considered as fractional generalization of the standard Poisson probability distribution. The presence of additional parameter \mu brings new features in comparison with the standard Poisson distribution.

Mean

The mean \overline{n}_\mu of the fractional Poisson process has been found in Ref.[1].


\overline{n}_\mu =\sum\limits_{n=0}^\infty nP_\mu (n,t)=\frac{\nu t^\mu }{
\Gamma (\mu +1)}.

The second order moment

The second order moment of the fractional Poisson process \overline{n^2}_\mu has been found for the first time by Nick Laskin (see, Ref.[1])


\overline{n_\mu ^2}=\sum\limits_{n=0}^\infty n^2P_\mu (n,t)=\overline{n}_\mu
+\overline{n}_\mu ^2\frac{\sqrt{\pi }\Gamma (\mu +1)}{2^{2\mu -1}\Gamma (\mu
+\frac 12)}.

Variance

The variance of the fractional Poisson process is (see, Ref.[1])


\sigma _\mu =\overline{n_\mu ^2}-\overline{n}_\mu ^2=\overline{n}_\mu + \overline{n}_\mu ^2\left\{ \frac{\mu B(\mu ,\frac 12)}{2^{2\mu -1}} - 1\right\},

where B(\mu ,\frac 12) is the Beta-function.

Characteristic function

The characteristic function of the fractional Poisson process has been found for the first time in Ref.[1],


C_\mu (s,t)=\sum\limits_{n=0}^\infty e^{isn}P_\mu (n,t)=E_\mu (\nu t^\mu (e^{is}-1)).

or in a series form


C_\mu (s,t)=\sum\limits_{m=0}^\infty \frac 1{\Gamma (m\mu +1)}\left( \nu
t^\mu (e^{is}-1)\right) ^m,

with the help of the Mittag-Leffler function series representation.

Then, for the moment of k^{{\rm th}} order we have


\overline{n_\mu ^k}= (1/i^k)\frac{\partial ^kC_\mu (s,t)}{\partial s^k}|_{s=0}.

Generating function

The generating function G_\mu(s,t) of the fractional Poisson probability distribution function is defined as (see, Ref.[1]).


G_\mu (s,t)=\sum\limits_{n=0}^\infty s^nP_\mu (n,t).

The generating function of the fractional Poisson probability distribution was obtained for the first time by Nick Laskin in Ref.[1].


G_\mu (s,t)=E_\mu (\nu t^\mu (s-1)),

where E_\mu (z) is the Mittag-Leffler function given by its series representation


E_\mu (z)=\sum\limits_{m=0}^\infty \frac{z^m}{\Gamma (\mu m+1)}.

Moment generating function

The equation for the moment of any integer order of the fractional Poisson can be easily found by means of the moment generating function H_\mu (s,t) which is defined as


H_\mu (s,t)=\sum\limits_{n=0}^\infty e^{-sn}P_\mu (n,t).

For example, for the moment of k^{{\rm th}} order we have


\overline{n_\mu ^k}=(-1)^k\frac{\partial ^kH_\mu (s,t)}{\partial s^k}|_{s=0}.

The moment generating function H_\mu (s,t) is (see, Ref.[1])


H_\mu (s,t)=E_\mu (\nu t^\mu (e^{-s}-1)),

or in a series form


H_\mu (s,t)=\sum\limits_{m=0}^\infty \frac 1{\Gamma (m\mu +1)}\left( \nu
t^\mu (e^{-s}-1)\right) ^m,

with the help of the Mittag-Leffler function series representation.

Waiting time distribution function

A time between two successive arrivals is called as waiting time and it is a random variable. The waiting time probability distribution function is an important attribute of any arrival or counting random process.

Waiting time probability distribution function \psi _\mu (\tau) of the fractional Poisson process is defined as (see, Refs.[1,3])


\psi _\mu (\tau )=-\frac d{d\tau }P_\mu (\tau ),

where P_\mu (\tau ) is the probability that a given interarrival time is greater or equal to \tau


P_\mu (\tau )=1-\sum\limits_{n=1}^\infty P_\mu (n,\tau )=E_\mu (-\nu \tau^\mu ),

and P_\mu (n,\tau ) is the fractional Poisson probability distribution function.

The waiting time probability distribution function \psi _\mu (\tau) of the fractional Poisson process was found for first time by Nick Laskin in Ref.[1],


\psi _\mu (\tau )=\nu \tau ^{\mu -1}E_{\mu ,\mu }(-\nu \tau ^\mu ),\qquad
t\geq 0,\qquad 0<\mu \leq 1,

here E_{\alpha ,\beta }(z) is the generalized two-parameter Mittag-Leffler function


E_{\alpha ,\beta }(z)=\sum\limits_{m=0}^\infty \frac{z^m}{\Gamma (\alpha
m+\beta )},\qquad E_{\alpha ,1}(z)=E_\alpha (z).

Waiting time probability distribution function \psi _\mu (\tau ) has the following asymptotic behavior (see, Ref.[1])


\psi _{\mu }(\tau )\simeq 1/\nu \tau ^{\mu +1},\qquad \tau \rightarrow
\infty ,

and


\psi _{\mu }(\tau )\simeq \nu \tau ^{\mu -1},\qquad \tau \rightarrow 0.

Fractional compound Poisson process

Fractional compound Poisson process has been introduced and developed for the first time by Nick Laskin (see, Ref.[1]). The fractional compound Poisson process \{X(t), t\geq 0\} is represented by


X(t)=\sum\limits_{i=1}^{N(t)}Y_i,

where \{N(t), t\geq 0\} is a fractional Poisson process, and \{Y_i, i=1,2,\ldots\} is a family of independent and identically distributed random variables with probability distribution function p(Y) for each Y_i. The process \{N(t), t\geq 0\} and the sequence \{Y_i, i=1,2,\ldots\} are assumed to be independent.

The fractional compound Poisson process is natural generalization of the compound Poisson process.

Applications of fractional Poisson probability distribution

The fractional Poisson probability distribution has physical and mathematical applications. Physical application is in the field of quantum optics. Mathematical applications are in the field of combinatorial numbers (see, Ref.[4]).

Physical application: New coherent states

A new family of quantum coherent states |\varsigma > has been introduced as[4]


|\varsigma \rangle=\sum\limits_{n=0}^{\infty }\frac{(\sqrt{\mu }\varsigma ^{\mu
})^{n}}{\sqrt{n!}}(E_{\mu }^{(n)}(-\mu |\varsigma |^{2\mu }))^{1/2}|n\rangle,

where |n\rangle is an eigenvector of the photon number operator, complex number \varsigma stands for labeling the new coherent states,


E_\mu ^{(n)}(-\mu |\varsigma |^{2\mu })=\frac{d^n}{dz^n}E_\mu (z)|_{z=-\mu
|\varsigma |^{2\mu }}

and E_\mu(x) is the Mittag-Leffler function.

Then the probability P_{\mu }(n) of detecting n photons is:


P_{\mu }(n)=|\langle n|\varsigma \rangle |^{2}=\frac{(\mu |\varsigma |^{2\mu })^{n}}{n!}
\left( E_{\mu }^{(n)}(-\mu |\varsigma |^{2\mu })\right) ,

which is recognized as fractional Poisson probability distribution.

In terms of photon field creation and annihilation operators a^{+} and a that satisfy the canonical commutation relation [a,a^{+}]=aa^{+}-a^{+}a=1, the average number of photons \bar n in a coherent state |\varsigma \rangle can be presented as (see, Ref.[4])


\bar n=\langle \varsigma |a^{+}a|\varsigma \rangle =\sum\limits_{n=0}^\infty nP_\mu
(n)=(\mu |\varsigma |^{2\mu })/\Gamma (\mu +1)
.

Mathematical applications: New polynomials and numbers

The fractional generalization of Bell polynomials, Bell numbers, Dobinski's formula and Stirling numbers of the second kind have been introduced and developed by Nick Laskin (see, Ref.[4]). The appearance of fractional Bell polynomials is natural if one evaluates the diagonal matrix element of the evolution operator in the basis of newly introduced quantum coherent states. Fractional Stirling numbers of the second kind have been applied to evaluate the skewness and kurtosis of the fractional Poisson probability distribution function. A new representation of the Bernoulli numbers in terms of fractional Stirling numbers of the second kind has been discovered (see, Ref.[4]).

In the limit case μ =1 when the fractional Poisson probability distribution becomes the Poisson probability distribution, all of the above listed applications turn into the well-known results of the quantum optics and the enumerative combinatorics.

Statistical application and inference

The point and interval estimators for the model parameters are developed by Cahoy et. al, (2010) (see, Ref.[5]). [5]

See also

References

  1. N. Laskin, (2003), http://dx.doi.org/10.1016/S1007-5704(03)00037-6 Fractional Poisson process, Communications in Nonlinear Science and Numerical Simulation, vol. 8 issue 3–4 September–December, 2003. pp. 201–213.
  2. A.I. Saichev and G.M. Zaslavsky, (1997), http://dx.doi.org/10.1063/1.166272 Fractional kinetic equations: solutions and applications, Chaos vol. 7 (1997) pp. 753–764.
  3. O. N. Repin and A. I. Saichev, (2000), http://www.springerlink.com/content/r88713p577701148 Fractional Poisson Law, Radiophysics and Quantum Electronics, vol 43, Number 9 (2000), 738-741.
  4. N. Laskin, (2009), Some applications of the fractional Poisson probability distribution, J. Math. Phys. 50, 113513 (2009) (12 pages), http://jmp.aip.org/resource/1/jmapaq/v50/i11/p113513_s1?bypassSSO=1. (also available online: http://arxiv.org/abs/0812.1193)
  5. D.O. Cahoy V.V. Uchaikin W.A. Woyczyński (2010). Parameter estimation for fractional Poisson processes. Journal of Statistical Planning and Inference. pp. 3106–3120.

Further reading

This article is issued from Wikipedia - version of the Wednesday, April 20, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.