---
abstract: 'The integro-differential wave equation for the probability density function for a classical one-dimensional Lévy walk with continuous sample paths has been derived. This equation involves a classical wave operator together with memory integrals describing the spatio-temporal coupling of the Lévy walk. It is valid for any running time PDF of a walker and it does not involve any long-time large-scale approximations. It generalizes the well-known telegraph equation obtained from the persistent random walk. Several non-Markovian cases are considered when the particle’s velocity alternates at the gamma and power-law distributed random times.'
author:
- Sergei Fedotov
title: 'Single integro-differential wave equation for Lévy walk '
---

Lévy walk is a fundamental notion in physics and biology [@KlSo] with numerous applications including the transport of light governed by Lévy statistics [@La], anomalous superdiffusion of cold atoms in optical lattices [@Ba], T-cell motility in the brain [@Le2], endosomal active transport within living cells [@Gr]. In the last few years, the interest in the Lévy walk models increases rapidly in the context of Lévy flight foraging hypothesis according to which a Lévy transport is optimal to search for randomly located objects [@Le4]. Lévy movement pattern has been observed in microorganisms, insects, molluscs, birds, etc. [@Men]. Recent publication *Lévy walk* in* Review of Modern Physics* provides a detailed discussion of existing applications and the most current status of Lévy walk theory [@Za].

The derivation of the governing equations for a Lévy walk from an underlying stochastic movement is a long-standing problem. Many contributions on this subject have been given since the pioneering works three decades ago (see, for example, [@Kl; @Blu]). The standard model for a Lévy walk is based on the continuos time random walk (CTRW) with coupled probability density function for the jump length and the waiting time between two successive jumps. Two integral equations for the probability density function (PDF) $p(x,t)$ for a walker position $x$ at time $t$ and arrival rate $j\left( x,t\right) $ can be formulated and solved by the Fourier-Laplace transform technique [@KlSo; @Za]. An equivalent single integral equation for $p(x,t)$ has been also formulated [@Zabur]. To describe a Lévy walk spatiotemporal coupling, two dynamical equations involving fractional material derivatives have been suggested in [@SokR]. Another approach to describe the superdiffusive behavior is based on the analysis of joint PDF $p(x,v,t)$ of the particle’s position $x$ and velocity $v$. Various fractional generalizations of Kramers-Fokker-Planck equation for $p(x,v,t)$  have been derived [@Kra1; @Kra2; @Kra3]. However, to the author’s knowledge, the single integro-differential equation for the Lévy walk PDF $p(x,t)$ is still not available in the literature. In this paper we obtain the following integro-differential wave equation $$\frac{\partial ^{2}p}{\partial t^{2}}-v^{2}\frac{\partial ^{2}p}{\partial
x^{2}}+\int_{0}^{t}\int_{V}K\left( \tau \right) \varphi (u)\left( \frac{%
\partial }{\partial t}-u\frac{\partial }{\partial x}\right) \times$$$$p\left( x-u\tau ,t-\tau \right) dud\tau =0,  \label{main}$$where $v$ is a constant speed of walker, $\varphi (u)$ is the velocity jump density:$$\varphi (u)=\frac{1}{2}\delta \left( u-v\right) +\frac{1}{2}\delta \left(
u+v\right)  \label{jump}$$in the velocity space $V$ and $K(\tau )$ is the standard memory kernel from the CTRW theory [@KlSo]. It is determined by its Laplace transform $\hat{%
K}(s)=\hat{\psi}(s)/\hat{\Psi}(s),$ where $\hat{\psi}(s)$ and $\hat{\Psi}(s)$ are the Laplace transforms of the running time density $\psi (\tau )$ andthe survival function $\Psi (\tau )$. From Eq. (\[main\]) with $%
p(x,0)=p_{0}\left( x\right) $ and $p_{t}(x,0)=0,$ one can obtain the well-known expression for the Fourier-Laplace transform of the PDF $p(x,t)$ $$\hat{p}(k,s)=\frac{\left[ \hat{\Psi}(s+ikv)+\hat{\Psi}(s-ikv)\right] \hat{p}%
_{0}(k)}{2-\hat{\psi}(s+ikv)-\hat{\psi}(s-ikv)},  \label{LF}$$where $\hat{p}_{0}(k)=\int_{\mathbb{R}}p_{0}\left( x\right) e^{ikx}dx$ [Za]{}.

*Derivation.* We consider the Lévy walk as the random particle’s motion with continuos sample paths (no jumps) along one-dimensional space. The particle starts to move with constant speed $v$ at time $t=0$ and after a random time (running time) it either continues the movement in the same direction with curtain probability or changes the direction and moves with the same constant speed. The random running time is defined by the switching rate $\gamma (\tau )$ or the running time PDF $\psi (\tau )=\gamma (\tau
)\exp \left[ -\int_{0}^{\tau }\gamma (s)ds\right] .$ To derive the governing equation for the Lévy walk, we start with Markovian model involving structural densities with the extra running time variable $\tau $ [MFH,Cox,Vlad,Fer,FedotAbby]{}. We define the structural PDF’s of walker, $%
n_{+}(x,t,\tau ),$ at point $x$ and time $t$ that moves in the right direction, $\left( +\right) ,$ with constant speed $v$ during time $\tau $ since the last switching. The probability density function $n_{-}(x,t,\tau )$ corresponds to the walker that moves in the negative direction, $\left(
-\right) .$ The balance equations for both structural PDF’s $n_{+}(x,t,\tau
) $ and $n_{-}(x,t,\tau )$ can be written as $$\frac{\partial n_{\pm }}{\partial t}\pm v\frac{\partial n_{\pm }}{\partial x}%
+\frac{\partial n_{\pm }}{\partial \tau }=-\gamma (\tau )n_{\pm },
\label{str}$$where the switching rate $\gamma (\tau )$ depends on the running time $\tau $. If the walker moves in the positive direction it can switch with rate $%
\gamma (\tau )$ either to the opposite direction with the probability $%
\alpha _{-}$ or keep the same direction with the probability $\alpha _{+}$ such that $\alpha _{+}+\alpha _{-}=1.$ For the walker moving in the negative direction corresponding characteristics are $\beta _{+}$ and $\beta
_{-}.$ The well-known *velocity model* and *two-state model* are just particular cases of this general two-state model. For example, the choice $\alpha _{+}=\beta _{-}=0$ and $\alpha _{-}=\beta _{+}=1$ corresponds to the two-state model [@Mas]. The probabilities $\alpha _{\pm }=\beta
_{\pm }=1/2$ correspond to the velocity model [@Zu]. We assume that at the initial time $t=0$ all walkers start to move with zero running time $%
\tau $ $$n_{\pm }(x,0,\tau )=\frac{1}{2}p_{0}(x)\delta (\tau ).  \label{initial}$$Here we consider the symmetrical initial conditions when walker start to move on the right with probability $1/2$ and on the left with the same probability. The boundary conditions at zero running time can be formulated as $$\begin{aligned}
n_{\pm }(x,t,0) &=&\alpha _{\pm }\int_{0}^{t}\gamma (\tau )n_{+}(x,t,\tau
)d\tau +  \notag \\
&&\beta _{\pm }\int_{0}^{t}\gamma (\tau )n_{-}(x,t,\tau )d\tau .
\label{zero2}\end{aligned}$$Our aim is to derive the master equations for the probability density functions $p_{+}(x,t)$ and $p_{-}(x,t)$ defined as$$p_{\pm }(x,t)=\int_{0}^{t}n_{\pm }(x,t,\tau )d\tau .  \label{den}$$By differentiating (\[den\]) with respect to time $t$ and using the balance equations (\[str\]) we obtain $$\begin{aligned}
\frac{\partial p_{\pm }}{\partial t} &=&n_{\pm }(x,t,t)\mp v\int_{0}^{t}%
\frac{\partial n_{\pm }}{\partial x}d\tau \\
&&-\int_{0}^{t}\frac{\partial n_{\pm }}{\partial \tau }d\tau
-\int_{0}^{t}\gamma (\tau )n_{\pm }(x,t,\tau )d\tau .\end{aligned}$$This equation can be rewritten as the master equation$$\frac{\partial p_{\pm }}{\partial t}\pm v\frac{\partial p_{\pm }}{\partial x}%
=j_{\pm }(x,t)-i_{\pm }(x,t),  \label{mas}$$where the rates of switching $i_{\pm }(x,t)$ are defined as $$i_{\pm }(x,t)=\int_{0}^{t}\gamma (\tau )n_{\pm }(x,t,\tau )d\tau  \label{i0}$$and the arrival rates $j_{\pm }(x,t)$ are $$j_{\pm }(x,t)=n_{\pm }(x,t,0).  \label{ar}$$By using the definitions of switching and arrival rates (\[i0\]) and ([ar]{}), Eq. (\[zero2\]) can be written in the compact form$$j_{\pm }(x,t)=\alpha _{\pm }i_{+}(x,t)+\beta _{\pm }i_{-}(x,t).  \label{ra}$$In what follows we consider only the simple case of a symmetric Lévy walk for which $\alpha _{\pm }=\beta _{\pm }=1/2.$ In general, the probabilities $\alpha _{\pm }$ and $\beta _{\pm }$ can be useful to formulate the impact of the external force or chemotactic substance. Substitution of (\[ra\]) with $\alpha _{\pm }=\beta _{\pm }=1/2$ into the master equation (\[mas\]) gives $$\frac{\partial p_{+}}{\partial t}+v\frac{\partial p_{+}}{\partial x}=-\frac{1%
}{2}i_{+}(x,t)+\frac{1}{2}i_{-}(x,t),  \label{mean3}$$$$\frac{\partial p_{-}}{\partial t}-v\frac{\partial p_{-}}{\partial x}=\frac{1%
}{2}i_{+}(x,t)-\frac{1}{2}i_{-}(x,t),  \label{mean4}$$where the switching rates $i_{\pm }(x,t)$ can be found as follows. By the method of characteristics, we find from (\[str\]) $$n_{\pm }(x,t,\tau )=n_{\pm }(x\mp v\tau ,t-\tau ,0)\Psi (\tau ),\quad \tau
<t,  \label{cha}$$where $\Psi (\tau )$ is the survival function $$\Psi (\tau )=e^{-\int_{0}^{\tau }\gamma (s)ds}.  \label{surva}$$Note that at $\tau =t$ we have a singularity due to the initial condition (\[initial\]). Substitution of (\[cha\]) into (\[i0\]) and (\[den\]) together with the initial condition (\[initial\]) gives $$i_{\pm }(x,t)=\int_{0}^{t}j_{\pm }(x\mp v\tau ,t-\tau )\psi (\tau )d\tau +%
\frac{1}{2}p_{0}(x\mp vt)\psi (t),$$$$p_{\pm }(x,t)=\int_{0}^{t}j_{\pm }(x\mp v\tau ,t-\tau )\Psi \left( \tau
\right) d\tau +\frac{1}{2}p_{0}(x\mp vt)\Psi \left( t\right) .$$Applying the Fourier-Laplace transform to these equations, we can find expressions for $i_{+}(x,t)$ and $i_{-}(x,t)$ in terms of $\rho _{+}(x,t)$ and $\rho _{-}(x,t)$ as [@FedotAbby] $$i_{\pm }(x,t)=\int_{0}^{t}K(\tau )p_{\pm }(x\mp v\tau ,t-\tau )d\tau .
\label{rate_ii}$$Note that four integral equations $i_{\pm }(x,t)$ and $p_{\pm }(x,t)$ above can be written in the standard form of two equations for the PDF:$$p(x,t)=p_{+}(x,t)+p_{-}(x,t).  \label{total}$$and the arrival rate $j=j_{+}+j_{-}$ with the jump density $w(z|\tau )=\frac{%
1}{2}\delta \left( z-v\tau \right) +\frac{1}{2}\delta \left( z+v\tau \right)
.$ For the symmetrical case $j=i_{+}+i_{-}$ we obtain the well-known equations$$\begin{aligned}
p(x,t) &=&\int_{0}^{t}\int_{\mathbb{R}}j(x-z,t-\tau )w(z|\tau )\Psi \left(
\tau \right) dzd\tau + \\
&&\int_{\mathbb{R}}p_{0}(x-z)w(z|t)dz\Psi \left( t\right) ,\end{aligned}$$$$\begin{aligned}
j(x,t) &=&\int_{0}^{t}\int_{\mathbb{R}}j(x-z,t-\tau )w(z|\tau )\psi (\tau
)dzd\tau + \\
&&\int_{\mathbb{R}}p_{0}(x-z)w(z|t)dz\psi (t).\end{aligned}$$Note that these two equations are the starting point for the various studies of the Lévy walk [@Za].

Our main purpose is to reduce the system (\[mean3\]), (\[mean4\]) with (\[rate\_ii\]) to a single governing equation for the PDF $p(x,t)$ defined by (\[total\]). First we introduce the flux [@Kac] $$J(x,t)=vp_{+}(x,t)-vp_{-}(x,t).  \label{flux}$$Then by adding (\[mean3\]) and (\[mean4\]) we obtain the standard conservation equation $$\frac{\partial p}{\partial t}+\frac{\partial J}{\partial x}=0.  \label{ro}$$Equation for the flux $J$ can be obtain by multiplication of (\[mean3\]) and (\[mean4\]) by $v$ and subtraction $$\frac{\partial J}{\partial t}+v^{2}\frac{\partial p}{\partial x}=-v\left[
i_{+}(x,t)-i_{-}(x,t)\right] .  \label{J}$$By differentiating (\[ro\]) with respect to $t$ and (\[J\]) with respect to $x$ and eliminating $\partial ^{2}J/\partial t\partial x,$ we obtain$$\frac{\partial ^{2}p}{\partial t^{2}}=v^{2}\frac{\partial ^{2}p}{\partial
x^{2}}+v\frac{\partial }{\partial x}\left[ i_{+}(x,t)-i_{-}(x,t)\right] .
\label{kkk}$$Now we need to express the last term in (\[kkk\]) in terms of $p(x,t)$ alone. From (\[total\]) and (\[flux\]) we find the expressions for $%
p_{+}(x,t)$ and $p_{-}(x,t)$ in terms of the PDF $p(x,t)$ and the flux $%
J(x,t)$ $$p_{\pm }(x,t)=\frac{p(x,t)}{2}\pm \frac{J(x,t)}{2v}.  \label{+}$$Substitution of (\[+\]) into (\[rate\_ii\]) gives $$v\frac{\partial }{\partial x}\left[ i_{+}(x,t)-i_{-}(x,t)\right] =$$$$\begin{aligned}
&&\frac{1}{2}\int_{0}^{t}K\left( \tau \right) [\frac{\partial J}{\partial x}%
(x-v\tau ,t-\tau )+\frac{\partial J}{\partial x}(x+v\tau ,t-\tau )  \notag \\
&&+v\frac{\partial p}{\partial x}(x-v\tau ,t-\tau )-v\frac{\partial p}{%
\partial x}(x+v\tau ,t-\tau )]d\tau .  \label{iiii}\end{aligned}$$By using (\[iiii\]) together with (\[ro\]) and (\[kkk\]), we obtain a single integro-differential equation for the PDF $p(x,t)$ $$\frac{\partial ^{2}p}{\partial t^{2}}=v^{2}\frac{\partial ^{2}p}{\partial
x^{2}}-$$$$\begin{aligned}
&&\frac{1}{2}\int_{0}^{t}K\left( \tau \right) \left[ \left( \frac{\partial }{%
\partial t}-v\frac{\partial }{\partial x}\right) p\left( x-v\tau ,t-\tau
\right) \right] d\tau -  \notag \\
&&\frac{1}{2}\int_{0}^{t}K\left( \tau \right) \left[ \left( \frac{\partial }{%
\partial t}+v\frac{\partial }{\partial x}\right) p\left( x+v\tau ,t-\tau
\right) \right] d\tau .  \label{cla}\end{aligned}$$This integro-differential wave equation for a *Lévy walk* is valid for any running time PDF $\psi \left( \tau \right) $. It generalizes the well-known telegraph equation obtained from the persistent random walk with the constant rate of switching $\gamma .$ By using the velocity jump density (\[jump\]), we rewrite Eq. (\[cla\]) in the compact form ([main]{}). It is instructive to compare the uncoupled CTRW with Lévy walk. When jumps and waiting times of the CTRW are independent random variables, one can convert the single integral equation for $p(x,t)$ into a master equation [@KlSo; @MFH]. For the Lévy walk involving spatio-temporal coupling the integro-differential wave equation (\[cla\]) plays the same role as the master equation for the uncoupled CTRW.

Let us now consider several examples of the running time PDF $\psi \left(
\tau \right) .$

*Exponential running time density.* In the Markovian case with the exponential running time PDF $$\psi \left( \tau \right) =\ \frac{1}{T}\exp \left( -\frac{\tau }{T}\right)
\label{expo}$$for which $\hat{\psi}(s)=\left( 1+Ts\right) ^{-1}$ and $K\left( \tau \right)
=T^{-1}\delta \left( \tau \right) $, we obtain from (\[cla\]) the classical Cattaneo or telegraph equation [@MFH] $$\frac{\partial ^{2}p}{\partial t^{2}}+\frac{1}{T}\frac{\partial p}{\partial t%
}-v^{2}\frac{\partial ^{2}p}{\partial x^{2}}=0.  \label{Cat}$$This hyperbolic equation ensures that the density profile propagates with finite speed $v.$

*Gamma PDF* $g(\tau ,2,\lambda ).$* *For the biological applications it is important to consider a running time PDF that takes the maximum value not at zero time as the exponential density (\[expo\]) [Men]{}. For example, it was found that the running time density for a single bacterium might deviate significantly from exponential approximation [Kor]{}. One example of such PDF is the *gamma density* * *$$\psi \left( \tau \right) =g(\tau ,2,\lambda )=\ \lambda ^{2}\tau \exp \left(
-\lambda \tau \right)  \label{Ga}$$with  $\hat{\psi}(s)=\lambda ^{2}$ $\left( s+\lambda \right) ^{-2}$ and $%
\hat{K}(s)=\lambda ^{2}$ $\left( 2\lambda +s\right) ^{-1}.$ The memory kernel in Eq. (\[cla\]) takes the form $$K\left( \tau \right) =\lambda ^{2}\exp \left( -2\lambda \tau \right) .$$The advantage of this exponential memory kernel is that one can localize integro-differential equation (\[cla\]) by direct differentiation of ([cla]{}) with respect to time twice: $$\frac{\partial ^{4}p}{\partial t^{4}}+4\lambda \frac{\partial ^{3}p}{%
\partial t^{3}}+5\lambda ^{2}\frac{\partial ^{2}p}{\partial t^{2}}+2\lambda
^{3}\frac{\partial p}{\partial t}=$$$$v^{2}\frac{\partial ^{2}}{\partial x^{2}}\left[ 2\frac{\partial ^{2}p}{%
\partial t^{2}}+4\lambda \frac{\partial p}{\partial t}+3\lambda ^{2}p-v^{2}%
\frac{\partial ^{2}p}{\partial x^{2}}\right] .  \label{four}$$Note that non-Markovian particle’s movement with velocities alternating at Erlang-distributed and gamma-distributed random times was considered in [Ital0,Ital]{}. Next we consider the anomalous case involving walker’s velocities alternating at power-law distributed random times [SokR,Fer,Mas]{}.

*Anomalous enhanced transport.* We consider two anomalous cases: (1) strong ballistic case for which the mean squared displacement: $<x^{2}>\sim
t^{2}$ and (2) subballistic superdiffusion with $<x^{2}>\sim t^{3-\mu },$ where $1<\mu <2$ [@Za]. In the ballistic case, we use the survival function $\Psi \left( \tau \right) =E\left[ -\left( \tau /\tau _{0}\right)
^{\mu }\right] $ with $0<\mu <1$ for which the mean running time is divergent [@Sca; @Steve]. In this case $$\hat{\psi}\left( s\right) =\frac{1}{1+\left( \tau _{0}s\right) ^{\mu }}%
,\qquad 0<\mu <1$$and the Laplace transform of the memory kernel $K\left( \tau \right) $ is $$\hat{K}(s)=\frac{s\hat{\psi}(s)}{1-\hat{\psi}(s)}=\frac{s^{1-\mu }}{\tau
_{0}^{\mu }}.$$For this kernel the main equation (\[main\]) can be rewritten in the different forms by using material fractional derivatives [SokR,Mer0,Bae,Mer1,Sib,Mag]{}. We write it in the form $$\frac{\partial ^{2}p}{\partial t^{2}}-v^{2}\frac{\partial ^{2}p}{\partial
x^{2}}+\frac{1}{\tau _{0}^{\mu }}\int_{V}\varphi (u)L_{u}^{1-\mu }pdu=0,
\label{lee}$$where the operator $L_{u}^{1-\mu }$ is defined by its Fourier-Laplace transform $$\mathcal{FL}\left\{ L_{u}^{1-\mu }p\right\} =\left( s-iku\right) ^{1-\mu }%
\left[ \left( s+iku\right) p(k,s)-p_{0}(k)\right] .  \label{O}$$In the subballistic superdiffusive case, one can obtain for small $s$ $$\hat{\psi}\left( s\right) \simeq 1-Ts+ATs^{\mu },\qquad 1<\mu <2  \label{A}$$for which the first moment $T=\int_{0}^{\infty }\Psi (\tau )d\tau $ is finite and the second moment is divergent. Then $$\hat{K}(s)\simeq \frac{1}{T}\left( 1+As^{\mu -1}\right)$$as $s\rightarrow 0.$ From (\[main\]) we obtain the following equation $$T\frac{\partial ^{2}p}{\partial t^{2}}+\frac{\partial p}{\partial t}-D\frac{%
\partial ^{2}p}{\partial x^{2}}+A\int_{V}\varphi (u)L_{u}^{1-\mu }pdu=0
\label{le}$$with the diffusion coefficient $D=Tv^{2}.$ It is clear from (\[lee\]) and (\[le\]) that ballistic and subballistic cases are fundamentally different. In the strong ballistic case ($0<\mu <1$), the integral term is in the balance with classical wave equation, while for the subballistic superdiffusive case ($1<\mu <2$), the memory term is in the balance with the Cattaneo (telegraph) equation. One can perform various asymptotic analysis of (\[lee\]) and (\[le\]), obtain the pseudo-differential equations for the walker’s PDF position [@Bae; @Mer1; @Sib; @Mag] and determine the shape of PDF profiles [@Bar].

As an illustration let us find the long-time asymptotic solution to Eq. ([lee]{}). In the limit $\tau _{0}\rightarrow 0$, the evolution of the PDF $%
p\left( x,t\right) $ is determined by the integral term, while the first two wave equation terms in (\[lee\]) are irrelevant. The PDF $p\left(
x,t\right) $ obeys $$\int_{V}\varphi (u)L_{1-\mu }^{u}pdu=0.  \label{La}$$By using (\[O\]), we obtain $$\left( s-ikv\right) ^{1-\mu }\left[ \left( s+ikv\right) \hat{p}(k,s)-\hat{p}%
_{0}(k)\right] +$$$$\left( s+ikv\right) ^{1-\mu }\left[ \left( s-ikv\right) \hat{p}(k,s)-\hat{p}%
_{0}(k)\right] =0.  \label{LL}$$In particular, for $\mu =1/2$ and $\hat{p}_{0}(k)=1$, we factorize Eq. ([LL]{}) and find $$\left( s^{2}+k^{2}v^{2}\right) ^{\frac{1}{2}}\hat{p}(k,s)-1=0.$$By using the inverse Fourier-Laplace transform we obtain the well-known self-similar profile [@KlSo] $$p\left( x,t\right) =\pi ^{-1}\left( v^{2}t^{2}-x^{2}\right) ^{-\frac{1}{2}}.$$For the arbitrary $\mu $ from the interval $0<\mu <1,$ one can find from (\[LL\]) that the solution to (\[La\]) is the Lamperti distribution [Za]{}. It would be interesting to use a new wave equation for the non-normalizable density problem for superdiffusive anomalous transport [Re]{}. Note that our equation can be generalized for the case when the particle moves with the random velocity: $\dot{x}_{\pm }(t)=\pm v+\dot{B}%
\left( t\right) ,$ where $\dot{B}\left( t\right) $ is the Gaussian white noise.

In summary, we derived the integro-differential wave equation for the probability density function for a position of a Lévy walker with continuous sample paths. This equation involves a classical wave operator together with memory integrals describing a spatio-temporal coupling of Lévy walk. It is valid for any running time PDF and it does not involve any long-time large-scale approximations. It generalizes the well-known telegraph equation obtained from the persistent random walk. For the Lévy walk the integro-differential wave equation (\[main\]) plays the same role as the master equation for the uncoupled CTRW [@KlSo; @MFH]. Our technique may lead to the significant advances in the extension of the linear Lévy walk models. Our approach may be particularly helpful when we deal with the *Lévy particles* interactions [@Int] when it would be difficult to take into account nonlinear terms within the standard approach involving integral equations (see the similar problems for subdiffusion [@F2; @F1]). The description of the *Lévy walk* in terms of the Markovian balance equations (\[str\]) proves to be very useful for the analysis of *Lévy walk* with a random death process [@FedotAbby]. Our theory can be useful to formulate the impact of the external force or the chemotactic substance on the random movement of particles with finite velocities. This theory can be also useful for the implementation of the non-linear reactions and development of the theory of wave propagation in reaction-transport systems involving enhanced diffusion and memory effects [@MFH; @Memory].

*Acknowledgement*. This work was funded by EPSRC Grant No. EP/J019526/1 *Anomalous reaction-transport equations*. The author gratefully acknowledges the hospitality of the Isaac Newton Institute, Cambridge and thanks Nickolay Korabel and Vicenc Mendez for very useful discussions during the programme *Coupling Geometric PDEs with Physics for Cell Morphology, Motility and Pattern Formation.* The author thanks Eli Barkai for useful suggestions.

[99]{} J. Klafter, and I. M. Sokolov, *First Steps in Random Walks. From Tools to Applications* (Oxford University Press, Oxford, 2011).

P. Barthelemy, J. Bertolotti, and D. S. Wiersma, Nature **453**, 495 (2008).

D. A. Kessler and E. Barkai, Phys. Rev. Lett. **108**, 230602 (2012).

T. H. Harris et all, Nature **486**, 545 (2012).

K. Chen, B. Wang and S. Granick. Nature Mater. **14,** 589 (2015).

G. M. Viswanathan et all, Nature **401**, 911 (1999); N. E. Humphries, et al. Nature **465**, 1066 (2010); M. de Jager et all, Science **332**, 1551 (2011).

V. Mendez, D. Campos, F. Bartumeus, *Stochastic foundations in movement ecology: anomalous diffusion, front propagation and random searches* (Springer Science & Business Media, 2013)

V. Zaburdaev, S. Denisov, and J. Klafter, Rev. Mod. Phys. **87**, 483 (2015).

J. Klafter, A. Blumen, and M. F. Shlesinger, Phys. Rev. A **35**, 3081 (1987).

A. Blumen, G. Zumofen, J. Klafter, Phys. Rev. A, **40,** 3964 (1989).

V. Yu. Zaburdaev and K. V. Chukbar, JETP **94,** 252 (2002).

I. M. Sokolov and R. Metzler, Phys. Rev. E **67**, 010101(R) (2003).

R. Metzler, Phys. Rev. E, **62**, 6233 (2000).

E. Barkai and R. Silbey, J. Phys. Chem. B **104,** 3866 (2000).

R. Friedrich, F. Jenko, A. Baule, S. Eule, Phys. Rev. Lett. **96**, 230601 (2006).

V. Mendez, S. Fedotov, and W. Horsthemke, *Reaction-Transport Systems* (Springer, Berlin, 2010).

D. R. Cox and H. D. Miller, The Theory of Stochastic Processes (Methuen, London, 1965).

M. O. Vlad, V. T. Popa, and E. Segal, Phys. Lett. A **100**, 387 (1984); M. O. Vlad and J. Ross, Phys. Rev. E **66**, 061908 (2002).

R. Ferrari, A. J. Manfroi, and W. R. Young. Physica D **154,** 111 (2001).

S. Fedotov, A. Tan, and A. Zubarev, Phys. Rev. E **91**, 042124 (2015).

J. Masoliver, K. Lindenberg, and G. H. Weiss, Physica A **157**, 891 (1989).

G. Zumofen, and J. Klafter, Phys. Rev. E, **47**, 851 (1993).

E. Korobkova et al., Nature (London) **428**, 574 (2004).

M. Kac. Rocky Mountain J. Math. **4**, 497 (1974).

A. Di Crescenzo. Advances in Applied Probability, **33** (2001) 690.

A. Di Crescenzo and B. Martinucci, Computer Aided Systems Theory–EUROCAST 2007, Lecture Notes in Computer Science, Vol. 4739 (Springer, Berlin, Heidelberg, 2007), pp. 163–170.

E. Scalas, R. Gorenflo, F. Mainardi, and M. Raberto, Fractals **11**, 281 (2003).

S. Fedotov and S. Falconer, Phys. Rev. E **87**, 052139 (2013).

M. M. Meerschaert, D. A. Benson, H. P. Scheffler, and P. Becker-Kern. Phys. Rev. E, **66**, 102R (2002).

B. Baeumer, M. M. Meerschaert, J. Mortensen, Proc. Amer. Math. Soc., **133**, 2273 (2005).

P. Becker-Kern, M. M. Meerschaert and H. P. Scheffler. Ann. Probab. **32**, 730 (2004).

V. V. Uchaikin, and R. T. Sibatov, J. Phys. A: Math. Theor. **44,** 145501 (2011).

M. Magdziarza, H.P. Scheffler, P. Straka, P. Zebrowski, Stoch. Proc. and their Appl. (2015, in press).

D. Froemberg, M. Schmiedeberg, E. Barkai, and V. Zaburdaev, Phys. Rev. E **91**, 022131 (2015).

A. Rebenshtok, S. Denisov, P. Hanggi, E. Barkai, Phys. Rev. Lett, **112,** 110601 (2014).

M. de Jager et all, Proc. R. Soc. B** 281,** 20132605 (2014).

S. Fedotov, Phys. Rev. E **88**, 032104 (2013).

P. Straka and S. Fedotov, J. Theor. Biology** 366**, 71 (2015).

A. Yadav, S. Fedotov, V. Mendez, and W. Horsthemke, Phys. Lett. A **371**, 374 (2007).
