TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averages and their statistical variation for the Ornstein-Uhlenbeck process BT - Role of initial particle distributions and relaxation to stationarity JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and outof-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022134 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Thapa, Samudrajit A1 - Seno, Flavio A1 - Liu, Xianbin A1 - Metzler, Ralf T1 - Fractional Brownian motion with random diffusivity BT - emerging residual nonergodicity below the correlation time JF - Journal of physics : A, Mathematical and theoretical N2 - Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments. KW - stochastic processes KW - anomalous diffusion KW - fractional Brownian motion KW - diffusing diffusivity KW - weak ergodicity breaking Y1 - 2020 U6 - https://doi.org/10.1088/1751-8121/aba467 SN - 1751-8113 SN - 1751-8121 VL - 53 IS - 47 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes JF - New journal of physics : the open-access journal for physics N2 - We demonstrate the non-ergodicity of a simple Markovian stochastic process with space-dependent diffusion coefficient D(x). For power-law forms D(x) similar or equal to vertical bar x vertical bar(alpha), this process yields anomalous diffusion of the form < x(2)(t)> similar or equal to t(2/(2-alpha)). Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time-averaged mean-squared displacement <(delta(2)(Delta))over bar> remains linear in the lag time Delta and thus differs from the corresponding ensemble average < x(2)(t)>. We analyse the non-ergodic behaviour of this process in terms of the time-averaged mean- squared displacement (delta(2)) over bar and its random features, i.e. the statistical distribution of (delta(2)) over bar and the ergodicity breaking parameters. The heterogeneous diffusion model represents an alternative approach to non- ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media. Y1 - 2013 U6 - https://doi.org/10.1088/1367-2630/15/8/083039 SN - 1367-2630 VL - 15 IS - 15 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ultraslow scaled Brownian motion JF - New journal of physics : the open-access journal for physics N2 - We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form D(t) similar or equal to 1/t. For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations. KW - anomalous diffusion KW - stochastic processes KW - ageing Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/6/063038 SN - 1367-2630 VL - 17 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Ageing and confinement in non-ergodic heterogeneous diffusion processes JF - Journal of physics : A, Mathematical and theoretical N2 - We study the effects of ageing-the time delay between initiation of the physical process at t = 0 and start of observation at some time t(a) > 0-and spatial confinement on the properties of heterogeneous diffusion processes (HDPs) with deterministic power-law space-dependent diffusivities, D(x) = D-0 vertical bar x vertical bar(alpha). From analysis of the ensemble and time averaged mean squared displacements and the ergodicity breaking parameter quantifying the inherent degree of irreproducibility of individual realizations of the HDP we obtain striking similarities to ageing subdiffusive continuous time random walks with scale-free waiting time distributions. We also explore how both processes can be distinguished. For confined HDPs we study the long-time saturation of the ensemble and time averaged particle displacements as well as the magnitude of the inherent scatter of time averaged displacements and contrast the outcomes to the results known for other anomalous diffusion processes under confinement. KW - stochastic processes KW - anomalous diffusion KW - ageing KW - weak ergodicity breaking Y1 - 2014 U6 - https://doi.org/10.1088/1751-8113/47/48/485002 SN - 1751-8113 SN - 1751-8121 VL - 47 IS - 48 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity JF - Soft matter N2 - We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells. Y1 - 2014 U6 - https://doi.org/10.1039/c3sm52846d SN - 1744-683X SN - 1744-6848 VL - 10 IS - 10 SP - 1591 EP - 1601 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna S. T1 - Ultraslow scaled Brownian motion JF - New journal of physics : the open-access journal for physics N2 - We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations. KW - anomalous diffusion KW - stochastic processes KW - ageing Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/6/063038 SN - 1367-2630 VL - 17 IS - 063038 PB - Dt. Physikalische Ges., IOP CY - Bad Honnef, London ER - TY - GEN A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna S. T1 - Ultraslow scaled Brownian motion T2 - New journal of physics : the open-access journal for physics N2 - We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 188 KW - anomalous diffusion KW - stochastic processes KW - ageing Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-78618 ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averaging, ageing and delay analysis of financial time series JF - New journal of physics N2 - We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics. KW - time averaging KW - diffusion KW - geometric Brownian motion KW - financial time series Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa7199 SN - 1367-2630 VL - 19 SP - 1 EP - 11 PB - IOP CY - London ER - TY - GEN A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averaging, ageing and delay analysis of financial time series N2 - We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 347 KW - diffusion KW - financial time series KW - geometric Brownian motion KW - time averaging Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-400541 ER - TY - JOUR A1 - Bodrova, Anna A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Quantifying non-ergodic dynamics of force-free granular gases JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Brownian motion is ergodic in the Boltzmann-Khinchin sense that long time averages of physical observables such as the mean squared displacement provide the same information as the corresponding ensemble average, even at out-of-equilibrium conditions. This property is the fundamental prerequisite for single particle tracking and its analysis in simple liquids. We study analytically and by event-driven molecular dynamics simulations the dynamics of force-free cooling granular gases and reveal a violation of ergodicity in this Boltzmann-Khinchin sense as well as distinct ageing of the system. Such granular gases comprise materials such as dilute gases of stones, sand, various types of powders, or large molecules, and their mixtures are ubiquitous in Nature and technology, in particular in Space. We treat-depending on the physical-chemical properties of the inter-particle interaction upon their pair collisions-both a constant and a velocity-dependent (viscoelastic) restitution coefficient epsilon. Moreover we compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behaviour of the ensemble mean squared displacement (MSD) and the velocity correlations in the limit of weak dissipation. Qualitatively, the reported non-ergodic behaviour is generic for granular gases with any realistic dependence of epsilon on the impact velocity of particles. Y1 - 2015 U6 - https://doi.org/10.1039/c5cp02824h SN - 1463-9076 SN - 1463-9084 VL - 17 IS - 34 SP - 21791 EP - 21798 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averaging, ageing and delay analysis of financial time series JF - New journal of physics : the open-access journal for physics N2 - We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black-Scholes-Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics. KW - time averaging KW - diffusion KW - geometric Brownian motion KW - financial time series Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa7199 SN - 1367-2630 VL - 19 SP - 135 EP - 147 PB - IOP Publ. Ltd. CY - Bristol ER - TY - GEN A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion BT - (non-)existence of the overdamped limit in anomalous diffusion N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 267 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-97158 ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion BT - (non-)existence of the overdamped limit in anomalous diffusion JF - Scientific reports N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. Y1 - 2016 U6 - https://doi.org/10.1038/srep30520 SN - 2045-2322 VL - 6 PB - Nature Publishing Group CY - London ER - TY - JOUR A1 - Safdari, Hadiseh A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Thiel, Felix A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Quantifying the non-ergodicity of scaled Brownian motion JF - Journal of physics : A, Mathematical and theoretical N2 - We examine the non-ergodic properties of scaled Brownian motion (SBM), a non-stationary stochastic process with a time dependent diffusivity of the form D(t) similar or equal to t(alpha-1). We compute the ergodicity breaking parameter EB in the entire range of scaling exponents a, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths T and short lag times Delta the EB parameter as function of the scaling exponent a has no divergence at alpha - 1/2 and present the asymptotes for EB in different limits. We generalize the analytical and simulations results for the time averaged and ergodic properties of SBM in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion. KW - scaled Brownian motion KW - anomalous diffusion KW - ageing Y1 - 2015 U6 - https://doi.org/10.1088/1751-8113/48/37/375002 SN - 1751-8113 SN - 1751-8121 VL - 48 IS - 37 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion JF - Scientific reports N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. Y1 - 2016 U6 - https://doi.org/10.1038/srep30520 SN - 2045-2322 VL - 6 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Safdari, Hadiseh A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna A1 - Metzler, Ralf T1 - Aging underdamped scaled Brownian motion BT - Ensemble- and time-averaged particle displacements, nonergodicity, and the failure of the overdamping approximation JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We investigate both analytically and by computer simulations the ensemble- and time-averaged, nonergodic, and aging properties of massive particles diffusing in a medium with a time dependent diffusivity. We call this stochastic diffusion process the (aging) underdamped scaled Brownian motion (UDSBM). We demonstrate how the mean squared displacement (MSD) and the time-averaged MSD of UDSBM are affected by the inertial term in the Langevin equation, both at short, intermediate, and even long diffusion times. In particular, we quantify the ballistic regime for the MSD and the time-averaged MSD as well as the spread of individual time-averaged MSD trajectories. One of the main effects we observe is that, both for the MSD and the time-averaged MSD, for superdiffusive UDSBM the ballistic regime is much shorter than for ordinary Brownian motion. In contrast, for subdiffusive UDSBM, the ballistic region extends to much longer diffusion times. Therefore, particular care needs to be taken under what conditions the overdamped limit indeed provides a correct description, even in the long time limit. We also analyze to what extent ergodicity in the Boltzmann-Khinchin sense in this nonstationary system is broken, both for subdiffusive and superdiffusive UDSBM. Finally, the limiting case of ultraslow UDSBM is considered, with a mixed logarithmic and power-law dependence of the ensemble-and time-averaged MSDs of the particles. In the limit of strong aging, remarkably, the ordinary UDSBM and the ultraslow UDSBM behave similarly in the short time ballistic limit. The approaches developed here open ways for considering other stochastic processes under physically important conditions when a finite particle mass and aging in the system cannot be neglected. Y1 - 2017 U6 - https://doi.org/10.1103/PhysRevE.95.012120 SN - 2470-0045 SN - 2470-0053 VL - 95 PB - American Physical Society CY - College Park ER -