arxiv: v1 [math.ap] 22 Apr 2015
|
|
- Maximilian Shelton
- 5 years ago
- Views:
Transcription
1 Densities of Lévy walks and the corresponding fractional equations Marcin Magdziarz, Tomasz Zorawik arxiv: v [math.ap] Apr 5 Abstract In this paper we derive eplicit formulas for the densities of Lévy walks. Our results cover both jump-first and wait-first scenarios. The obtained densities solve certain fractional differential equations involving fractional material derivative operators. In the particular case, when the stability inde is rational, the densities can be represented as an integral of Meijer G function. This allows to efficiently evaluate them numerically. Our results show perfect agreement with the Monte Carlo simulations. Introduction The Continuous Time Random Walk CTRW is a stochastic process determined uniquely by i.i.d. random variables J, J,... representing consecutive jumps of the random walker, and i.i.d. positive random variables T,T,... representing waiting times between jumps, []. The trajectories of CTRW are step functions with intervals T i and jumps J i. Lévy walk is a particular case of CTRW satisfying the following additional condition J i = T i for every i N length of the jump equal to the length of the preceding waiting time. Lévy walks were defined for the first time in [, ] in the framework of generalized master equations. Since then, they became one of the most popular models in statistical physics with large number of important applications. Lévy walks have been found to be ecellent models in the description of various real-life phenomena and comple anomalous systems. The most striking eamples are: transport of light in optical materials [], foraging patterns of animals [, 3, 6], epidemic spreading [4, 7], human travel [5, 8], blinking nano-crystals [6], and fluid flow in a rotating annulus [4]. The main idea underlying the Lévy walk is the spatial-temporal coupling, which is manifested by the condition J i = T i. Thus, even if we assume that the distribution of jumps Y i is heavy-tailed with diverging moments, the Lévy walk itself has finite moments of all orders intuitively, long jumps are penalized by requiring more time to be performed. This is very different from α-stable Lévy processes with α <, which have infinite second moment. These desirable properties make Lévy walk particularly attractive for physical applications. Let us now recall the formal definition of Lévy walk. Let T i, i =,,..., be the sequence of i.i.d. positive random variables, representing the waiting times of the walker. Assume that they belong to the domain of attraction of one-sided α-stable law, i.e. [nt] n /α i= T i d S α t as n. Here, S α t is the α-stable subordinator with the Laplace transform [9] Eep{ ss α t} = ep{ ts α }, < α <.
2 Denote by Nt = ma{k : k T i t} the corresponding counting process. Net, define the jumps of the walker as J i = I i T i, where I,I,... are i.i.d. random variables, which are assumed independent of the sequence of waiting times T i. Each I i governs the direction of the jump, i.e. i= PI i = = p, PI i = = p, p. Thus, the walker jumps up with probability p and down with probability p. Clearly, the condition J i = T i is satisfied. The sequence of jumps J i defined above belongs to the domain of attraction of α-stable distribution [3] [nt] n /α i= J i d L α t as n. The process L α t is the α-stable Lévy motion with the corresponding Fourier transform Eep{ikL α t} = ep{ t k α cosπα/ ip tanπα/sgnk}. Later on we will also use the corresponding left-limit process Finally, the process L def α t = liml α s. sրt Nt J i Rt = is called Lévy walk. Rt is also known as wait-first Lévy walk in the literature [5], since the walker at the beginning of its motion t = first waits and then performs the jump. As shown in [3], Rt obeys the following scaling limit in distribution Rnt n i= d Xt. as n. Here Xt is the right-continuous version of the process L αsα t. Moreover, Sα t is the inverse of S α t, i.e. Sα t = inf{τ : S α τ > t}. Additionally, the instants of jumps as well as the respective jump lengths of the processes L α t and S α t are eactly the same. Also, the probability density function PDF p t of Xt satisfies the following fractional equation [, 3] [ p t α + p t + α ] t α p t = δ Γ α..3 α t y are the fractional material derivatives introduced in [3]. In the Here, the operators Fourier-Laplace space they are given by { F y L t t y α } fy,t = s ik α fk,s.
3 In what follows, we will find the eplicit solution of.3 by determining the PDF p t of the limit process Xt given in.. We will also consider the so-called jump-first Lévy walk [5] Rt = Nt+ i= It has the following scaling limit in distribution [4, 5] Rnt n J i. d Yt = L α S α t.4 as n. Moreover, the PDF w t y of Yt satisfies the fractional equation of the form [, 4] [ p t α + p y t + α ] w t y = y α Γ α t δ y uu α du..5 In the net two sections we will derive eplicit formulas for the PDFs of the limit processes Xt and Yt, respectively. This way we will obtain solutions of fractional equations.3 and.5. Densities of wait-first Lévy walks Let us start with the wait-first scenario. In the result below we determine the PDF of the limit processes Xt from eq... Theorem. If p, and α, then the PDF p t of the process Xt = L α S α t equals p t = α pα p α Γ α t z α w + w t w αr α z r p α z dzdw α p,t, α. where r α is a density of a positive α-stable random variable Z α with the Laplace transform Ee uzα = e uα. Proof. We start with reminding the following equality of distributions see [3]: L α t,s α t d = p α S α t p α S α t,p α S α t+ p α S α t,. where S α t and S α t are independent copies of S α t. Using this fact it was shown in [3] that the Lévy measure ν Lα,S α of L α t,s α t equals νd,ds = pδ s dν Sα ds+ pδ s dν Sα ds,.3 3
4 where ν Sα ds = α Γ α α ds. Hence, an infinitesimal generator A of the process L α t,s α t equals Af,s = [f+y,s+w f,s]ν Lα,S αdy,dw R = [f+y,s+w f,s].4 R pδ w dyν Sα dw+ pδ w dyν Sα dw. In [7] the authors consider among other processes the -dimensional process Xt, Vt, where Vt is the age process counting time that passed since the last jump ofxt. Equation.4 together with Remark 4. in the mentioned paper provides us with the joint distribution of Xt,Vt PXt = d,vt = dv = ν Lα,S αr [v, Ud,t dv { v t},.5 where Ud,ds is the -potential measure of the process L α t,s α t, defined as Ud,ds = PL α u = d,s α u = dsdu. The processvt is beyond our interest, but we will later obtain the distribution ofxt as a marginal distribution of Xt,Vt. This is the reason why we turn now our attention to Eq..5. We have ν Lα,S αr [v, = pν Sα [v, + pν Sα [v, = ν Sα [v, α = Γ α α d = v α Γ α..6 v Furthermore, a density u,s of the potential measure U can be calculated as u,s = w u,sdu,.7 where w t,s is the PDF of the process L α t,s α t. From Eq.. we get PL α t = d,s α t = ds = P S α = ds+d,s p α = ds d. α p α The Jacobian determinant in the above formula for the linear transformation of d, ds equals. Taking into account the independence of S p /α p /α α t and S α t we epress w t,s in terms of q t s - the PDF of S α t: s+ s w t,s = q pα p t q α p t. α p α Moreover the self-similarity of S α t provides us with the equation q u = q u = r α u α u α, α u α 4
5 here r α is the density of the random variable S α. We take into account the two above equations in Eq..7 and substitute u = z α to calculate the integral: s+ s u,s = q pα p u q α p u du α p α = u s+ pα p α rα u s α p α r α u α p α du.8 α α = z α s+ s r pα p α z r α p α z dz. α p α Finally we combine eqs..5,.6,.8 and integrate with respect to dv to obtain α p t = Γ α pα p α t v α z α t v + t v r α z r p α z dzdv α p α α = Γ α pα p α t z α w + w t w αr α z r p α z dzdw α p α We used here the fact that r α vanishes outside the positive half-line and substituted v = t w. One should mention here that the PDF of wait-first Lévy walk in the etreme cases p = and p = was already derived in [3]. For the special case α = we get the following simple epression for p t: Corollary. When α = and p, the PDF p t can be epressed as p t = π p p t p t+ pt+p + p+,t Proof. We have see [] r / = π 3 ep 4 Taking into account this formula and using the property of Gamma function Γ3/ = Γ/ = π we can calculate the following integral: z w + w r / p z r / p z dz = 4π w +w 3. 4p p z 5 ep p w + p w + w +w = 3 p 3 p 3. π p w + p w z dz.9
6 Substituting Eq..9 into Eq.. we arrive at p t = π p p t t w p +p w+ p p 3 dw,t. One can notice that below a,b and c are constants d c w = c w dw a+bca+bw a+bw 3, which finally implies p t = π p p t p t+ pt+p + p+,t. It is worth to mention here, that in a special case α = and p = we recover the known result [3] p t = t t w w 3 dw π,t = t πt,t. Figures and present the densities obtained from the above Corollary. 3.5 t=5 t= t= p t Figure : Plot of PDF p t of the process L αs α t for different values of t, α =.5, p =.. The drawback of using Theorem. to calculate p t in the general case is the necessity of knowing values of r α. One can use some known algorithms to approimate r α, however results are not perfect. Computations are time-consuming and inaccurate. They differ from those obtained via Monte-Carlo methods. However, as the net Corollary shows, we can epress p t in the form of an integral from Meijer G function see [9]. This special function is implemented in most of numerical packages, including Mathematica and Matlab. This representation is valid for rational α and can be used to approimate irrational α with a rational one. 6
7 6 5 p= p= p= Figure : For α =.5 we compare the obtained densities p for different values of p red solid lines with densities estimated using Monte Carlo methods blue pluses. Corollary. If p, and α = l k, where l,k N, then the PDF p t of the process Xt = L α S α t equals p t = α k+l k [ α t w + α pπ k l l α Γ α t w α w ] G k,k p α w+ l k,k α l, l, l+k,l+k p w k,, l,l α dw,t, l where G k,k l,l k,a l,b is a special list of k elements. is the Meijer G function see [9] and k,a = a k, a+ k a+ a+k,a+ k,..., k Proof. For α = l/k where k, l are positive integers with k > l we can epress r α in terms of the Meijer G function see [8]: r l/k = kl l l π k l/ Gk, l,k k k l l, k,. 7
8 This gives us z α w + w q α z q p α z dz α p α kl 4 p αp α l = z α π k l G k, ll w + l, w w + l,k z k k p k, α l G k, ll w l, l,k z k k p dz k, α k 4 p αp l α = u α π k l l G k, ll w+ w w + l,k u l, k k p k, α l G k, ll w l,k u l, k k p du k, α k 4 p αp α l l α/l α w + = π k l w w + k k p α G k,k p α w + l k,k α l, l, l+k,l+k p w k,, l,l α. l We substituted in the above equations z l = u. Moreover we used properties of Meijer G functions see [9] to calculate the integral from multiplication of these functions. To end the proof we substitute the calculated above integral into Eq... Figure 3 presents densities calculated here. Another method for calculating r α, where is not necessarily rational, was presented in []. α=.5 α=.5 α= Figure 3: For p =.5 we compare the densities p obtained in Corollary. for different values of α red solid lines with densities estimated using Monte Carlo methods blue pluses. 3 Densities of first-jump Lévy walks This section is devoted to the first-jump scenario. Below we find the PDF of the process Yt from.4. At the same time this PDF solves the fractional equation.5. 8
9 Theorem 3. If p, and α, then the PDF w t y of the process Yt = L α S α t equals: i if y < t, then ii if y t, then y t α p α w t y = z α p αγ α t t y+ w + w r α z r p α z y α dzdwd α p α + α p α p αγ α t t z α y t+y w + w r α z r p α z y α dzdwd, α p α w t y = α p α t t z α p αγ α t w + w r α z r p α z y α dzdwd, α p α iii if y t, then w t y = α p α p αγ α t t z α t w + w r α z r p α z y α dzdwd. α p α Here r α y is a density of a positive α-stable random variable Z α with the Laplace transform and a b = maa,b Proof. From Remark 4. in [7] we get PYt = dy,rt = dr = Ee uzα = e uα R w [,t] Ud,dwν Lα,S α dy,dr +t w where Rt is a process that counts how long the process Yt remains constant from the moment t. After applying Eq..3 the above equation reads PYt = dy,rt = dr [ = Ud,dw pδ dr+t w dy ν Sα dr +t w R w [,t] ] + pδ dr t+w dy ν Sα dr +t w 9
10 Now we can recover the distribution of Yt as a marginal distribution. We calculate the density w t y applying Tonelli s theorem to change the order of integration: w t y = α Γ α r [, R w [,t] [ u, w pδ r+t w y dr +t w α + pδ dr t+w y dr +t w α] dwddr α t t [ = u, w p Γ α,y t y+, wy α t + p y, t+y, w y α] dwd. The area of integration in the above integral depends on the relation of y with t - we have three cases. The first case is described by the condition y < t. Then If y t then w t y = w t y = and finally if y t then pα y t Γ α t + pα Γ α t y+ t t y t+y u,wy α dwd u,w y α dwd. pα y t u,wy α dwd Γ α t t y+ w t y = pα Γ α t t y t+y u,w y α dwd. Substituting Eq..8 for u,w in all three cases ends the proof. In the special case α = we get the following simple epression for w ty: Corollary 3. When α = the PDF w ty can be epressed as if y < t and as w t y = p pt π w t y = p π t t y t+y +y t y +t+y yt y t+y p pt+ py t y py t + py +t if y t and as w t y = p π t y pt y + p y t if y t.
11 Proof. The density w t y is symmetric in a following sense: w t y p = w t y p, where w t y c denotes the PDF of the process Yt with parameter p = c. Thus it is enough to calculate w t y for y > - we assume now y >. There are two cases: y < t and y t. In the first case from Theorem 3. and Eq..9 we have w t y = p p 3/ π + p p 3/ π y t t t y+ p w y+ p w +y 3/ y 3/ dwd t t y t+y p w y+ p w +y 3/ y 3/ dwd Notice that in the above integrals t y+ = t y+ when y t and t y+ = = if y t y+t. Similarly t+y = t+y when and t+y = = if y+t. Thus w ty can be epressed as a sum of four integrals, each of them having a simple region of integration. We can calculate all of them using standard integration techniques. After tedious computations we finally get the desired result for w t y in case when y,t. In the second case, that is y t, from Theorem 3. and Eq..9 we have w t y = p p 3/ π t p w y+ p w +y 3/ y 3/ dwd + t t t p w y+ p w +y 3/ y 3/ dwd Again, applying standard integration techniques we obtain the desired result for w t y when y > t. In a special case α = and p = the PDF w ty has the form if y < t and w t y = t π t t+y t y y t y +t+y w t y = π yt y t+y t t y y t + y +t if y t. Figure 4 presents w y for different values of p. Notice that in all cases we have a characteristic sharp peak at y = t and y = t. However when p the left peak gets bigger and the right one diminishes. The opposite situation appears when p.
12 p= p= p=.5 y y.4 y Figure 4: For α =.5 we compare the densities w y obtained in Corollary 3. for different values of p red solid lines with densities estimated using Monte Carlo methods blue pluses..5. t=5 t= t=.5 w t y y Figure 5: Plot of PDF w t y of the process L α S α t for different values of t, α =.5, p =.5. Acknowledgements This research was partially supported by NCN Maestro grant no. /6/A/ST/58. References [] P. Barthelemy, P.J. Bertolotti, D.S. Wiersma, A Lévy flight for light. Nature 453 8, [] W.J. Bell Searching Behaviour. Chapman & Hall, London 999. [3] H.C. Berg Random Walks in Biology. Princeton University Press, Princeton 983. [4] D. Brockmann, Human mobility and spatial disease dynamics. Chap., pp. 4 in Reviews of Nonlinear Dynamics and Compleity, Vol., edited by H. Shuster, Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim. [5] D. Brockmann, L. Hufnagel, T. Geisel, The scaling laws of human travel. Nature 439 6, [6] M. Buchanan, Ecological modelling: The mathematical mirror to animal nature. Nature 453 8,
13 [7] B. Dybiec, Random strategies of contact tracking. Physica A 387 8, [8] M.C. Gonzales, C.A. Hidalgo, A.L. Barabási, Understanding individual human mobility patterns. Nature 453 8, [9] A. Janicki, A. Weron, Simulation and Chaotic Behavior of α-stable Stochatic Processes. Dekker, New York 994. [] A. Jurlewicz, P. Kern, M.M. Meerschaert, H.P Scheffler, Fractional governing equations for coupled random walks. Comput. Math. Appl. 64, 3Ű-336. [] J. Klafter, A. Blumen, M.F. Shlesinger, Stochastic pathway to anomalous diffusion. Phys. Rev. A , [] J. Klafter, I.M. Sokolov, First Steps in Random Walks. From Tools to Applications. Oford University Press, Oford. [3] M. Magdziarz, W. Szczotka, P. Zebrowski, Langevin Picture of Lévy Walks and Their Etensions. J. Stat. Phys. 7,, [4] M. Magdziarz, M. Teuerle, Asymptotic properties and numerical simulation of multidimensional Lévy walks. Commun. Nonlinear Sci. Numer. Simul. 5, [5] M. Magdziarz, M. Teuerle, P. Zebrowski, Scaling limits of overshooting Levy walks, Acta Phys. Polon. B 43, 3. [6] G. Margolin, E. Barkai, Nonergodicity of blinking nanocrystals and other Lévy-walk processes. Phys. Rev. Lett. 94 5, 86. [7] M.M. Meerchaert, P. Straka: Semi-Markov approach to continuous time random walk limit processes. The Annals of Probability, Vol. 4, No. 4 4, [8] K.A. Penson, K. Gorska, Eact and Eplicit Probability Densities for One-Sided Lévy Stable Distributions. Phys. Rev. Lett. 5, 64. [9] A.P. Prudnikov, Yu. A. Brychkov A., O. I. Marichev O.I., Integrals and Series. Gordon and Breach, Amsterdam 998. [] A. Saa, R. Venegeroles, Alternative numerical computation of one-sided LÃľvy and Mittag- Leffler distributions, Phys.Rev.E 84, 67. [] K.I. Sato, Lévy Processes and Infinitely Divisible Distributions. Cambridge University Press, Cambridge 999. [] M.F. Shlesinger, J. Klafter, Y.M Wong, Random walks with infinite spatial and temporal moments, J. Stat. Phys. 7 98, [3] I.M. Sokolov, R. Metzler, Towards deterministic equations for Lévy walks: The fractional material derivative. Phys. Rev. E 67 3, R. [4] T.H. Solomon, E.R. Weeks, H.L. Swinney, Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow. Phys. Rev. Lett , [5] V. Zaburdaev, S. Denisov, J. Klafter, Lévy Walks. arxiv:
14 Hugo Steinhaus Center, Department of Mathematics, Wroclaw University of Technology, Wyspianskiego 7, 5-37 Wroclaw, Poland. 4
arxiv: v1 [math.pr] 25 May 2018
Scaling limits for Lévy walks with rests arxiv:1805.10027v1 [math.pr] 25 May 2018 Abstract Marek Teuerle Faculty of Pure and Applied Mathematics, Hugo Steinhaus Center, Wroc law University of Science and
More informationSUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION
Vol. 4 (29) ACTA PHYSICA POLONICA B No 5 SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION Joanna Janczura, Agnieszka Wyłomańska Institute of Mathematics and Computer Science Hugo Steinhaus
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationarxiv: v1 [cond-mat.stat-mech] 23 Apr 2010
Fractional Fokker-Planck Equations for Subdiffusion with Space-and-Time-Dependent Forces. B.. Henry Department of Applied Mathematics, School of Mathematics and Statistics, University of New South Wales,
More informationITÔ FORMULA FOR SUBORDINATED LANGEVIN EQUATION. A CASE OF TIME DEPENDENT FORCE
Vol. 4 (29) ACTA PHYSICA POLONICA B No 5 ITÔ FORMULA FOR SUBORDINATED LANGEVIN EQUATION. A CASE OF TIME DEPENDENT FORCE Alesander Weron Hugo Steinhaus Center, Institute of Mathematics and Computer Science
More informationRandom Walk Model with Waiting Times Depending on the Preceding Jump Length
Journal of Statistical Physics, Vol. 123, No. 4, May 2006 ( C 2006 ) DOI: 10.1007/s10955-006-9104-z Random Walk Model with Waiting Times Depending on the Preceding Jump Length Vasily Yu. Zaburdaev 1 Received
More informationWeak Ergodicity Breaking WCHAOS 2011
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,
More informationCovariance structure of continuous time random walk limit processes
Covariance structure of continuous time random walk limit processes Alla Sikorskii Department of Statistics and Probability Michigan State University Fractional Calculus, Probability and Non-local Operators:
More informationWeak Ergodicity Breaking. Manchester 2016
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Burov, Froemberg, Garini, Metzler PCCP 16 (44), 24128 (2014) Akimoto, Saito Manchester 2016 Outline Experiments: anomalous diffusion of single molecules
More informationTheory of fractional Lévy diffusion of cold atoms in optical lattices
Theory of fractional Lévy diffusion of cold atoms in optical lattices, Erez Aghion, David Kessler Bar-Ilan Univ. PRL, 108 230602 (2012) PRX, 4 011022 (2014) Fractional Calculus, Leibniz (1695) L Hospital:
More informationLecture 25: Large Steps and Long Waiting Times
Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,
More informationWeak Ergodicity Breaking
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t
More informationarxiv: v1 [q-bio.nc] 5 Dec 2015
Modelling axon growing using CTRW arxiv:1512.02603v1 [q-bio.nc] 5 Dec 2015 Xavier Descombes, Inria Sophia Antipolis, France Elena Zhizhina, Sergey Komech Institute for Information Transmission Problems,
More informationLecture notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT September 24, Anomalous diffusion
Lecture notes for 12.086/12.586, Modeling Environmental Complexity D. H. Rothman, MIT September 24, 2014 Contents 1 Anomalous diffusion 1 1.1 Beyond the central limit theorem................ 2 1.2 Large
More informationIntroduction: the plague. The Scaling Laws Of Human Travel. Introduction: SARS (severe acute respiratory syndrom)
The Scaling Laws Of Human Travel Verena Zuber Department of Statistics University of Munich 7. July 2006 Introduction: the plague highly infectious disease several pandemics have influenced the human history
More informationVolatility and Returns in Korean Futures Exchange Markets
Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,
More informationFull characterization of the fractional Poisson process
epl draft Full characterization of the fractional Poisson process Mauro Politi,3, Taisei Kaizoji and Enrico Scalas 2,3 SSRI & Department of Economics and Business, International Christian University, 3--2
More informationCONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES
CONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES Aleksei V. Chechkin, 1,2 Michael Hofmann 3, and Igor M. Sokolov 3 1 School of Chemistry, Tel Aviv University, Ramat Aviv, 69978 Tel Aviv, Israel
More informationSemi-Markov approach to continuous time random walk limit processes
Semi-Markov approach to continuous time random walk limit processes Joint Mathematics Meetings San Antonio, TX January 13, 2015 Mark M. Meerschaert Michigan State University East Lansing, Michigan, USA
More informationArea coverage of radial Lévy flights with periodic boundary conditions
PHYSICAL REVIEW E 87, 4236 (23) Area coverage of radial Lévy flights with periodic boundary conditions Mahsa Vahabi,,2 Johannes H. P. Schulz, Babak Shokri, 2,3 and Ralf Metzler 4,5 Physics Department,
More informationTime fractional Schrödinger equation
Time fractional Schrödinger equation Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 The Schrödinger equation is considered with the first order time
More informationCauchy Problems in Bounded Domains and Iterated Processes
Cauchy Problems in Bounded Domains and Iterated Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J
More informationAnomalous wave transport in one-dimensional systems with. Lévy disorder
Anomalous wave transport in one-dimensional systems with Lévy disorder Xujun Ma,2 and Azriel Z. Genack,2 Department of Physics, Queens College of the City University of New York, Flushing, NY, 367 USA
More informationCLUSTER CONTINUOUS TIME RANDOM WALKS
1 CLUSTER CONTINUOUS TIME RANDOM WALKS A. JURLEWICZ, M. M. MEERSCHAERT, AND H.-P. SCHEFFLER Abstract. In a continuous time random walk (CTRW), a random waiting time precedes each random jump. The CTRW
More informationCTRW Limits: Governing Equations and Fractal Dimensions
CTRW Limits: Governing Equations and Fractal Dimensions Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 19-23, 2013 Joint work with Z-Q. Chen, M. D Ovidio, M.M. Meerschaert,
More informationLévy Walks and scaling in quenched disordered media
Lévy Walks and scaling in quenched disordered media A. Vezzani, - CNR Modena - Parma L. Caniparoli - Sissa Trieste S.Lepri - Firenze CNR - ISC Raffaella Burioni - Parma Sperlonga - September 2010 Work
More informationThe Portevin-Le Chatelier effect in the Continuous Time Random Walk. framework
The Portevin-Le Chatelier effect in the Continuous Time Random Walk framework A. Sarkar * and P. Barat Variable Energy Cyclotron Centre 1/AF Bidhan Nagar, Kolkata 764, India PACS numbers: 62.2.Fe, 5.4.Fb,
More informationCauchy Problems Solved by Running Subordinate Processes
Cauchy Problems Solved by Running Subordinate Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J 4
More informationCluster continuous time random walks
STUDIA MATHEMATICA 25 (1) (211) Cluster continuous time random walks by Agnieszka Jurlewicz (Wrocław), Mark M. Meerschaert (East Lansing, MI) and Hans-Peter Scheffler (Siegen) Abstract. In a continuous
More informationA short story about FRAP (fluourescence recovery after photobleaching)
1 A short story about FRAP (fluourescence recovery after photobleaching) Bob Pego (Carnegie Mellon) work involving team leader: James McNally (NCI Lab of Receptor Biology & Gene Expression) Biophysicists:
More informationEvaluation of integrals by differentiation with respect to a parameter
December 8 Evaluation of integrals by differentiation with respect to a parameter Khristo N Boyadzhiev Department of Mathematics and Statistics, Ohio Northern University, Ada, OH 458, USA k-boyadzhiev@onuedu
More informationPERSISTENT RANDOM WALK EFFECT IN THE SUBDIFFUSION-REACTION PROCESS
Vol. 45 (214) ACTA PHYSICA POLONICA B No 9 PERSISTENT RANDOM WALK EFFECT IN THE SUBDIFFUSION-REACTION PROCESS Tadeusz Kosztołowicz, Mateusz Piwnik Institute of Physics, Jan Kochanowski University Świętokrzyska
More informationAnomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm
Anomalous Lévy diffusion: From the flight of an albatross to optical lattices Eric Lutz Abteilung für Quantenphysik, Universität Ulm Outline 1 Lévy distributions Broad distributions Central limit theorem
More informationarxiv: v1 [cond-mat.stat-mech] 15 Jan 2012
arxiv:2.32v [cond-mat.stat-mech] 5 Jan 22 The nonextensive entropy approach versus the stochastic in describing subdiffusion Tadeusz Koszto lowicz Institute of Physics, Jan Kochanowski University, ul.
More informationA Stochastic Paradox for Reflected Brownian Motion?
Proceedings of the 9th International Symposium on Mathematical Theory of Networks and Systems MTNS 2 9 July, 2 udapest, Hungary A Stochastic Parado for Reflected rownian Motion? Erik I. Verriest Abstract
More informationAnomalous diffusion of volcanic earthquakes
Anomalous diffusion of volcanic earthquakes SUMIYOSHI ABE 1,2 and NORIKAZU SUZUKI 3 1 Department of Physical Engineering, Mie University, Mie 514-8507, Japan 2 Institute of Physics, Kazan Federal University,
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
More informationAnomalous diffusion in cells: experimental data and modelling
Anomalous diffusion in cells: experimental data and modelling Hugues BERRY INRIA Lyon, France http://www.inrialpes.fr/berry/ CIMPA School Mathematical models in biology and medicine H. BERRY - Anomalous
More informationarxiv: v1 [cond-mat.stat-mech] 15 Nov 2018
arxiv:1811.06265v1 [cond-mat.stat-mech] 15 Nov 2018 Multimodal stationary states in symmetric single-well potentials driven by Cauchy noise Karol Capała, Bartłomiej Dybiec Marian Smoluchowski Institute
More informationNoisy Lévy walk analog of two-dimensional DNA walks for chromosomes of S. cerevisiae
PHYSICAL REVIEW E VOLUME 58, NUMBER 1 JULY 1998 Noisy Lévy walk analog of two-dimensional DNA walks for chromosomes of S. cerevisiae Guillermo Abramson, 1, * Pablo A. Alemany, 1,2 and Hilda A. Cerdeira
More informationA SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION
THERMAL SCIENCE: Year 7, Vol., Suppl., pp. S-S7 S A SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION by Wei XU a, Wen CHEN a*, Ying-Jie LIANG a*, and Jose WEBERSZPIL b a State Key Laboratory
More informationarxiv: v2 [cond-mat.stat-mech] 12 Apr 2011
epl draft arxiv:8.762v2 cond-mat.stat-mech] 2 Apr 2 Raman Research Institute, Bangalore 568, India PACS 2.5.-r Probability theory, stochastic processes, and statistics PACS 5.4.-a Fluctuation phenomena,
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998
Multiplicative processes and power laws arxiv:cond-mat/978231v2 [cond-mat.stat-mech] 23 Feb 1998 Didier Sornette 1,2 1 Laboratoire de Physique de la Matière Condensée, CNRS UMR6632 Université des Sciences,
More informationNormal and Anomalous Diffusion: A Tutorial
Normal and Anomalous Diffusion: A Tutorial Loukas Vlahos, Heinz Isliker Department of Physics, University of Thessaloniki, 54124 Thessaloniki, Greece Yannis Kominis, and Kyriakos Hizanidis School of Electrical
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationCoupled continuous time random walks in finance
Physica A 37 (26) 114 118 www.elsevier.com/locate/physa Coupled continuous time random walks in finance Mark M. Meerschaert a,,1, Enrico Scalas b,2 a Department of Mathematics and Statistics, University
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationCover times of random searches
Cover times of random searches by Chupeau et al. NATURE PHYSICS www.nature.com/naturephysics 1 I. DEFINITIONS AND DERIVATION OF THE COVER TIME DISTRIBUTION We consider a random walker moving on a network
More informationTransform Techniques - CF
Transform Techniques - CF [eview] Moment Generating Function For a real t, the MGF of the random variable is t t M () t E[ e ] e Characteristic Function (CF) k t k For a real ω, the characteristic function
More informationLARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*
LARGE EVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILE EPENENT RANOM VECTORS* Adam Jakubowski Alexander V. Nagaev Alexander Zaigraev Nicholas Copernicus University Faculty of Mathematics and Computer Science
More informationA Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions
Int. J. Contemp. Math. Sci., Vol. 1, 2006, no. 4, 155-161 A Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions Tomasz J. Kozubowski 1 Department of Mathematics &
More informationRandom walk through Anomalous processes: some applications of Fractional calculus
Random walk through Anomalous processes: some applications of Fractional calculus Nickolay Korabel Fractional Calculus Day @ UCMerced 12 June 2013 Fractional Calculus Mathematics Engineering Physics Levy
More informationRandom walks in a 1D Levy random environment
Random Combinatorial Structures and Statistical Mechanics, University of Warwick, Venice, 6-10 May, 2013 Random walks in a 1D Levy random environment Alessandra Bianchi Mathematics Department, University
More informationFractional Cauchy Problems for time-changed Processes
Fractional Cauchy Problems for time-changed Processes Erkan Nane Department of Mathematics and Statistics Auburn University February 17, 211 Erkan Nane (Auburn University) Fractional Cauchy Problems February
More informationContinuous time random walks in finite domains and general boundary conditions: some formal considerations
IOP PUBISHING JOURNA OF PHYSICS A: MATHEMATICA AND THEORETICA J. Phys. A: Math. Theor. 41 (008) 15004 (11pp) doi:10.1088/1751-8113/41/1/15004 Continuous time random walks in finite domains and general
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationTransform Techniques - CF
Transform Techniques - CF [eview] Moment Generating Function For a real t, the MGF of the random variable is t e k p ( k) discrete t t k M () t E[ e ] e t e f d continuous Characteristic Function (CF)
More informationFRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS
FRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS L. Boyadjiev*, B. Al-Saqabi** Department of Mathematics, Faculty of Science, Kuwait University *E-mail: boyadjievl@yahoo.com **E-mail:
More informationPhysica A. Extremal behavior of a coupled continuous time random walk
Physica A 39 (211) 55 511 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa Extremal behavior of a coupled continuous time random walk Rina Schumer a,,
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationLévy-Brownian motion on finite intervals: Mean first passage time analysis
PHYSICAL REVIEW E 73, 046104 006 Lévy-Brownian motion on finite intervals: Mean first passage time analysis B. Dybiec* and E. Gudowska-Nowak Marian Smoluchowski Institute of Physics and Mark Kac Center
More informationWeak Limits for Multivariate Random Sums
Journal of Multivariate Analysis 67, 398413 (1998) Article No. MV981768 Weak Limits for Multivariate Random Sums Tomasz J. Kozubowski The University of Tennessee at Chattanooga and Anna K. Panorska - The
More informationWeak chaos, infinite ergodic theory, and anomalous diffusion
Weak chaos, infinite ergodic theory, and anomalous diffusion Rainer Klages Queen Mary University of London, School of Mathematical Sciences Marseille, CCT11, 24 May 2011 Weak chaos, infinite ergodic theory,
More informationContinuous Time Random Walk Limits in Bounded Domains
Continuous Time Random Walk Limits in Bounded Domains Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY May 16, 2012 Acknowledgement Boris Baeumer, University of Otago, New Zealand.
More informationThe story of the film so far... Mathematics for Informatics 4a. Continuous-time Markov processes. Counting processes
The story of the film so far... Mathematics for Informatics 4a José Figueroa-O Farrill Lecture 19 28 March 2012 We have been studying stochastic processes; i.e., systems whose time evolution has an element
More informationTransform Techniques - CF
Transform Techniques - CF [eview] Moment Generating Function For a real t, the MGF of the random variable is t t M () t E[ e ] e Characteristic Function (CF) k t k For a real ω, the characteristic function
More informationarxiv: v2 [physics.flu-dyn] 31 Aug 2016
arxiv:608.0724v2 [physics.flu-dyn] 3 Aug 206 ANALYSIS OF FLUID VELOCITY VECTOR FIELD DIVERGENCE u IN FUNCTION OF VARIABLE FLUID DENSITY ρ( x,t) const AND CONDITIONS FOR VANISHING VISCOSITY OF COMPRESSIBLE
More informationSearch for Food of Birds, Fish and Insects
Search for Food of Birds, Fish and Insects Rainer Klages Max Planck Institute for the Physics of Complex Systems, Dresden Queen Mary University of London, School of Mathematical Sciences Diffusion Fundamentals
More informationReaction, Lévy flights, and quenched disorder
PHYSICAL REVIEW E, VOLUME 65, 1119 Reaction, Lévy flights, and quenched disorder Ligang Chen and Michael W. Deem Department of Chemical Engineering, University of California, Los Angeles, California 995
More informationTaylor Series and Asymptotic Expansions
Taylor Series and Asymptotic Epansions The importance of power series as a convenient representation, as an approimation tool, as a tool for solving differential equations and so on, is pretty obvious.
More informationAnomalous diffusion in biology: fractional Brownian motion, Lévy flights
Anomalous diffusion in biology: fractional Brownian motion, Lévy flights Jan Korbel Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Minisymposium on fundamental aspects behind cell systems
More informationReturn probability on a lattice
Return probability on a lattice Chris H. Rycroft October 24, 2006 To begin, we consider a basic example of a discrete first passage process. Consider an unbiased Bernoulli walk on the integers starting
More informationLecture 4a: Continuous-Time Markov Chain Models
Lecture 4a: Continuous-Time Markov Chain Models Continuous-time Markov chains are stochastic processes whose time is continuous, t [0, ), but the random variables are discrete. Prominent examples of continuous-time
More informationSpectra of Large Random Stochastic Matrices & Relaxation in Complex Systems
Spectra of Large Random Stochastic Matrices & Relaxation in Complex Systems Reimer Kühn Disordered Systems Group Department of Mathematics, King s College London Random Graphs and Random Processes, KCL
More informationTopic 5.1: Line Element and Scalar Line Integrals
Math 275 Notes Topic 5.1: Line Element and Scalar Line Integrals Textbook Section: 16.2 More Details on Line Elements (vector dr, and scalar ds): http://www.math.oregonstate.edu/bridgebook/book/math/drvec
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More informationSome properties of the moment estimator of shape parameter for the gamma distribution
Some properties of the moment estimator of shape parameter for the gamma distribution arxiv:08.3585v [math.st] 7 Aug 0 Piotr Nowak Mathematical Institute, University of Wrocław Pl. Grunwaldzki /4, 50-384
More informationFractional Quantum Mechanics and Lévy Path Integrals
arxiv:hep-ph/9910419v2 22 Oct 1999 Fractional Quantum Mechanics and Lévy Path Integrals Nikolai Laskin Isotrace Laboratory, University of Toronto 60 St. George Street, Toronto, ON M5S 1A7 Canada Abstract
More informationRELIABILITY FOR SOME BIVARIATE BETA DISTRIBUTIONS
RELIABILITY FOR SOME BIVARIATE BETA DISTRIBUTIONS SARALEES NADARAJAH Received 24 June 24 In the area of stress-strength models there has been a large amount of work as regards estimation of the reliability
More informationarxiv: v2 [cond-mat.stat-mech] 23 Jan 2015
Lévy walks V. Zaburdaev, 1, S. Denisov, 2, 3, 4, and J. Klafter 5, 1 Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, D-1187 Dresden, Germany 2 Institute of Physics, University
More informationGeneralized Post-Widder inversion formula with application to statistics
arxiv:5.9298v [math.st] 3 ov 25 Generalized Post-Widder inversion formula with application to statistics Denis Belomestny, Hilmar Mai 2, John Schoenmakers 3 August 24, 28 Abstract In this work we derive
More informationsuperdiffusion blocking Lagrangian Coh. Structures ocean currents internal waves and climate Coriolis force leads to barriers to transport (KAM torus)
Oceanic & Atmospheric Dynamics YESTERDAY Coriolis force leads to 2D flow (distances > ~100 km) jets and vortices barriers to transport (KAM torus) TODAY superdiffusion blocking Lagrangian Coh. Structures
More informationFirst-Passage Statistics of Extreme Values
First-Passage Statistics of Extreme Values Eli Ben-Naim Los Alamos National Laboratory with: Paul Krapivsky (Boston University) Nathan Lemons (Los Alamos) Pearson Miller (Yale, MIT) Talk, publications
More informationState Space and Hidden Markov Models
State Space and Hidden Markov Models Kunsch H.R. State Space and Hidden Markov Models. ETH- Zurich Zurich; Aliaksandr Hubin Oslo 2014 Contents 1. Introduction 2. Markov Chains 3. Hidden Markov and State
More informationA REFINED FACTORIZATION OF THE EXPONENTIAL LAW P. PATIE
A REFINED FACTORIZATION OF THE EXPONENTIAL LAW P. PATIE Abstract. Let ξ be a (possibly killed) subordinator with Laplace exponent φ and denote by I φ e ξs ds, the so-called exponential functional. Consider
More informationA source model for ISDN packet data traffic *
1 A source model for ISDN packet data traffic * Kavitha Chandra and Charles Thompson Center for Advanced Computation University of Massachusetts Lowell, Lowell MA 01854 * Proceedings of the 28th Annual
More informationDiffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking
Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking David Kessler Bar-Ilan Univ. E. Barkai (Bar-Ilan) A. Dechant (Augsburg) E. Lutz (Augsburg) PRL, 105 120602 (2010) arxiv:1105.5496
More informationHow does a diffusion coefficient depend on size and position of a hole?
How does a diffusion coefficient depend on size and position of a hole? G. Knight O. Georgiou 2 C.P. Dettmann 3 R. Klages Queen Mary University of London, School of Mathematical Sciences 2 Max-Planck-Institut
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationarxiv: v1 [math.ap] 26 Mar 2013
Analytic solutions of fractional differential equations by operational methods arxiv:134.156v1 [math.ap] 26 Mar 213 Roberto Garra 1 & Federico Polito 2 (1) Dipartimento di Scienze di Base e Applicate per
More informationQMUL, 10 Jul 2009 A Short Introduction to the Theory of Lévy Flights
QMUL, 10 Jul 2009 A Short Introduction to the Theory of Lévy Flights A. Chechkin School of Chemistry, Tel Aviv University, ISRAEL and Akhiezer Institute for Theoretical Physics, National Science Center
More informationJordan Journal of Mathematics and Statistics (JJMS) 7(4), 2014, pp
Jordan Journal of Mathematics and Statistics (JJMS) 7(4), 2014, pp.257-271 ON RELATIONS FOR THE MOMENT GENERATING FUNCTIONS FROM THE EXTENDED TYPE II GENERALIZED LOGISTIC DISTRIBUTION BASED ON K-TH UPPER
More informationNumerical evaluation of the upper critical dimension of percolation in scale-free networks
umerical evaluation of the upper critical dimension of percolation in scale-free networks Zhenhua Wu, 1 Cecilia Lagorio, 2 Lidia A. Braunstein, 1,2 Reuven Cohen, 3 Shlomo Havlin, 3 and H. Eugene Stanley
More informationA rational approximation of the arctangent function and a new approach in computing pi
arxiv:63.33v math.gm] 4 Mar 6 A rational approimation of the arctangent function and a new approach in computing pi S. M. Abrarov and B. M. Quine March 4, 6 Abstract We have shown recently that integration
More informationFractional Laplacian
Fractional Laplacian Grzegorz Karch 5ème Ecole de printemps EDP Non-linéaire Mathématiques et Interactions: modèles non locaux et applications Ecole Supérieure de Technologie d Essaouira, du 27 au 30 Avril
More informationStability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk
Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk ANSAPW University of Queensland 8-11 July, 2013 1 Outline (I) Fluid
More informationarxiv: v1 [math.pr] 30 Mar 2014
Binomial discrete time ruin probability with Parisian delay Irmina Czarna Zbigniew Palmowski Przemys law Świ atek October 8, 2018 arxiv:1403.7761v1 [math.pr] 30 Mar 2014 Abstract. In this paper we analyze
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Pacet The Group Velocity Phys 375 Overview and Motivation: Last time we looed at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More information