over the last fifteen years there has been a surge of studies in anomalous diffusion , characterized by the deviation of the mean squared displacement ( msd ) @xmath5 of a stochastic process with the probability density function @xmath1 to find the particle at position @xmath6 at time @xmath7 , from the linear time dependence @xmath8 of ordinary brownian motion @xcite .
anomalous diffusion is usually characterized in terms of the power - law form @xmath9 with the anomalous diffusion coefficient @xmath10 of physical dimension @xmath11 and the anomalous diffusion exponent @xmath12 .
depending on the value of @xmath12 we distinguish subdiffusion ( @xmath13 ) and superdiffusion ( @xmath14 ) @xcite .    since the milestone discoveries of superdiffusion in turbulence already in 1926 @xcite and of subdiffusion in amorphous semiconductors @xcite the recent vast increase of interest in anomalous diffusion is due to its discovery in numerous microscopic systems , in particular in biological contexts .
the cytoplasm of biological cells is heavily crowded with various obstacles , including proteins , nucleic acids , ribosomes , the cytoskeleton , as well as internal membranes compartmentalizing the cell @xcite .
diffusion of natural and artificial tracers in this complex environment is often subdiffusive .
similar situations are encountered in cell membranes @xcite .
the experimental evidence for subdiffusion in the crowded cytoplasm of living cells ranges from the motion of small labeled proteins @xcite over mrna molecules and chromosomal loci @xcite , lipid and insulin granules @xcite , virus particles @xcite , as well as chromosomal telomeres @xcite and cajal bodies @xcite inside the nucleus .
subdiffusion was observed for membrane resident proteins experimentally @xcite and for membrane lipid molecules in computer simulations @xcite . in
controlled _ in vitro _
experiments with artificial crowders , anomalous diffusion was consistently observed @xcite . on larger scales anomalous diffusion
was observed , for instance , for the motion of bacteria in a biofilm @xcite .
the observed subdiffusion was ascribed to various physical mechanisms @xcite .
apart from the apparent transient anomalous diffusion caused by a crossover from free normal diffusion to the plateau value of the msd @xcite , typically , three main families of anomalous diffusion processes are considered : ( i ) diffusion in a fractal environment where dead ends and bottlenecks slow down the motion on all scales @xcite ; ( ii ) motion in a viscoelastic environment , in which the effective anomalous motion of the tracer particle in the correlated many - body environment shows long - ranged antipersistent motion @xcite .
the latter process is associated with fractional brownian motion ( fbm ) and generalized langevin equation motion with a power - law memory form of the friction kernel @xcite .
( iii ) and continuous time random walk ( ctrw ) models , in which the moving particle is successively trapped by binding events to the environment or caging effects for waiting times @xmath15 distributed like a power - law @xmath16 with @xmath13 @xcite .
all three mechanisms lead to the power - law msd ( [ msd ] ) and they were indeed identified as processes generating the motion of different tracers in different cellular environments @xcite or in colloidal systems _ in vitro _ @xcite .
these three anomalous diffusion mechanisms characterized by anomalous diffusion with a constant in time generalized diffusivity were identified in experiments involving fairly large endogenous as well as artificial tracers .
recent experiments on eukaryotic cells @xcite using considerably smaller tracer proteins sampling over much larger subvolumes of the cell indicated a systematic variation of the cytoplasmic diffusivity with the separation from the cell nucleus .
these spatial diffusivity gradients are due partly to the non - uniform distribution of crowders in the cytoplasm .
this distribution can non - trivially affect the diffusion of tracers of different sizes in the cell cytoplasm . _ in vitro _ , fast gradients of the diffusivity can be realized , for instance , via a local variation of the temperature in thermophoresis experiments @xcite , as sketched in fig .
[ fig - model ] .
the diffusion of brownian particles in explicit solvents with temperature gradients was recently studied by multi - particle collision dynamics @xcite .
fluctuation - dissipation relations and spurious drift effects in systems with spatially varying friction coefficient were recently treated theoretically in ref .
@xcite . on larger scales ,
the diffusion of water molecules monitored by diffusive magnetic resonance imaging in the brain white matter was demonstrated to be heterogeneous and strongly anisotropic @xcite .
the anisotropy is due to the presence of some spatially - oriented structures in the tissue and obstacles which give rise to a tensorial character of the apparent diffusion coefficient .
the existence of a population splitting into two pools of water molecules with slow and fast diffusivities was shown for brain white matter @xcite .
finally , spatial heterogeneities are also abundant in the completely different context of anomalous diffusion in subsurface hydrology @xcite .    .
the diffusivity of the random walker is coupled to the spatial variation of the temperature of the environment , as symbolized by the speed of the runner.,width=302 ]    in the present study we examine the effects of the strength of the diffusivity gradient as defined by the scaling exponent @xmath4 and of the initial particle position @xmath17 in heterogeneous diffusion processes ( hdps ) with space - dependent diffusion coefficient @xmath18 .
we pay particular attention to a phenomenon , that recently received considerable attention for its immediate relevance to the surging field of single particle tracking experiments in microscopic systems , namely , the so - called _ weak ergodicity breaking_. this is the distinct disparity between physical observables depending on whether they are evaluated in the conventional ensemble sense or , from measured time series @xmath19 of the particle position , in terms of time averages @xcite .
we find that the hdp process gives rise to weakly non - ergodic behavior with a pronounced amplitude scatter of the time averaged msd of individual trajectories .
we obtain details of the distribution of this scatter as well as the frequently used ergodicity breaking and non - gaussianity parameters .
our analysis uncovers remarkable similarities of the hdp process with those of ctrw motion .
in particular , we study the behavior of the hdp at the critical value of the scaling exponent of the diffusivity .
the paper is organized as follows .
we introduce the hdp model with @xmath6-dependent diffusivity in sec .
[ sec - model ] .
the main results for the evolution of the msd , the time averaged msd , the probability density function , the ergodicity breaking and non - gaussianity parameters in the whole range of the model parameters are then presented in sec .
[ sec - main - results ] .
we discuss our results and point out the directions for future research in sec .
[ sec - discussion ] .
hdps are defined in terms of the multiplicative yet markovian langevin equation @xcite @xmath20 where @xmath18 is the position - dependent diffusion coefficient and @xmath21 represents white gaussian noise . in what follows we concentrate on the power - law form @xmath22 for the diffusivity , where the amplitude @xmath23 has dimension @xmath24 .
logarithmic and exponential forms for @xmath18 were considered in ref .
the offset @xmath25 in eq .
( [ eq - dx ] ) avoids either divergencies of @xmath18 ( @xmath26 ) or stalling of the particle ( @xmath27 ) around @xmath28 in the simulations . in the following calculations we use the scaling form @xmath29 .
we interpret the langevin equation ( [ lang ] ) in the stratonovich sense @xcite .
the msd following from the stochastic equation ( [ lang ] ) with diffusivity ( [ eq - dx ] ) takes on the power - law form @xcite @xmath30 where we introduced the scaling exponent @xmath31 which denotes superdiffusion for @xmath32 and subdiffusion for @xmath26 . for @xmath27
the diffusivity grows away from the origin , leading to a progressive acceleration of the particle as it ventures into more distant regions from the origin , and vice versa for @xmath26 . in the special case @xmath3 , the theory developed in ref .
@xcite breaks down as the scaling exponent ( [ pexp ] ) diverges . in refs .
@xcite an exponential growth of the msd was found , see the discussion below . for even larger values of @xmath4
the anomalous diffusion exponent @xmath33 becomes negative , i.e. , we observe a strong localization , see the discussion below .
the probability density function ( pdf ) of the hdp is given by the stretched or compressed gaussian @xcite @xmath34 for positive @xmath4 it exhibits a cusp , while for negative @xmath4 the pdf features a dip to zero at the origin .    in single particle tracking experiments one
measures the time series @xmath19 of the particle position for a time span @xmath35 .
it is usually evaluated in terms of the time averaged msd @xcite @xmath36 ^ 2dt , \label{eq - tamsd}\ ] ] where @xmath37 is the lag time .
brownian motion is ergodic and for sufficiently long measurement times @xmath35 the equality @xmath38 holds @xcite . anomalous
diffusion ( [ msd ] ) described by fbm or fractional langevin equation motion is asymptotically ergodic @xcite but may feature transiently non - ergodic behavior @xcite in confinement as well as transient aging @xcite , the explicit dependence on the time span elapsing between initial preparation of the system and start of the recording of the position , @xmath19 .
ctrw processes with diverging time scales of the waiting time distribution @xmath39 show weakly non - ergodic behavior for all @xmath37 .
namely , despite the scaling ( [ msd ] ) of the msd , @xmath40 scales _ linearly _ with @xmath37 .
more precisely , if we average over sufficiently many trajectories , the quantity @xmath41 for ctrws scales like @xmath42 @xcite .
this linear dependence on @xmath37 is preserved for aging ctrws @xcite .
apart from this linear @xmath37-dependence we also observe the dependence on the process time @xmath35 , a signature of aging @xcite .
interestingly , the hdp with diffusivity ( [ eq - dx ] ) displays weak ergodicity breaking of the form @xcite @xmath43 this behavior is analogous to that of scale - free ctrw motion , despite of the fact that the increment correlation function of hdps are ( anti)persistent in analogy to the ergodic fbm @xcite .
we note that also other processes such as aging and correlated ctrws @xcite as well as scaled brownian motion with time - dependent diffusivity @xcite exhibit the duality between ( [ msd ] ) and a linear @xmath37-dependence of @xmath40 .
( [ eq - tamsd - scaling ] ) shows the @xmath44-scaling as function of the process time @xmath35 .
for subdiffusion with @xmath45 , that is , the effective diffusivity of the process decays over time , as the particle ventures into low - diffusivity areas .
conversely , for @xmath46 the diffusivity increases over time , the particle discovers areas with increasingly higher @xmath18 .    in the current paper
we perform a detailed analysis of the msd and the time averaged msd in the entire range of @xmath4 , including the critical value @xmath3 .
we are particularly interested in the extent of the weakly non - ergodic behavior , especially the fluctuations of the time averaged msd around the mean value @xmath47 characterized by the ergodicity breaking parameter .
moreover we analyze the non - gaussianity of the process . in our analysis
we study effects of the scaling exponent @xmath4 of @xmath18 as well as the initial position of the particle .
the latter is known to affect the time - scales at which anomalous diffusion becomes significant @xcite .
we combine analytical and numerical approaches .
the simulations scheme for hdps was introduced in ref .
in this section we start with the analysis of the msd and the time averaged msd with its amplitude fluctuations , the latter being quantified by the corresponding scatter distribution .
we then analyze the ergodicity breaking and non - gaussianity parameters .
finally , we study the pdf of the function that quantifies the degree of particle dispersion .      fig .
[ fig - msd - tamsd - alfa ] shows the results from computer simulations for the msd @xmath48 and the time averaged msd @xmath40 from individual realizations along with the average @xmath49 taken over all single time traces .
the values for the scaling exponent @xmath4 of the diffusivity ( [ eq - dx ] ) studied in fig .
[ fig - msd - tamsd - alfa ] cover both the sub- and superdiffusive domains and include , in particular , the critical value @xmath3 where the scaling exponent @xmath33 of the msd ( [ msd ] ) diverges . for each @xmath4
we show @xmath50 single trajectories . in all cases ,
apart from the critical point , the scaling of the msd @xmath51 and both individual ( @xmath40 ) and mean ( @xmath47 ) time averaged msds agree well with the expected analytical behavior shown by the dashed lines : the scaling exponent of the msd varies consistently with @xmath4 , while that of the time averaged msd remains unity throughout .
the initial discrepancy between theory and msds is due to the choice for the initial position @xmath52 , whose influence relaxes on a time scale depending on both @xmath17 and @xmath4 .
the deviations of individual time traces @xmath40 at long lag times from the predicted behavior is due to unavoidable , bad statistics when the lag time gets close to the overall length @xmath35 of the time series .
irreproducibility of time averages of physical observables such as the msd is an intrinsic property of weakly non - ergodic processes @xcite . for ctrw processes with scale - free waiting time distribution
individual traces @xmath19 contain different , few dominating waiting time events that cause the amplitude scatter between different trajectories or , in other words , fluctuations of the apparent effective diffusion constant .
this phenomenon occurs no matter how long the measurement time @xmath35 is taken . for hdps
the scatter is due to the difference in the extent of excursions to regions of significantly different diffusivity .    the amplitude scatter between individual realizations @xmath40 in our hdp shown in fig .
[ fig - msd - tamsd - alfa ] varies significantly with the value of the scaling exponent @xmath4 : away from the brownian value @xmath53 in both subdiffusive ( @xmath26 ) and superdiffusive ( @xmath27 ) cases the fluctuations of @xmath40 become more pronounced when the modulus of @xmath4 increases . at the critical point
the fluctuations of @xmath40 appear to diverge , while beyond this critical point , the fluctuations decrease again .
moreover , a population splitting in a faster ( steeper slope of @xmath54 ) and slower ( shallower slope ) fraction of trajectories appears @xcite , see especially the panel for @xmath55 . in terms of the dimensionless variable
@xmath56 the amplitude scatter distribution @xmath57 reflects the randomness of individual time averages of the msd . for a sub- and a superdiffusive @xmath4
it was analyzed in ref .
as shown here for a whole spectrum of @xmath4 values , there is a clear trend towards extreme fluctuations at the critical point @xmath3 , but even for considerably smaller values such as @xmath58 the fluctuations around the mean @xmath59 are significant .
the width of the fluctuations of @xmath40 in each panel varies only marginally with the lag time @xmath37 , apart from the behavior at @xmath60 . as studied in ref .
@xcite , the relative amplitude scatter distribution @xmath57 can be fitted with the generalized gamma - distribution . in particular
, it tends to zero at @xmath61 , in contrast to subdiffusive ctrw processes , for which @xmath62 is always positive , indicating completely stalled trajectories @xcite .
we note that when the exponent @xmath4 approaches the critical value @xmath3 , the number of steps @xmath63 necessary to approach the theoretically predicted asymptote ( [ eq - msd - scaling ] ) increases significantly .
thus , as seen from fig .
[ fig - msd - tamsd - alfa ] , for @xmath52 only few simulation steps , @xmath64 are needed for negative @xmath4 with larger modulus .
it increases to @xmath65 steps for @xmath66 , @xmath67 for @xmath68 , and already @xmath69 for @xmath58 .
once the anomalous diffusion exponent @xmath33 becomes negative , that is , for values of @xmath4 larger than the critical value @xmath3 , the msd ( [ eq - msd - scaling ] ) becomes a decreasing function of time .
this agrees with our simulations if @xmath17 is chosen sufficiently large to enable the relaxation to the theoretical asymptote .
for instance , for the extreme value @xmath70 and values of the initial position of @xmath71 and above the msd indeed follows the theoretical prediction @xmath72 ( not shown ) . in this region
@xmath73 the diffusivity grows fast away from the origin the decreasing msd corresponds to the localization of particles in regions of slow diffusivity , see the detailed discussion in sec .
[ sec - spread ] below .    in the limit @xmath74 , due to the huge spread one or few extremely large amplitudes @xmath40 of the time averaged
msd may substantially affect the mean @xmath47 .
the msd in this limit follows an exponential growth @xmath75 , as indicated by the dashed curve in the panel for @xmath3 in fig .
[ fig - msd - tamsd - alfa ] .
such a fast increase of the msd is consistent with the divergence of the scaling exponent @xmath33 as function of @xmath4 in eq .
( [ eq - msd - scaling ] ) as well as with the exponential msd growth predicted for a parabolically space - varying diffusivity in refs .
this property can be straightforwardly inferred from the diffusion equation , @xmath76.\ ] ] multiplying both sides with @xmath77 and integrating over @xmath6 one arrives at @xcite @xmath78.\ ] ] this also rationalizes the observation that the scatter of @xmath40 is maximal for @xmath3 as the particles perform extremely far - reaching excursions relative to other @xmath4-values .      as introduced in refs .
@xcite the fluctuations of the time averaged msd @xmath40 can be quantified by their variance , the ergodicity breaking parameter @xmath79 when @xmath80 , it means that the process is perfectly reproducible and all time averages over sufficiently long trajectories give the same value .
this case corresponds to the sharp scatter distribution @xmath81 @xcite .
for the canonical brownian motion the ergodicity breaking parameter reaches the zero value in the form @xmath82 at finite ratio @xmath83 .
weakly non - ergodic processes have a positive value of @xmath84 . an alternative , weaker condition for ergodicity
is when the msd ( [ msd ] ) and the time averaged msd ( [ eatamsd ] ) coincide . to measure the relative deviations from ergodicity
, the @xmath85 parameter was introduced @xcite @xmath86 for ergodic dynamics its value is unity .     on the lag time @xmath37 .
the brownian asymptote @xmath87 is shown as the dashed curve in panel a ) .
the asymptotes ( [ eq - eb2 ] ) are shown in panel b ) as dashed lines of the corresponding color .
the initial position of the particle was @xmath52 and @xmath88 traces were used for each value of @xmath4 , all other parameters are the same as those used in fig .
[ fig - msd - tamsd - alfa ] .
the symbols in both panels correspond to the same parameters.,width=302 ]    , eq .  ( [ pexp ] ) , of @xmath89 in panels a ) and the magnification in panel b ) as well as @xmath90 in panel c ) .
the theoretical predictions from ref .
@xcite correspond to the dashed black curve in panels a ) and b ) . the ergodicity breaking parameter of the subdiffusive ctrw given by eq .  ( [ eq - eb - ctrw ] )
is shown by the dot - dashed red curve in panel b ) . to compute each point in the graphs for @xmath50 traces of @xmath91 steps
takes some two hours on a standard 3 ghz core workstation .
different sets of the shown points correspond to varying particle initial conditions @xmath17 indicated by different colors . note the alternative scale for the power exponent @xmath92 on the top axes .
parameters are the same as in fig .
[ fig - msd - tamsd - alfa].,width=264 ]    to calculate the @xmath84 parameter analytically is not always an easy task . to see this , note that the msd follows from the change of the stochastic variable to the standard wiener process in the form @xmath93^{-1/2}$ ]
for the time averaged msd @xmath40 the calculation is already more complicated , as it involves the two - point position correlation function .
the latter is expressed via fox @xmath94-functions for the hdp @xcite .
the analytical derivation of the ergodicity breaking parameter @xmath84 , however , requires fourth - order moments , whose calculation is a formidable task .
so far only approximate methods are known for hdps @xcite . for subdiffusive
ctrw processes it is possible to obtain the @xmath84 parameter more easily from the conjectured and numerically proven equivalence @xmath95 of @xmath96 with the ratio of the number of steps @xmath97 in an individual realization and the average @xmath98 @xcite . as for the limiting distribution for @xmath99
is known @xcite , this allows straightforward calculation of @xmath57 and its moments @xcite . for the multiplicative process
studied here such a scheme does not work .
we also note that another parameter involving the fourth moment of the particle displacement is the non - gaussianity parameter @xmath100 @xcite discussed below . due to
the lack of an analytical theory , a major reason for the current simulations study is to explore the behavior of these two parameters in the whole range of the model parameters .    figs .
[ fig - eb - lag ] and [ fig - eb - p ] display the dependence of the ergodicity breaking parameters on the lag time @xmath37 , the initial position @xmath17 , as well as the scaling exponent @xmath33 .
we find that for standard brownian motion with @xmath101 the ergodicity breaking parameter eb follows the known asymptote @xmath102 from above .
concurrently @xmath103 as expected for ergodic motion , except for very small @xmath37 values because of the initial relaxation of the influence of the initial position @xmath17 .
as we depart from the value @xmath53 of brownian motion , the magnitude of @xmath84 grows and its functional dependence on the lag time becomes less pronounced , see fig .
[ fig - eb - lag]a . approaching the critical point @xmath104 , due to the huge fluctuations of @xmath40 and the ensuing possibility of extreme events
the ergodicity breaking parameter also explicitly depends on the number of traces @xmath105 used for the averaging and reaches values of @xmath106 and higher . on both sides of the critical point , corresponding to large positive or negative values of @xmath33 , the values of the parameter @xmath107 approach one another , as seen in fig .
[ fig - eb - p]a , while @xmath108 exhibits a jump , fig .  [ fig - eb - p]b .
the ratio of the different msds given by @xmath85 also depends on the power exponent @xmath4 and the initial value @xmath17 . as anticipated already from fig .
[ fig - msd - tamsd - alfa ] , for the chosen initial condition @xmath52 the magnitude of @xmath85 decreases as @xmath4 gets progressively negative , while @xmath85 grows as @xmath4 increases from 0 to 2 , reaching very large values @xmath109 at @xmath3 .
this is illustrated in fig .
[ fig - eb - lag]b which is also consistent with the scaling @xmath110 with the lag time at different @xmath111 values .
because of the initial condition @xmath17 , similar to our statements for @xmath47 this asymptote is approached later when @xmath74 .    in fig .
[ fig - eb - p ] we also analyze the effects of the initial position @xmath17 and show the dependence of the ergodicity breaking parameters for long traces or short lag times , i.e. , when @xmath112 which is practically equivalent to @xmath113 , versus the anomalous diffusion exponent @xmath33 and the scaling exponent @xmath4 .
we observe that in the region @xmath46 the values of @xmath107 obtained from simulations is in good agreement with the analytical estimate @xmath113 from ref .
@xcite , represented by the dashed curve in fig .  [ fig - eb - p]a and b. these approximations are based on the asymptotic scaling of @xmath18 and do not consider the effect of the initial positions @xmath17 , assumed to be relaxed in the relevant @xmath114 limit . to see their actual impact on the dynamics we study the ergodicity breaking parameters numerically .
we observe that the value of @xmath107 grows from the small pre - asymptotic brownian value at @xmath115 to progressively larger values at larger modulus of @xmath33 , as @xmath4 approaches the critical value @xmath3 from below and above .
the functional dependence of @xmath116 obtained from our simulations follows quite well the predictions from the approximate calculation of @xmath84 @xcite , particularly in the range @xmath117 , see fig .
[ fig - eb - p]b .
the deviations for even more pronounced variation of @xmath18 as @xmath4 tends to 2 are likely due to insufficient statistics . as we show in fig .
[ fig - eb - p]a by the circles and triangles , the value @xmath118 reveals measurable deviations for @xmath119 as compared to @xmath50 traces used for the averaging , while for less extreme @xmath4 values the two sets yield nearly identical results .
this is a further indication towards a divergence of the fluctuations at the critical point .
namely , the chance to find an even larger amplitude of @xmath40 in a given realization grows with the number of simulated trajectories .    as already anticipated in ref .
@xcite the hdp process considered here is in some aspects reminiscent of subdiffusive ctrw processes with scale - free , power - law waiting time distributions with regard to the scaling of the msd and the time averaged msd , compare also sec .
[ sec - discussion ] .
apparently , the ergodicity breaking properties of these two processes also reveal similar features : the fluctuations of @xmath40 increase when the anomalous diffusion exponent @xmath33 becomes successively smaller than the brownian value 1 . indeed , as shown in fig .
[ fig - eb - p]b for the relevant region @xmath45 there exists even a close agreement of the quantitative behavior of ctrw and hdp .
we show the ergodicity breaking parameter obtained for the ctrw process @xcite @xmath120 where @xmath13 is the exponent in the pdf @xmath16 of the waiting times @xmath15 @xcite . for the ctrw process
the exponent @xmath12 also occurs in the msd ( [ msd ] ) .
we note that in contrast to hpds in superdiffusive ctrw processes the nonergodic behavior is merely ultraweak , effecting a different prefactor in the time averaged msd compared to the msd @xcite .    for @xmath121
we find that @xmath122 for initial positions @xmath17 close to the origin , in agreement with the results presented in ref .
@xcite . for @xmath73 or @xmath123 the analytical model of ref .
@xcite no longer applies but the simulations yield another region of growth for the value of @xmath107 . note that @xmath107 also depends on the trace length @xmath35 ( not shown ) , as discussed for two - dimensional hdps @xcite . as we show in fig .
[ fig - eb - p]c , the magnitude of @xmath108 strongly varies with @xmath17 . in the region
@xmath124 the ratio of the msds grows as @xmath125 indicated by the dashed asymptote in fig .
[ fig - eb - p]c , while for large negative @xmath33 we observe @xmath126 .
we finally note that @xmath108 reveals a much weaker dependence on the number @xmath105 of simulated traces , see the circles and triangles in fig .
[ fig - eb - p]c .
the reason is that @xmath85 involves only the second moment of the time averaged msd , while the @xmath84 is defined in terms of the fourth order moment , which is more sensitive to large variations .       as function of the lag time ,
plotted for the parameters of fig .
[ fig - eb - lag ] with @xmath50 trajectories for each @xmath4 value .
the approximately zero - valued trace @xmath127 for brownian motion ( @xmath53 ) is only shown partially .
color coding is the same as in fig .
[ fig - eb - lag ] .
* b ) * the value of @xmath128 evaluated for varying power exponent @xmath4 and initial position @xmath17 .
the color scheme is the same as in fig .
[ fig - eb - p].,width=317 ]    the non - gaussianity parameter @xmath100 is a sensitive experimental indicator that often enables one to distinguish the type of particular diffusion processes observed in single - particle tracking experiments @xcite .
it is related to the stationarity of increments of the diffusion process and involves the fourth moment along the time - averaged trajectory . for the diffusion process in an embedding space of dimension
@xmath129 it is defined via the experimentally relevant time averaged quantities as @xcite @xmath130 where , in analogy to eq .
( [ eq - tamsd ] ) , the fourth moment is defined via @xmath131 ^ 4~dt.\ ] ] for brownian and fractional brownian motion , both gaussian processes , one consistently finds @xmath132 . for diffusion processes revealing a transient anomalous diffusion behavior ,
this parameter deviates substantially from zero , see , e.g. , the discussion in ref .
@xcite .
we systematically examine the behavior of @xmath127 in fig .
[ fig - non - gauss]a . corroborating the results for the ergodicity breaking parameter , for brownian motion we obtain approximately zero values .
as the power @xmath4 in @xmath133 deviates from zero , the non - gaussianity parameter reveals a rich behavior as a function of the lag time @xmath37 and the initial position @xmath17 .
more specifically , in the region @xmath26 the non - gaussianity parameter progressively grows and reaches considerably large values for large @xmath2 , see the curve for @xmath134 in fig .
[ fig - non - gauss]a .
the value @xmath127 systematically decreases with the lag time @xmath37 along the trace for negative @xmath4 . in the region @xmath135
the value of @xmath100 also grows with @xmath4 but , in contrast to the case @xmath26 , the function @xmath127 stays rather constant with @xmath37 . as we approach the critical @xmath3 the non - gaussianity parameter reaches very high values . in the region @xmath73 the value of @xmath127 decreases again .
these features of the functional behavior of @xmath136 correspond with the properties of @xmath137 , compare the curves in fig .
[ fig - non - gauss]a and [ fig - eb - lag]a .
note that similar to @xmath84 the statistics required to obtain a smooth curve for @xmath127 strongly depends on @xmath4 : the result from @xmath50 traces in fig .
[ fig - non - gauss]a acquires pronouncedly higher fluctuations when we approach the critical value @xmath3 .
we find that variations in the initial positions @xmath17 only have a marginal effect on @xmath138 . in fig .
[ fig - non - gauss]b we illustrate the range of @xmath100 values for exponents in the range @xmath139 . even for the wide range from @xmath71 to @xmath140 the value of @xmath128
does not reveal any appreciable variation , in stark contrast to the strong @xmath17-dependence of @xmath141 .
this indicates that for the hdp the non - gaussianity parameter is more robust than @xmath84 with respect to the choice of the initial condition . in fig .
[ fig - non - gauss]b the dashed line indicates the proportionality to @xmath142 , which nicely matches the measured shape of @xmath138 .
we finally analyze the time evolution of the pdf of the particle position in fig .
[ fig - spread ] for three different values of the scaling exponent @xmath4 , distinguishing superdiffusion ( top row for @xmath68 , i.e. , ballistic motion with @xmath143 ) , normal diffusion ( middle row for @xmath53 and @xmath115 ) , and subdiffusion ( bottom row for @xmath144 , i.e. , @xmath145 ) .
in each row the leftmost panel represents a very early evolution of the pdf close to the initial condition , while in the rightmost panel the pdf in some of the cases has almost reached the diffusion limit , in which the analytical asymptote ( [ pdf - hdp ] ) is valid .
we observe that for hdps with positive @xmath4 , i.e. , when the diffusivity @xmath18 grows with the modulus of the position , the pdf for particles with an initial position far away from the origin , an asymmetric shape of the pdf is effected .
namely , they progressively move towards regions of small diffusivity and accumulate there .
this focusing due to the quenched nature of the diffusivity erases any memory of the initial condition and the common asymptote ( [ pdf - hdp ] ) of stretched gaussian shape ( @xmath135 ) is approached . in the opposite case with negative
@xmath4 the focusing of particles in lower diffusivity regions applies again , albeit at longer times than shown here .
this time , however , instead of the cusp at the origin for positive @xmath4 , the pdf drops down to zero at the origin and acquires the bimodal shape predicted by ( [ pdf - hdp ] ) , a compressed gaussian . for the brownian case with vanishing @xmath4 ,
no focusing takes place . in the final panel
the distribution is already so broad that the individual , shifted pdfs appear on top of each other .
similarly to the evolution of the msd shown in fig .
[ fig - msd - tamsd - alfa ] , the relaxation time required for the system of diffusing particles to approach this long - time limit , @xmath1 depends on the power @xmath4 of the diffusivity in eq .
( [ eq - dx ] ) . for more negative @xmath4 values
the traces are not yet relaxed to the theoretical long - time shape even after @xmath91 steps .
the particles starting at larger @xmath17 remain trapped for the whole length of the simulation .
conversely , when the diffusivity increases fast away from the origin ( moderate positive @xmath4 values ) , the equilibration to the long - time pdf is relatively fast , see the panel for @xmath68 in fig .
[ fig - spread ] .
for hpd processes whose diffusivity @xmath18 varies in power - law form with the distance @xmath6 from the origin and which give rise to anomalous diffusion we analyzed in detail the weakly non - ergodic behavior for varying power exponents @xmath4 and initial positions @xmath17 of the particles . in particular
, we examined the functional dependencies and magnitudes of the variation of the ensemble and time averaged characteristics of such hdps with these parameters .
the fluctuations of the amplitude of the time averaged msd of individual time traces were shown to systematically grow with increasing departure of the anomalous diffusion exponent @xmath33 from the brownian value @xmath115 , corresponding to the scaling exponent @xmath53 in the power - law form for @xmath18 .
the fluctuations increase dramatically when the scaling exponent approaches its critical value @xmath3 . at that critical point
we corroborated the turnover from the power - law scaling in time of the msd to an exponential growth .
beyond the critical point the fluctuations of the time averaged msd become smaller again .
concurrently we observe a distinct population splitting into a faster and slower fraction of time averaged msds .
we paid particular emphasis on several parameters used to classify the departure from ergodic behavior , namely the two ergodicity breaking parameters @xmath84 and @xmath85 as well as the non - gaussianity parameter @xmath100 .
while @xmath146 is simply defined as the ratio of the mean time averaged msd versus the msd , both @xmath84 and @xmath100 are based on the fourth order moments of the particle position and are known analytically only from approximate theories . a detailed numerical analysis of their properties was therefore used to obtain more concrete information on their behavior for different values of the scaling exponent @xmath4 and the initial particle position @xmath17 in the heterogeneous environment @xcite .
within the analyzed parameter range we find that the behavior of @xmath147 is indeed in agreement with the heuristic theoretical analysis from ref .
the parameter portraits for the ergodicity breaking and non - gaussianity parameters obtained from our numerical analysis will be useful for actual data evaluation
. we also explored the detailed behavior of the hdp dynamics at the critical value @xmath3 and its vicinity , in particular with respect to the dramatic values reached by @xmath84 reflecting the dramatic fluctuations of the time averaged msd .
a systematic numerical analysis of the particle pdf for varying initial conditions @xmath17 sheds additional light on hdps with different exponent @xmath4 , in particular , the visualization of the particle focusing in low diffusivity regions .    in hdps
the non - ergodic behavior arises due to the heterogeneity of the environment .
physically , this represents a space - dependent mobility or temperature .
it is not a property of the particle , and each time the particle revisits a given point @xmath6 in space it has the same diffusivity @xmath18 . in a random walk sense , this scenario could also be translated into a local dependence of the waiting time for a jump event . in that interpretation
the hdp corresponds to a motion in a quenched energy landscape @xcite , albeit a deterministic ( in contrast to random ) one .
as such , the process is intrinsically different from renewal ctrw processes , which correspond to the motion in an annealed environment @xcite .
nevertheless we observed that subdiffusive hdps share a number of features with subdiffusive ctrws with scale - free , power - law waiting time distribution .
these features include the scaling laws for the ensemble and time averaged msds , and the form of the ergodicity breaking parameter @xmath85 ( in the relevant range @xmath45 ) . from the sole analysis of the msds @xmath51 and @xmath40 as well as the ergodicity breaking parameters , a significant distinction between ctrw and hdp
is therefore difficult . yet there
exist some crucial differences between the hdp and ctrw processes , that can be measured .
thus , in hdps with @xmath45 the distribution @xmath57 of the relative amplitude @xmath40 of individual realizations decays to zero at @xmath61 , in contrast to the ctrw s finite fraction of immobile particles reflected in the finite value @xmath148 .
indeed the scatter distribution @xmath57 was previously advocated as a good diagnosis tool for different stochastic processes @xcite , complementing other stochastic analysis methods @xcite .
hdps are physically meaningful alternatives to the standard stochastic models for anomalous diffusion processes , namely , the ctrw process with long - tailed waiting time distribution , fbm and fractional langevin equation motion based on gaussian yet long - ranged in time correlations , diffusion in fractal environments , or their combinations @xcite .
hpds are weakly non - ergodic , sharing some yet not all properties of ctrw processes .
it is therefore important to explore the properties of hdps in more detail in the future as well as to develop more sophisticated tools to interpret measured time series and with confidence identify the underlying stochastic mechanism . in particular
, we will study the properties of confined and aged hdps as well as spatially varying diffusivities with a random component , and the coupling of hdps with active motion .
i. bronstein , y. israel , e. kepten , s. mai , y. shav - tal , e. barkai , and y. garini , phys .
lett . * 103 * , 018102 ( 2009 ) ; k. burnecki , e. kepten , j. janczura , i. bronshtein , y. garini , and a. weron , biophys .
j. * 103 * , 1839 ( 2012 ) .
g. r. kneller , k. baczynski , and m. pasenkiewicz - gierula , j. chem . phys . * 135 * , 141105 ( 2011 ) ; j .-
jeon , h. martinez - seara monne , m. javanainen , and r. metzler , phys .
* 109 * , 188103 ( 2012 ) ; m. javanainen , h. hammaren , l. monticelly , j .- h . jeon , m. s. miettinen , h. martinez - seara monne , r. metzler , and i. vattulainen , faraday disc .
* 161 * , 397 ( 2013 ) .
s. c. kou , ann .
* 2 * , 501 ( 2008 ) ; e. lutz , phys .
e * 64 * , 051106 ( 2001 ) ; r. kupferman , j. stat .
phys . * 114 * , 291 ( 2004 ) ; w. min , g. luo , b. j. cherayil , s. c. kou , and x. s. xie , phys .
lett . * 94 * , 198302 ( 2005 ) .      q. xu , l. feng , r. sha , n. c. seeman , and p. m. chaikin , phys .
106 * , 228102 ( 2011 ) ; i. y. wong , m. l. gardel , d. r. reichman , e. r. weeks , m. t. valentine , a. r. bausch , and d. a. weitz , _ ibid . _ * 92 * , 178101 ( 2004 ) .
bouchaud , j. phys .
i ( paris ) * 2 * , 1705 ( 1992 ) ; g. bel and e. barkai , phys .
lett . * 94 * , 240602 ( 2005 ) ; a. rebenshtok and e. barkai , _ ibid . _ * 99 * , 210601 ( 2007 ) ; m. a. lomholt , i. m. zaid , and r. metzler , phys .
lett . * 98 * , 200603 ( 2007 ) .
v. tejedor and r. metzler , j. phys .
a * 43 * , 082002 ( 2010 ) ; m. magdziarz , r. metzler , w. szczotka , and p. zebrowski , phys . rev .
e * 85 * , 051103 ( 2012 ) .
see also the alternative approach in a. v. chechkin , m. hofmann , and i. m. sokolov , phys .
e * 80 * , 031112 ( 2009 ) ; j. h. p. schulz , a. v. chechkin , and r. metzler , j. phys .
a. * 46 * , 475001 ( 2013 ) .
g. zumofen and j. klafter , physica d * 69 * , 436 ( 1993 ) ; a. godec and r. metzler , phys .
e * 88 * , 012116 ( 2013 ) ; d. froemberg and e. barkai , _ ibid .
_ * 87 * , 030104(r ) ( 2013 ) ; _ ibid . _ * 88 * , 024101 ( 2013 ) ; t. akimoto , phys .
* 108 * , 164101 ( 2012 ) .
m. magdziarz , a. weron , k. burnecki , and j. klafter , phys .
. lett . * 103 * , 180602 ( 2009 ) ; m. magdziarz and j. klafter , phys .
e * 82 * , 011129 ( 2010 ) ; y. meroz , i. eliazar , and j. klafter , j. phys .
a * 42 * , 434012 ( 2009 ) ; y. meroz , i. m. sokolov , and j. klafter , phys . rev . lett . * 110 * , 090601 ( 2013 ) ; s. condamin , v. tejedor , r. voituriez , o. bnichou , and j. klafter , proc .
usa * 105 * , 5675 ( 2008 ) ; v. tejedor , o. bnichou , r. voituriez , r. jungmann , f. simmel , c. selhuber - unkel , l. oddershede , and r. metzler , biophys .
j. * 98 * , 1364 ( 2010 ) ; m. bauer , r. valiullin , g. radons , and j. kaerger , j. chem .
phys . * 135 * , 144118 ( 2011 ) ; j .- h .
jeon , e. barkai , and r. metzler , j. chem . phys .
* 139 * , 121916 ( 2013 ) .