arxiv: v1 [cond-mat.stat-mech] 25 Aug 2016

Transcription

1 Influence of eernal poenials on heerogeneous diffusion processes Ryis Kazakevičius and Julius Ruseckas Insiue of Theoreical Physics and Asronomy, Vilnius Universiy, Saulėekio 3, LT-1 Vilnius, Lihuania arxiv: v1 [cond-ma.sa-mech] 5 Aug 16 In his paper we consider heerogeneous diffusion processes wih he power-law dependence of he diffusion coefficien on he posiion and invesigae he influence of eernal forces on he resuling anomalous diffusion. The heerogeneous diffusion processes can yield subdiffusion as well as superdiffusion, depending on he behavior of he diffusion coefficien. We assume ha no only he diffusion coefficien bu also he eernal force has a power-law dependence on he posiion. We obain analyic epressions for he ransiion probabiliy in wo cases: when he power-law eponen in he eernal force is equal o η 1, where η is he power-law eponen in he dependence of he diffusion coefficien on he posiion, and when he eernal force has a linear dependence on he posiion. We found ha he power-law eponen in he dependence of he mean square displacemen on ime does no depend on he eernal force, his force changes only he anomalous diffusion coefficien. In addiion, he eernal force having he power-law eponen differen from η 1 limis he ime inerval where he anomalous diffusion occurs. We epec ha he resuls obained in his paper may be relevan for a more complee undersanding of anomalous diffusion processes. PACS numbers: 5.4.-a,.5.-r, 5.1.Gg I. INTRODUCTION There are many sysems and processes where he ime dependence of he cenered second momen is no linear as in he classical Brownian moion. Such family of processes is called anomalous diffusion. In one dimension he anomalous diffusion is characerized by he power-law ime dependence of he mean square displacemen MSD [1] = K α α. 1 Here K α is he anomalous diffusion coefficien. When α 1, his ime dependence deviaes from he linear funcion of ime characerisic for he Brownian moion. If α < 1, he phenomenon is called subdiffusion. Occurrence of subdiffusion has been eperimenally observed, for eample, in he behavior of individual colloidal paricles in wodimensional random poenial energy landscapes []. I has been heoreically shown ha acive paricles moving a consan speed in a heerogeneous wo-dimensional space eperience self rapping leading o subdiffusion [3]. Usually i is assumed ha subdiffusive behavior is caused by he paricle being rapped in local minima for prolonged imes before i escapes o a neighboring minima. Coninuous ime random walks CTRWs wih on-sie waiing-ime disribuions falling slowly as α 1 predic a subdiffusive behavior [4, 5]. Superdiffusion processes, characerized by he nonlinear dependence 1 of he MSD on ime wih he power-law eponen in he range 1 < α <, consiue anoher subclass of anomalous diffusion processes. Superdiffusion is observed, for eample, in vibraed dense granular media [6]. Theoreical models sugges ha supperdiffusion can be caused by Lévy flighs [4]. Lévy flighs resuling in a superdiffusion can be modeled by fracional Fokker-Planck equaions [7] or Langevin equaions wih an addiive Lévy sable noise. In many eperimenal sudies i is only possible o show ha signal inensiy disribuion has Lévy law-ails: disribuion funcion of urbulen magneized plasma emiers [8] and sep-size disribuion of phoons in ho vapors of aoms [9] have Lévy ails. This indirecly shows ha in hese sysems superdiffusion could be found. Anomalous diffusion does no uniquely indicae he processes occurring in he sysem, because here are differen sochasic processes sharing he behavior of he MSD 1. The physical mechanisms leading o he deviaions from he linear ime dependence of he MSD can depend on he sysem or on he emporal and spaial ranges under consideraion. For eample, diffusion described by CTRW has been observed for sub-micron racers in biological cells [1 1], srucured coloidal sysems [13] and for charge carrier moion in amorphous semiconducors [5, 14]. Fracional Brownian moion and fracional Langevin equaions has been used o model he dynamics in membranes [15, 16], moion of polymers in cells [17], racer moion in comple liquids [18, 19]. Diffusion of even smaller racers in biological cells has been described by a spaially varying diffusion coefficien []. ryis.kazakevicius@fai.vu.l

2 Recenly, in Refs. [1 4] i was suggesed ha he anomalous diffusion can be a resul of heerogeneous diffusion process HDP, where he diffusion coefficien depends on he posiion. Spaially dependen diffusion can occur in heerogeneous sysems. For eample, heerogeneous medium wih seep gradiens of he diffusiviy can be creaed in hermophoresis eperimens using a local variaion of he emperaure [5, 6]. Mesoscopic descripion of ranspor in heerogeneous porous media in erms of space dependen diffusion coefficiens is used in hydrology [7, 8]. In urbulen media he Richardson diffusion has been described by heerogeneous diffusion processes [9]. Power-law dependence of he diffusion coefficien on he posiion has been proposed o model diffusion of a paricle on random fracals [3, 31]. In bacerial and eukaryoic cells he local cyoplasmic diffusiviy has been demonsraed o be heerogeneous [, 3]. Moion of a Brownian paricle in an environmen wih a posiion dependen emperaure has been invesigaed in Ref. [33]. In random walk descripion he spaially varying diffusion coefficien can be included via posiion dependence of he waiing ime for a jump even [34], he posiion dependence occurs because in he heerogeneous medium he properies of a rap can reflec he medium srucure. This is he case for diffusion on fracals and mulifracals [35]. Inhomogeneous versions of coninuous ime random walk models for waer permeaion in porous ground layers were proposed in Ref. [36]. Heerogeneous diffusion process migh be applicable o describe anomalous diffusion in such sysems. The goal of his paper is o consider HDPs wih he power-law dependence of he diffusion coefficien on he posiion and o analyically invesigae he influence of eernal forces on he resuling anomalous diffusion. The influence of he eernal forces on HDPs has no been sysemaically analyzed. In his paper we assume ha he forces are characerized by a power-law dependence on he posiion. Such forces can arise in various sysems: In many cases he poenials causing he deerminisic forces are power-law funcions of posiion, for eample linear and harmonic poenials. Power-law poenial wih arbirary power-law eponen acing on a nanoparicle have been creaed in Ref. [37]. Logarihmic poenials yielding forces behaving as 1 have been applied o describe dynamics of paricles near a long, charged polymer [38], momenum diffusion in dissipaive opical laices [39], long-range ineracing sysems [4], bubbles in DNA molecules [41]. As we demonsrae, he eernal force having a cerain value of he power-law eponen does no resric he region of diffusion and change only he anomalous diffusion coefficien wihou changing he scaling eponen α. Oher values of he power-law eponen in he eernal force can lead o he eponenial cu-off of he probabiliy densiy funcion PDF and resric he region of diffusion, limiing he ime inerval when he anomalous diffusion occurs. We epec ha he resuls obained in his paper may be relevan for a more complee undersanding of anomalous diffusion processes. The paper is organized as follows: In Sec. II we summarize he main properies of HDPs wih he power-law dependence of he diffusion coefficien on he posiion. The influence of an eernal force wih a paricular value of he power-law eponen is invesigaed in Sec. III. This value of he power-law eponen is he same as in he drif correcion for ransformaion from he Sraonovich o he Ió sochasic equaion. We consider eernal forces having oher values of he power-law eponen in Sec. IV. Sec. V summarizes our findings. II. FREE HETEROGENEOUS DIFFUSION PROCESS Heerogeneous diffusion process wih he power-law dependence of he diffusion coefficien on he posiion, inroduced in Ref. [1], is described by he Langevin equaion d = σ η dw. Here η is he power-law eponen of muliplicaive noise, σ is he ampliude of noise and W is a sandard Wiener process Brownian moion. This sochasic differenial equaion SDE is inerpreed in Sraonovich sense. For mahemaical convenience and for furher generalizaion we ransform Eq. o he Iô convenion: The Fokker-Planck equaion corresponding o he SDE is [4] d = 1 σ η η 1 d+σ η dw. 3 P, = 1 σ [ η η P, ]. 4 Wih he reflecive boundaries a small posiive = min and large = ma, Eq. 4 leads o he seady-sae PDF P η. Wihou suchboundariesheime-dependen soluionofeq.4wihhe iniialcondiionp, = δ is given by a sreched when η > or compressed when η < Gaussian [1] P, = η πσ ep 1 η 1 η σ. 5

3 3 When η < 1, his soluion describes eponenial cu-off a large values of and power-law behavior a small values of. The posiion of he cu-off moves owards large values of wih increase of ime. In conras, when η > 1 he soluion 5 describes eponenial cu-off a small values of and power-law behavior a large values of. The posiion of cu-off moves owards smaller values of wih increase of ime. In Ref. [1] i has been demonsraed ha Eq. or, equivalenly, Eq. 3 leads o he power-law ime dependence of he MSD σ 1 1 η. 6 This behavior of he MSD means ha HDP described by Eq. yields superdiffusion for 1 > η > and subdiffusion for η <. For η > he diffusiviy increases wih increasing, leading o a progressive acceleraion of he diffusing paricle. In conras, he diffusiviy decreases wih increasing when η <. For η > 1 he paricle becomes localized. When η = 1, he MSD grows no as a power-law of ime, bu eponenially [43, 44]. The HDP displays weak non-ergodiciy, ha is he scaling of ime and ensemble averages is differen. Specifically, in Ref. [1] i has been shown ha he average over he rajecories of he he ime-averaged MSD δ = 1 T δ = 1 N T N i=1 δi 7 [+ ] d 8 scales as δ T η η 1. 9 Thus ime-averaged MSD depends on he ime difference linearly, in conras o he power-law behavior of MSD in Eq. 6. Anoher ineresing propery of HDPs is he behavior of he disribuion of he ime-averaged MSD δ of individual realizaions. When η <, he disribuion of δ decays o zero a δ = [4]. This behavior of he disribuion in HDPs is differen han he behavior in CTRWs, where here is a finie fracion of immobile paricles resuling in he finie value of he disribuion a δ =. This difference allows o disinguish beween differen origins of anomalous diffusion. III. EXTERNAL FORCE THAT DOES NOT LIMIT THE ANOMALOUS DIFFUSION In general, no only a random force leading o he diffusion bu also a deerminisic drif force can be presen in he Langevin equaion. Therefore, he quesion arises how eernal force influences he anomalous diffusion described by heerogeneous diffusion process. To invesigae his quesion we will consider Eq. 3 wih an addiional drif erm. As he resuls obained in his Secion show, for he appearance of he anomalous diffusion i is sufficien o consider only posiive values of. Therefore from now on we will wrie jus insead of he absolue value, assuming he presence of he boundary ha does no allow for he diffusing paricle o ener he region. Due o he power-law dependence of he diffusion coefficien on he posiion, he drif erm also being a power-law funcion of he posiion allows us o obain analyical epressions. The eernal forces characerized by power-law dependence on he posiion can be creaed eperimenally, for eample, a power-law poenial wih arbirary powerlaw eponen acing on a nanoparicle have been creaed in Ref. [37]. Firs of all, in his Secion we will consider he case where he eernal force has he same dependence on coordinae as he drif correcion for he ransformaion from he Sraonovich o he Ió sochasic equaion. Such a drif erm can arise no only due o an eernal force bu can also represen a noise-induced drif [45]. Thus we will generalize he SDE 3 o ake he form d = σ η ν η 1 d+σ η dw. 1 Hereν isanew parameerdescribinghe addiionaldrif erm. The meaning ofhe parameerν isasfollows: when he reflecive boundaries a small posiive = min and large = ma are presen, he seady-sae PDF is a power-law

4 funcion of posiion wih he power-law eponen ν, P ν. Comparison of Eq. 1 wih Eq. 3 shows ha he free HDP is obained when ν = η. TheSDE1hasbeen proposedinrefs.[46,47]ogeneraesignalshaving1/f β noiseinawiderangeoffrequencies. Such nonlinear SDEs have been applied o describe signals in socio-economical sysems [48, 49] and as a model of neuronal firing [5]. According o Ref. [47], he power-law eponen β in he power specral densiy Sf f β is relaed o he parameers of SDE 1 as β = 1+ ν 3 η For some values of he parameer ν he SDE 1 can be obained from Eq. by changing he prescripion for calculaing sochasic inegrals. Le us consider he SDE ogeher wih he inerpreaion of sochasic inegrals as [45, 51] T d = σ η γ dw, 1 N 1 f γ dw = lim f n W n, N n= 4 n = n+γ N T. 13 The parameer γ, γ 1, defines he prescripion for calculaing sochasic inegrals. Commonly used values of he parameer γ are γ = corresponding o pre-poin Iô convenion, γ = 1/ corresponding o mid-poin Sraonovich convenion and γ = 1 corresponding o pos-poin Hänggi-Klimonovich [5, 53], kineic or isohermal convenion [54 56]. The inegraion convenion should be deermined from he eperimenal daa or derived from anoher model [57]. The SDE 1 also has been invesigaed in Ref. [58]. Transformaion of Eq. 1 o he Iô equaion yields [51, 59, 6] This equaion has he form of SDE 1 wih he parameer ν being d = σ γη 1 η 1 d+σ η dw. 14 ν = [γ +1 γη]. 15 Since γ 1, he range of possible values of he parameer ν obained by changing he prescripion in Eq. is limied: η ν when η < 1 and ν η when η > 1. In his Secion we do no place such resricion on he possible values of ν. The values of ν ouside of his range can be obained due o he acion of an eernal force. The Fokker-Planck equaion corresponding o he SDE 1 is [4] P, = σ ν η [η 1 P,]+ σ [η P,]. 16 The ime-dependen PDF of he process given by Eq. 1 can be obained as follows: Transformaion of he variable o a new variable y = 1 η assuming ha η 1 leads o he SDE where dy = 1 σ ν 1 y d+σ dw, 17 ν = η ν η 1, σ = η 1 σ. 18 Equaion 17 has he form of a Bessel process [61]. This connecion wih he Bessel process has also been poined ou in Ref. [58]. The known analyic form of he soluion of he Fokker-Planck equaion corresponding o SDE 17 is [61 63] P y = 1 σ ν y y 1 P y + 1 σ y P y 19 Py, y, = y 1 ν y 1+ν σ ep y +y yy σ I ν +1 σ.

5 Here I n z is he modified Bessel funcion of he firs kind. This PDF obeys he iniial condiion Py, = y, = δy y. For compleeness, one possible way of obaining his soluion of Eq. 19 is presened in Appendi A. The PDF can be normalized and represens an accepable soluion only when ν +1 > 1, ha is when ν < 1. When his inequaliy is no saisfied, he Bessel process leads o a oal absorpion a he origin in a finie ime [64]. Transforming he variables back we obain he ime-dependen PDF saisfying he iniial condiion P, =, = δ : P,, = 1 η ν 1 η+ν η 1 σ ep 1 η + 1 η η 1 σ Iν+1 η η η 1 η η 1 σ. 1 Similarly as Eq., he soluion 1 is valid when ν+1 η η 1 > 1. This condiion is equivalen o ν > 1 and η > 1 or ν < 1 and η < 1. Similar epression of he PDF 1 has been obained in Ref. [58] by considering he HDP wih various values of he prescripion parameer γ. Using Eq. 1 we can calculae he ime-dependen average of a power of : a = = a P,,d a η 1 η 1 η 1 σ a 1 η 1 F 1 a η 1 ; ν 1 η 1 1 η ; η 1 σ. Here 1 F 1 a;b;z is he Kummer confluen hypergeomeric funcion. The inegral is finie when: i η > 1 and ν > 1+a wih a > or ν > 1 wih a < ; ii η < 1 and ν < 1+a wih a < or ν < 1 wih a >. Eq. for he a-h momen of has been derived also in Ref. [58]. In paricular, he average of is equal o = ν η 1 η 1 η 1 σ 1 1 η 1 F 1 1 η 1 ; ν 1 η 1 1 η ; η 1 σ 3 and is finie when ν > and η > 1 or ν < 1 and η < 1. The average of he square of is equal o ν 3 η 1 = η 1 σ η 1 F 1 η 1 ; ν 1 1 η ; η 1 η 1 σ η 1 4 and is finie when ν > 3 and η > 1 or ν < 1 and η < 1. When ime is large, ha is when 1 η η 1 σ 1, 5 he hypergeomeric funcion in Eq. is approimaely equal o 1. Thus for large ime he average a does no depend on he iniial posiion and is a power-law funcion of ime: a a η 1 η 1 σ a 1 η. 6 η 1 As a consequence we obain ha for large ime saisfying he condiion 5 he average of he square of depends on ime as 1/1 η. In addiion, using Eqs. 3 and 4 we ge ha he he variance = has he same dependence on ime: 1/1 η. We can conclude, ha no only he original HDP equaion, bu also he equaion 1 wih addiional power-law force ehibi anomalous diffusion. The power-law eponen in he ime dependence is he same as in Eq. 6, he eernal force considered in his Secion does no change he characer of he anomalous diffusion as long as he condiion 5 is saisfied. The power-law eponen 1/1 η depends only on he random force in he SDE 1 and does no depend on he parameer ν characerizing he deerminisic eernal force. The eernal force changes only he anomalous diffusion coefficien. We check he analyic resuls obained in his Secion by comparing hem wih numerical simulaions. The powerlaw form of he coefficiens in SDE 1 allows us o inroduce an operaional ime τ in addiion o he physical ime

6 6 a 1 1 b c d e f FIG. 1. Dependence of he mean a,b,c and variance d,e,f on ime for various values of he parameers η and ν when he posiion of he diffusing paricle changes according o Eq. 1. Solid gray lines show numerical resul, dashed black lines are calculaed using Eqs. 3 and 4, doed red lines show he power-law dependence on ime 1/[1 η] for a,b,c and 1/1 η for d,e,f. The parameers are σ = 1 and η = 1, ν = 1 for a,d; η = 1, ν = for b,c; η = 3, ν = 5 for c,f. The iniial posiion is = 1. so ha he diffusion coefficien in he operaional ime becomes consan [65]. The relaion beween he physical ime and he operaional ime τ is specified by he equaion d = σ η dτ. For he numerical soluion of SDE 1 we use he Euler-Maruyama scheme wih a variable ime sep k = τ/σ η k which is equivalen o he inroducion of he operaional ime [65]. Thus he numerical mehod of soluion of SDE 1 is given by he equaions k+1 = k + η ν τ + τε k, k 7 k+1 = k + τ σ η. k 8 Here τ 1 is he ime sep in he operaional ime and ε k are normally disribued uncorrelaed random variables wih a zero epecaion and uni variance. To avoid he divergence of he diffusion and drif coefficiens a = in he numerical simulaion, we inser a reflecive boundary a = 1 3. This modificaion is analogous o he regularizaion of he diffusiviy a =, performed in Refs. [1, 4]. When η > 1 and ν > 1 he PDF of increases for small values of. In his case he operaional ime inroduced in such a way ha he coefficien before noise becomes proporional o he firs power of can lead o a more efficien numerical mehod. Such an operaional ime is inroduced by he equaion d = σ η 1 dτ and he numerical mehod of soluion becomes k+1 = k [1+ τ η ν + ] τε k, 9 τ k+1 = k + σ η 1 k. 3 We have calculaed he mean and he variance by averaging over 1 5 rajecories. Comparison of he analyic epressions 3, 4 wih numerically obained ime-dependen mean and variance is shown in Fig. 1. For numerical simulaion we have chosen hree differen values of he eponen η: η = 1, η = 1, and η = 3 corresponding, respecively, o subdiffusion, superdiffusion and o he localizaion of he paricle. For each case we have chosen a value of he parameer ν ha differs from η of he free HDP. For all numerical simulaions he iniial posiion is = 1. We see a good agreemen of he numerical resuls wih analyic epressions. As we can see in Fig. 1, he ime dependence of he mean and he variance becomes a power-law for large imes and he iniial posiion is forgoen for

7 7 a P b P c P FIG.. Time-dependen PDF P,, corresponding o imes = 1 ligh gray and = 1 dark gray for various values of he parameers η and ν when he posiion of he diffusing paricle changes according o Eq. 1. Dashed black lines are calculaed using Eq. 1. The doed line shows he slope ν. The parameers are σ = 1 and a η = 1, ν = 1; b η = 1, ν = ; c η = 3, ν = 5. The iniial posiion is = 1. he parameers used in Fig. 1b and e he difference beween he eac soluion and he power-law approimaion remains consan, bu he relaive difference is decreasing. Comparison of he analyic epression 1 for he ime-dependen PDF P,, wih he resuls of numerical simulaion is shown in Fig.. To illusrae how he PDF changes wih ime, he PDF is shown a wo differen ime momens = 1 and = 1. We see a good agreemen of he numerical resuls wih he analyic epression for boh ime momens. Wih increasing ime he PDF shifs o he larger values of when η < 1 and o he smaller values of when η > 1. IV. EXTERNAL FORCE LEADING TO AN EXPONENTIAL RESTRICTION OF THE DIFFUSION Now le us consider he eernal deerminisic force having a power-law dependence on bu wih he power-law eponen differen han η 1. When such a force is posiive if he power-law eponen is smaller han η 1 and negaive if he power-law eponen is larger han η 1, he SDE describing he HDP can be wrien as d = σ η ν η 1 + m 1 m1 min m d+σ η 1 m1 η 1+m η dw m. 31 Here m 1,m > and min, ma are he parameers of he eernal force. In Eq. 31 we included hree erms in he drif. Each erm has power-law dependence on posiion, bu wih differen power-law eponens: equal o η 1, smaller han η 1 and larger han η 1. The seady-sae PDF corresponding o Eq. 31 is { P ν min m } m1 ep. 3 We see ha he addiional erms in Eq. 31 lead o eponenial resricion of he diffusion. When min ma, he seady-sae PDF has he power-law form P ν. Confined HDP has been invesigaed in Refs. [66, 67]. Analysis of confined HDP can be relevan for he descripion of he racer paricles moving in he confinemen of cellular comparmens or for he paricle raced wih opical weezers ha eer a resoring force on he paricle [1]. We can mahemaically obain one of he erms in Eq. 31 leading o he eponenial resricion of he diffusion by ransforming he ime in he iniial equaion. Indeed, if we sar wih he Iô SDE ma ma d = ad+bdw 33 and inroduce a new sochasic variable z = gc hen he SDE for he sochasic variable z becomes dg z z dc z dz = d g +gdc a d+g b dw. 34 d g d g In Eq. 34 a new erm ha is proporional o z and o he derivaive of g appears in he drif. The PDF of he sochasic variable z is relaed o he PDF of via he equaion P z z, = 1 z g P g,c. 35

8 Thus, o inroduce a new erm ino Eq. 1, le us sar wih a sochasic variable y obeying he SDE 1 and consider a new sochasic variable 1 = e µ y κ eκ Here he funcions g and c are g = e µ and c = κ 1 e κ 1. From Eq. 34 follows ha he equaion for he sochasic variable is d = µ+σ η ν e κ µη 1 η 1 d+σe 1 κ µη 1 η dw. 37 We can obain an equaion wih ime-independen coefficiens by requiring ha κ = µη Using his value for he parameer κ he SDE for becomes d = µ+σ η ν η 1 d+σ η dw. 39 When µ has he same sign as η 1, his SDE can be wrien in he form similar o Eq. 31: d = σ η ν +η 1 m η 1 η 1 d+σ η dw, 4 where he parameer m is defined by he equaion µ = σ η 1 η 1 m. 41 Comparing Eq. 4 and Eq. 31 we see ha he ime ransformaion considered in his Secion inroduces an eponenial resricion of he diffusion a small values of when η > 1 and a large values of when η < 1. Using Eqs. 1 and 35 we obain he ime-dependen PDF for he sochasic variable obeying SDE 4: 8 η 1 η 1 m 1 ν η 1+ν η µη 1 e 1+ν η µ P,, = 1 e η 1 m ep 1 η + 1 η 1 e η 1µ 1 η 1 η I 1+ν η η 1 η 1 m sinhη 1µ e η 1µ 4 The condiions of validiy of his epression is he same as for Eq. 1. Tha is, he epression for he PDF given by Eq. 4 is valid when ν > 1 and η > 1 or ν < 1 and η < 1. The average of a power of, calculaed using Eq. 4, is a = = a P,,dy a η 1 η 1 a m 1 e η 1µ a 1 F 1 η 1 a η 1 ; ν 1 ; η 1 m η 1 1 η e η 1µ This average is finie under he same condiions as Eq.. In paricular, he average of is equal o ν η 1 = m 1 1 e η 1µ 1 1F 1 η 1 η 1 ; ν 1 ; η 1 m 1 η η 1 e η 1µ 1 η 1 and is finie when ν > and η > 1 or ν < 1 and η < 1. The average of he square of is equal o ν 3 η 1 = m 1 1F 1 e η 1µ 1 1 η 1 η 1 ; ν 1 ; η 1 m 1 η η 1 e η 1µ 1 η

9 and is finie when ν > 3 and η > 1 or ν < 1 and η < 1. When µ has he same sign as η 1 and hen he PDF 4 ends o he seady-sae PDF η 1 m P = ν ep η 1 m The seady-sae PDF 46 leads o he seady-sae averages of and s = s = ν η 1 η 1 ν 3 η 1 η 1 η m, 47 m. 48 Now le us consider he ime evoluion of he average, given by Eq. 45. In he case when he iniial posiion is far from he cu-off boundary m ha is m when η < 1 and m when η > 1, he ime evoluion of he average can be separaed ino hree pars. Firs, for small imes 9 1 η η 1 σ he influence of he iniial posiion is significan and he diffusion is approimaely normal, depends linearly on ime. For he inermediae imes 1 η η 1 σ 1 η 1µ = 1 η m η 1 σ he eponen e η 1µ in Eq. 45 differs from 1 only slighly, however he las argumen of he hypergeomeric funcion is already small. Approimaing he hypergeomeric funcion by 1 and epanding he eponen e η 1µ ino power series and keeping only he linear erm we obain ha he average depends on ime as a power-law, 1/1 η. Thus for his inermediae range of ime he anomalous diffusion occurs. Finally, for large imes 1 η 1µ he cu-off posiion m sars o influence he diffusion and approaches he seady-sae value 48. We can conclude ha he inroducion of he boundary via an eponenial cu-off does no change he anomalous diffusion when he saring posiion is far from he boundary an he ime is no oo large. Comparison of he analyic epressions 44, 45 wih numerically obained ime-dependen mean and variance is shown in Fig. 3. As in Sec. III, for numerical soluion we use he Euler-Maruyama scheme wih a variable ime sep, equivalen o he inroducion of he operaional ime in addiion o he physical ime. When he diffusion coefficien in he operaional ime τ does no depend on posiion, he numerical mehod is given by he equaions k+1 = k + [ η ν 1 +η 1 η 1 m 1 η k k ] τ + τε k, 49 k+1 = k + τ σ η. 5 k Whenη > 1,amoreefficiennumericalmehodisobainedwhenhechangeofhevariableinonesepisproporional o he value of he variable. Then he numerical mehod of soluion is described by he equaions k+1 = k [1+ τ η ν η 1 m + τη 1 + ] τε k, 51 k+1 = k + τ σ η 1 k k. 5

10 1 a 1 1 b 1 c d e f FIG. 3. Dependence of he mean a,b,c and variance d,e,f on ime for various values of he parameers η and ν when he posiion of he diffusing paricle changes according o Eq. 4. Solid gray lines show numerical resul, dashed black lines are calculaed using Eqs. 44 and 45, doed red lines show he power-law dependence on ime 1/[1 η] for a,b,c and 1/1 η for d,e,f. The parameers are σ = 1 and η = 1, ν = 1, m = 5 for a,d; η = 1, ν =, m = 1 for b,c; η = 3, ν = 5 m =.1 for c,f. The iniial posiion is = 1. a b c 1 P P P FIG. 4. Time-dependen PDF P,, corresponding o imes = 1 ligh gray and = 1 dark gray for various values of he parameers η and ν when he posiion of he diffusing paricle changes according o Eq. 4. Dashed black lines are calculaed using Eq. 4. The doed line shows he seady-sae PDF 46. The parameers are σ = 1 and a η = 1, ν = 1, m = 5; b η = 1, ν =, m = 1; c η = 3, ν = 5 m =.1. The iniial posiion is = 1. When η < 1, o avoid he divergence of he diffusion and drif coefficiens a = in he numerical simulaion, we inser a reflecive boundary a = 1 3. This is no necessary when η > 1 because he addiional erm in he drif creaes an eponenial cu-off a small values of m. For numerical simulaion we used he same values of he parameers η and ν as in Fig. 1, he iniial posiion is = 1. We have chosen he parameer m of he eernal force so ha he iniial posiion is far from he boundary m. In Fig. 3 we see a good agreemen of he numerical resuls wih analyic epressions. The numerical calculaion and he analyic epressions confirm he presence of a ime inerval where he mean and he variance have a power-law dependence on ime, as can be seen in Fig. 3. The upper limi of his ime inerval is deermined by he he posiion m of he eponenial cu-off. Comparison of he analyic epression 4 for he ime-dependen PDF P,, wih he resuls of numerical simulaion is shown in Fig. 4. We see a good agreemen of he numerical resuls wih he analyic epression. Wih increasing ime he PDF shifs o he larger values of when η < 1 and o he smaller values of when η > 1. However, in conras o he siuaion in he previous Secion where he resricing force was no presen, his shif of he PDF now is limied, a large imes he ime-depend PDF approaches he seady-sae PDF 46.

11 11 V. CONCLUSIONS In summary, we have obained analyic epressions 1 and 4 for he ransiion probabiliy of he heerogeneous diffusion process whose diffusiviy has a power-law dependence on he disance. In he descripion of he HDP we have included an addiional deerminisic force ha also has a power-law dependence on he posiion. The drif erm having a power-law dependence on he posiion can arise no only due o an eernal force bu can also represen a noise-induced drif [45]. Such a drif erm appears in a Langevin equaion describing overdamped flucuaions of he posiion of a paricle in nonhomogeneous medium [68]. Sochasic differenial equaions wih power-law drif and diffusion erms have been used o model random flucuaions of he amospheric forcing on he ocean circulaion [69] and pressure ime series rouinely used o define he inde characerizing he Norh Alanic Oscillaion [7]. The Brownian moion of a colloidal paricle in waer subjeced o he graviaional force and wih a space-dependen diffusiviy due o he presence of he boom wall of he sample cell has been invesigaed in Ref. [71]. Force causing eponenial cu-off in he PDF of he paricle posiion can describe HDP process in confined regions; such a descripion can be relevan, for eample, for he racer paricles moving in he confinemen of cellular comparmens [1, ]. A sysem obeying SDEs wih power-law drif and diffusion erms can be eperimeally realized as an elecrical circui driven by a muliplicaive noise, similarly as in Ref. [7]. The equaion describing he overdamped moion of he Brownian paricle akes he form of Eq. 1 when a emperaure gradien is presen in he medium and he paricle is subjeced o he eernal poenial ha is proporional o he emperaure profile [33]. In paricular, seady sae hea ransfer due o he emperaure difference beween he beginning and he end of he sysem corresponds o η = 1/ in Eq. 1 [33]. For a charged paricle he eernal force can be inroduced by applying he elecric poenial difference a he ends of he sysem. When he power-law eponen in he deerminisic force is equal o η 1, where η is he power-law eponen in he dependence of he diffusion coefficien on he posiion, he eernal force does no limi he region of diffusion. Oher values of he power-law eponen in he deerminisic force can cause he eponenial cu-off in he PDF of he paricle posiions. Such an eponenial resricion of he diffusion appears when he eernal force is posiive if he power-law eponen is smaller han η 1 and negaive if he power-law eponen is larger han η 1. We obained an analyic epression 4 for he ransiion probabiliy in a paricular case when he eernal resricing force has a linear dependence on he posiion. Using analyic epression for he ransiion probabiliy we calculaed he ime dependence of he momens of he paricle posiion, Eqs. and 43. We found ha he power-law eponen in he dependence of he MSD on ime does no depend on he eernal force, his force changes only he anomalous diffusion coefficien. This conclusion is valid for sufficienly large imes saisfying he condiion 5, ha is when he iniial posiion of he paricle is forgoen and he anomalous diffusion occurs. In addiion, he eernal force having he power-law eponen differen from η 1 limis he ime inerval where he anomalous diffusion occurs. The conclusions remain valid also when he eernal force can be represened as a sum of several erms, each erm being a power-law funcion of posiion wih a differen power-law eponen. Thus, our resuls indicae ha he anomalous diffusion caused by diffusiviy being a power-law funcion of he posiion is robus wih respec o an eernal perurbaion, he eponen α in Eq. 1 is deermined only by he diffusion coefficien. In addiion, he resuls of Sec. III show ha he characer of he anomalous diffusion does no depend on he inerpreaion of he Langevin equaion: he scaling eponen α in Eq. 1 is he same for boh Sraonovich and Iô convenions. This is because differen inerpreaions correspond o he differen values of he parameer ν in Eq. 1: Sraonovich convenion resuls in ν = η, Iô convenion in ν = η, and he scaling eponen α does no depend on ν. The same conclusion ha he eponen of he anomalous diffusion does no depend on he prescripion has been obained in Refs. [58, 73, 74] for equaions describing diffusion wihou he presence of an eernal force. Appendi A: Soluion of he Fokker-Planck equaion for he Bessel process Le us consider he Fokker-Planck equaion P = ν 1 P + 1 P and search for he ime-dependen soluion P,, wih he iniial condiion P,, = δ. The unnormalized ime-independen soluion of Eq. A1 is ν. The boundary condiion a = for Eq. A1 can be epressed using he probabiliy curren [75] A1 S, = ν P, 1 P,. A

12 We consider he boundary condiion corresponding o he vanishing probabiliy curren a =, S, =. One of he possible ways o obain he soluion of Eq. A1 is o use a Laplace ransform [61]. Here we solve Eq. A1 using he mehod of eigenfuncions. This mehod has been used in Ref. [76] for an equaion, similar o Eq. A1. An ansaz of he form 1 P, = P λ e λ A3 leads o he equaion ν 1 P λ + 1 P λ = λp λ, where P λ are he eigenfuncions and λ are he corresponding eigenvalues. The eigenfuncions obey he orhonormaliy relaion [75] A4 ν P λ P λ d = δλ λ. A5 Epansion of he ransiion probabiliy densiy P,, in erms of he eigenfuncions has he form [75] P,, = P λ ν P λ e λ dλ. A6 For he soluion of Eq. A4 i is convenien o wrie he eigenfuncions P λ as P λ = 1 γ u λ, A7 where γ = 1+ν Similar anzaz has been used in Refs. [6, 63]. The funcions u λ obey he equaion. A8 d d u λ + d d u λ +ρ γ u λ =, A9 where ρ = λ. A1 The probabiliy curren S λ, Eq. A, rewrien in erms of funcions u λ, becomes S λ = 1 γ γu λ + d d u λ. A11 The orhonormaliy of eigenfuncions A5 yields he orhonormaliy for funcions u λ u λ u λ d = δλ λ. A1 The general soluion of Eq. A9 is u λ = c 1 J γ ρ+c Y γ ρ, A13 where J γ and Y γ are he Bessel funcions of he firs and second kind, respecively. The coefficiens c 1 and c needs o be deermined from he boundary and normalizaion condiions for he funcions u λ. Using Eqs. A13 and A11 we ge he probabiliy curren S λ = 1 ρ1 γ [c 1 J γ 1 ρ+c Y γ 1 ρ]. A14

13 The requiremen S λ = leads o he condiion c = c 1 anπγ. Taking ino accoun his relaion beween c and c 1 we obain 13 u λ = c λ J γ ρ. A15 In addiion, he condiion S λ = implies c = when γ 1. Thus, when γ 1, ν 1 boh coefficiens c 1 and c are zero and he soluion A6 is no valid. I is known ha for a Bessel process wih such parameers a oal absorpion a he origin occurs in a finie ime [64]. Orhonormaliy condiion A1 leads o he equaion c λ c λ Since for he Bessel funcions he equaliy J γ ρj γ ρ d = δλ λ. A16 ρ J γ ρj γ ρ d = δρ ρ is valid, we obain c λ = 1. Using Eqs. A6, A7 and A15 he soluion of he Fokker-Planck equaion can be epressed as Inegraion yields P,, = 1 γ γ P,, = 1 γ γ Here I α is he modified Bessel funcion of he firs kind. J γ ρj γ ρ e 1 ρ ρdρ. A17 A18 ep + I γ. A19 [1] J.-P. Bouchaud and A. Georges, Phys. Rep. 195, [] F. Evers, C. Zunke, R. D. L. Hanes, J. Bewerunge, I. Ladadwa, A. Heuer, and S. U. Egelhaaf, Phys. Rev. E 88, [3] O. Chepizhko and F. Peruani, Phys. Rev. Le. 111, [4] R. Mezler and J. Klafer, Phys. Rep. 339, 1. [5] M. Schuber, E. Preis, J. C. Blakesley, P. Pingel, U. Scherf, and D. Neher, Phys. Rev. B 87, [6] C. Scallie, A. Gnoli, A. Puglisi, and A. Vulpiani, Phys. Rev. Le 114, [7] H. C. Fogedby, Phys. Rev. Le. 73, [8] Y. Marande, H. Capes, L. Godber-Moure, R. Guirle, M. Koubii, and R. Samm, Commun. Nonlinear. Sci. Commun. 8, [9] N. Mercadier, W. Guerin, M. Chevrollier, and R. Kaiser, Na. Phys. 5, 6 9. [1] S. M. A. Tabei, S. Burov, H. Y. Kim, A. Kuznesov, T. Huynh, J. Jureller, L. H. Philipson, A. R. Dinner, and N. F. Scherer, Proc. Nal Acad. Sci. USA 11, [11] A. V. Weigel, B. Simon, M. M. Tamkun, and D. Krapf, Proc. Nal Acad. Sci. USA 18, [1] J.-H. Jeon, V. Tejedor, S. Burov, E. Barkai, C. Selhuber-Unkel, K. Berg-Sørensen, L. Oddershede, and R. Mezler, Phys. Rev. Le. 16, [13] I. Y. Wong, M. L. Gardel, D. R. Reichman, E. R. Weeks, M. T. Valenine, A. R. Bausch, and D. A. Weiz, Phys. Rev. Le. 9, [14] H. Scher and E. W. Monroll, Phys. Rev. B 1, [15] J.-H. Jeon, H. Marinez-Seara Monne, M. Javanainen, and R. Mezler, Phys. Rev. Le. 19, [16] G. R. Kneller, K. Baczynski, and M. Pasienkewicz-Gierula, J. Chem. Phys 135, [17] E. Kepen, I. Bronshein, and Y. Garini, Phys. Rev. E 87, [18] J. Szymanski and M. Weiss, Phys. Rev. Le. 13, [19] J.-H. Jeon, N. Leijnse, L. B. Oddershede, and R. Mezler, New J. Phys. 15, [] T. Kühn, T. O. Ihalainen, J. Hyväluoma, N. Dross, S. F. Willman, J. Langowski, M. Vihinen-Rana, and J. Timonen, PLoS ONE 6, e [1] A. G. Chersvy, A. V. Chechkin, and R. Mezler, New J. Phys. 15, [] A. G. Chersvy and R. Mezler, Phys. Chem. Chem. Phys 15, 13. [3] A. G. Chersvy, A. V. Chechkin, and R. Mezler, Sof Maer 1,

14 [4] A. G. Chersvy and R. Mezler, Phys. Rev. E 9, [5] Y. T. Maeda, T. Tlusy, and A. Libchaber, Proc. Nal. Acad. Sci. 19, [6] C. B. Mas, S. Schink, U. Gerland, and D. Braun, Proc. Nal. Acad. Sci. 11, [7] R. Haggery and S. M. Gorelick, Waer Resources Res. 31, [8] M. Denz and D. Bolser, Phys. Rev. Le. 15, [9] L. F. Richardson, Proc. R. Soc. Lond. A 11, [3] B. O Shaughnessy and I. Procaccia, Phys. Rev. Le. 54, [31] C. Loverdo, O. Bénichou, R. Voiuriez, A. Biebricher, I. Bonne, and P. Desbiolles, Phys. Rev. Le. 1, [3] B. P. English, V. Hauryliuk, A. Sanamrad, S. Tankov, N. H. Dekker, and J. Elf, Proc. Nal Acad. Sci. USA 18, E [33] R. Kazakevičius and J. Ruseckas, J. Sa. Mech. 15, P1 15. [34] T. Srokowski, Phys. Rev. E 89, 31R 14. [35] D. Scherzer, M. Larchevêque, J. Duan, V. V. Yanovsky, and S. Lovejoy, J. Mah. Phys. 4, 1. [36] M. Denz, P. Gouze, A. Russian, J. Dweik, and F. Delay, Adv. Waer Res. 49, [37] A. E. Cohen, Phys. Rev. Le. 94, [38] G. S. Manning, J. Chem. Phys. 51, [39] Y. Sagi, M. Brook, I. Almog, and N. Davidson, Phys. Rev. Le. 18, [4] F. Bouche and T. Dauois, Phys. Rev. E 7, 4513R 5. [41] L. A. Wu, S. S. Wu, and D. Segal, Phys. Rev. E 79, [4] C. W. Gardiner, Handbook of Sochasic Mehods for Physics, Chemisry and he Naural Sciences Springer-Verlag, Berlin, 4. [43] A. W. C. Lau and T. C. Lubensky, Phys. Rev. E 76, [44] A. Fuliński, Phys. Rev. E 83, [45] G. Volpe and J. Wehr, Rep. Prog. Phys. 79, [46] B. Kaulakys and J. Ruseckas, Phys. Rev. E 7, 11R 4. [47] B. Kaulakys, J. Ruseckas, V. Gonis, and M. Alaburda, Physica A 365, [48] V. Gonis, J. Ruseckas, and A. Kononovicius, Physica A 389, 1 1. [49] J. Mahiesen, L. Anghelua, P. T. H. Ahlgren, and M. H. Jensen, Proc. Nal. Acad. Sci. 11, [5] R. Ton and A. Daffershofer, NeuroImage 16, doi: 1.116/j.neuroimage [51] I. Karazas and S. Shreve, Brownian Moion and Sochasic Calculus Springer, New York, 1. [5] P. Hänggi and H. Thomas, Phys. Rep. 88, [53] Y. L. Klimonovich, Phys. Usp. 37, [54] G. Pesce, A. McDaniel, S. Hoovy, J. Wehr, and G. Volpe, Na. Commun. 4, [55] B. C. dos Sanos and C. Tsallis, Phys. Rev. E 8, [56] I. M. Sokolov, Chem. Phys. 375, [57] N. G. van Kampen, J. Sa. Phys. 4, [58] M. Heidernäsch, On he diffusion in inhomogeneous sysems, Ph.D. hesis, Technische Universiä Chemniz, Faculy of Sciences, Insiue of Physics, Comple Sysems and Nonlinear Dynamics 15. [59] S. F. Kwok, Ann. Phys. 37, [6] Z. G. Arenas, D. G. Barci, and C. Tsallis, Phys. Rev. E 9, [61] M. Jeanblanc, M. Yor, and M. Chesney, Mahemaical Mehods for Financial Markes Springer, London, 9. [6] A. J. Bray, Phys. Rev. E 6, 13. [63] E. Marin, U. Behn, and G. Germano, Phys. Rev. E 83, [64] S. Karlin and H. M. Taylor, A Second Course in Sochasic Processes, nd ed. Academic, New York, [65] J. Ruseckas, R. Kazakevicius, and B. Kaulakys, J. Sa. Mech. 16, [66] A. G. Chersvy, A. V. Chechkin, and R. Mezler, J. Phys. A: Mah. Theor. 47, [67] A. G. Chersvy and R. Mezler, J. Sa. Mech. 15, P [68] J. M. Sancho, M. San Miguel, and D. Dürr, J. Sa. Phys., [69] P. D. Dilevsen, Geophys. Res. Le. 6, [7] P. G. Lind, A. Mora, J. A. C. Gallas, and M. Haase, Phys. Rev. E 7, [71] G. Volpe, L. Helden, T. Breschneider, J. Wehr, and C. Bechinger, Phys. Rev. Le. 14, [7] J. Smyhe, F. Moss, and P. V. E. McClinock, Phys. Rev. Le. 51, [73] K. S. Fa and E. K. Lenzi, Phys. Rev. E 67, [74] K. S. Fa, Phys. Rev. E 7, [75] H. Risken, The Fokker-Planck Equaion: Mehods of Soluion and Applicaions Springer-Verlag, Berlin, [76] J. Ruseckas and B. Kaulakys, Phys. Rev. E 81,