1
|
Bebon R, Godec A. Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime. PHYSICAL REVIEW LETTERS 2023; 131:237101. [PMID: 38134782 DOI: 10.1103/physrevlett.131.237101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 10/18/2023] [Accepted: 10/31/2023] [Indexed: 12/24/2023]
Abstract
We derive general bounds on the probability that the empirical first-passage time τ[over ¯]_{n}≡∑_{i=1}^{n}τ_{i}/n of a reversible ergodic Markov process inferred from a sample of n independent realizations deviates from the true mean first-passage time by more than any given amount in either direction. We construct nonasymptotic confidence intervals that hold in the elusive small-sample regime and thus fill the gap between asymptotic methods and the Bayesian approach that is known to be sensitive to prior belief and tends to underestimate uncertainty in the small-sample setting. We prove sharp bounds on extreme first-passage times that control uncertainty even in cases where the mean alone does not sufficiently characterize the statistics. Our concentration-of-measure-based results allow for model-free error control and reliable error estimation in kinetic inference, and are thus important for the analysis of experimental and simulation data in the presence of limited sampling.
Collapse
Affiliation(s)
- Rick Bebon
- Mathematical bioPhysics Group, Max Planck Institute for Multidisciplinary Sciences, 37077 Göttingen, Germany
| | - Aljaž Godec
- Mathematical bioPhysics Group, Max Planck Institute for Multidisciplinary Sciences, 37077 Göttingen, Germany
| |
Collapse
|
2
|
Singh P. Extreme value statistics and arcsine laws for heterogeneous diffusion processes. Phys Rev E 2022; 105:024113. [PMID: 35291128 DOI: 10.1103/physreve.105.024113] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Accepted: 01/26/2022] [Indexed: 06/14/2023]
Abstract
Heterogeneous diffusion with a spatially changing diffusion coefficient arises in many experimental systems such as protein dynamics in the cell cytoplasm, mobility of cajal bodies, and confined hard-sphere fluids. Here, we showcase a simple model of heterogeneous diffusion where the diffusion coefficient D(x) varies in a power-law way, i.e., D(x)∼|x|^{-α} with the exponent α>-1. This model is known to exhibit anomalous scaling of the mean-squared displacement (MSD) of the form ∼t^{2/2+α} and weak ergodicity breaking in the sense that ensemble averaged and time averaged MSDs do not converge. In this paper, we look at the extreme value statistics of this model and derive, for all α, the exact probability distributions of the maximum spatial displacement M(t) and arg-maximum t_{m}(t) (i.e., the time at which this maximum is reached) till duration t. In the second part of our paper, we analyze the statistical properties of the residence time t_{r}(t) and the last-passage time t_{ℓ}(t) and compute their distributions exactly for all values of α. Our study unravels that the heterogeneous version (α≠0) displays many rich and contrasting features compared to that of the standard Brownian motion (BM). For example, while for BM (α=0), the distributions of t_{m}(t),t_{r}(t), and t_{ℓ}(t) are all identical (á la "arcsine laws" due to Lévy), they turn out to be significantly different for nonzero α. Another interesting property of t_{r}(t) is the existence of a critical α (which we denote by α_{c}=-0.3182) such that the distribution exhibits a local maximum at t_{r}=t/2 for α<α_{c} whereas it has minima at t_{r}=t/2 for α≥α_{c}. The underlying reasoning for this difference hints at the very contrasting natures of the process for α≥α_{c} and α<α_{c} which we thoroughly examine in our paper. All our analytical results are backed by extensive numerical simulations.
Collapse
Affiliation(s)
- Prashant Singh
- International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bengaluru 560089, India
| |
Collapse
|
3
|
Singh P, Pal A. Extremal statistics for stochastic resetting systems. Phys Rev E 2021; 103:052119. [PMID: 34134348 DOI: 10.1103/physreve.103.052119] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2021] [Accepted: 04/28/2021] [Indexed: 11/07/2022]
Abstract
While averages and typical fluctuations often play a major role in understanding the behavior of a nonequilibrium system, this nonetheless is not always true. Rare events and large fluctuations are also pivotal when a thorough analysis of the system is being done. In this context, the statistics of extreme fluctuations in contrast to the average plays an important role, as has been discussed in fields ranging from statistical and mathematical physics to climate, finance, and ecology. Herein, we study extreme value statistics (EVS) of stochastic resetting systems, which have recently gained significant interest due to its ubiquitous and enriching applications in physics, chemistry, queuing theory, search processes, and computer science. We present a detailed analysis for the finite and large time statistics of extremals (maximum and arg-maximum, i.e., the time when the maximum is reached) of the spatial displacement in such system. In particular, we derive an exact renewal formula that relates the joint distribution of maximum and arg-maximum of the reset process to the statistical measures of the underlying process. Benchmarking our results for the maximum of a reset trajectory that pertain to the Gumbel class for large sample size, we show that the arg-maximum density attains a uniform distribution independent of the underlying process at a large observation time. This emerges as a manifestation of the renewal property of the resetting mechanism. The results are augmented with a wide spectrum of Markov and non-Markov stochastic processes under resetting, namely, simple diffusion, diffusion with drift, Ornstein-Uhlenbeck process, and random acceleration process in one dimension. Rigorous results are presented for the first two setups, while the latter two are supported with heuristic and numerical analysis.
Collapse
Affiliation(s)
- Prashant Singh
- International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bengaluru 560089, India
| | - Arnab Pal
- School of Chemistry, Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
4
|
Höll M, Wang W, Barkai E. Extreme value theory for constrained physical systems. Phys Rev E 2020; 102:042141. [PMID: 33212632 DOI: 10.1103/physreve.102.042141] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 10/12/2020] [Indexed: 06/11/2023]
Abstract
We investigate extreme value theory for physical systems with a global conservation law which describes renewal processes, mass transport models, and long-range interacting spin models. As shown previously, a special feature is that the distribution of the extreme value exhibits a nonanalytical point in the middle of the support. We expose exact relationships between constrained extreme value theory and well-known quantities of the underlying stochastic dynamics, all valid beyond the midpoint in general, i.e., even far from the thermodynamic limit. For example, for renewal processes the distribution of the maximum time between two renewal events is exactly related to the mean number of these events. In the thermodynamic limit, we show how our theory is suitable to describe typical and rare events which deviate from classical extreme value theory. For example, for the renewal process we unravel dual scaling of the extreme value distribution, pointing out two types of limiting laws: a normalizable scaling function for the typical statistics and a non-normalized state describing the rare events.
Collapse
Affiliation(s)
- Marc Höll
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Wanli Wang
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Eli Barkai
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
5
|
Bar A, Majumdar SN, Schehr G, Mukamel D. Exact extreme-value statistics at mixed-order transitions. Phys Rev E 2016; 93:052130. [PMID: 27300852 DOI: 10.1103/physreve.93.052130] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Indexed: 11/07/2022]
Abstract
We study extreme-value statistics for spatially extended models exhibiting mixed-order phase transitions (MOT). These are phase transitions that exhibit features common to both first-order (discontinuity of the order parameter) and second-order (diverging correlation length) transitions. We consider here the truncated inverse distance squared Ising model, which is a prototypical model exhibiting MOT, and study analytically the extreme-value statistics of the domain lengths The lengths of the domains are identically distributed random variables except for the global constraint that their sum equals the total system size L. In addition, the number of such domains is also a fluctuating variable, and not fixed. In the paramagnetic phase, we show that the distribution of the largest domain length l_{max} converges, in the large L limit, to a Gumbel distribution. However, at the critical point (for a certain range of parameters) and in the ferromagnetic phase, we show that the fluctuations of l_{max} are governed by novel distributions, which we compute exactly. Our main analytical results are verified by numerical simulations.
Collapse
Affiliation(s)
- Amir Bar
- Department of Complex Systems, Weizmann Institute, Rehovot, Israel
| | - Satya N Majumdar
- Université Paris-Sud, CNRS, LPTMS, UMR 8626, Orsay F-91405, France
| | - Grégory Schehr
- Université Paris-Sud, CNRS, LPTMS, UMR 8626, Orsay F-91405, France
| | - David Mukamel
- Department of Complex Systems, Weizmann Institute, Rehovot, Israel
| |
Collapse
|
6
|
Bologna M, Vanni F, Krokhin A, Grigolini P. Memory effects in fractional Brownian motion with Hurst exponent H<1/3. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:020102. [PMID: 20866763 DOI: 10.1103/physreve.82.020102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2010] [Indexed: 05/29/2023]
Abstract
We study the regression to the origin of a walker driven by dynamically generated fractional Brownian motion (FBM) and we prove that when the FBM scaling, i.e., the Hurst exponent H<1/3 , the emerging inverse power law is characterized by a power index that is a compelling signature of the infinitely extended memory of the system. Strong memory effects leads to the relation H=θ/2 between the Hurst exponent and the persistent exponent θ , which is different from the widely used relation H=1-θ . The latter is valid for 1/3<H<1 and is known to be compatible with the renewal assumption.
Collapse
Affiliation(s)
- Mauro Bologna
- Instituto de Alta Investigación, Universidad de Tarapacá, Chile
| | | | | | | |
Collapse
|
7
|
García-García R, Rosso A, Schehr G. Longest excursion of fractional Brownian motion: numerical evidence of non-Markovian effects. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:010102. [PMID: 20365309 DOI: 10.1103/physreve.81.010102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2009] [Indexed: 05/29/2023]
Abstract
We study, using exact numerical simulations, the statistics of the longest excursion l(max)(t) up to time t for the fractional Brownian motion with Hurst exponent 0<H<1. We show that in the large t limit, <l(max)(t)> proportional to variantQ(infinity)t, where Q(infinity) identical with Q(infinity)(H) depends continuously on H. These results are compared with exact analytical results for a renewal process with an associated persistence exponent theta=1-H. This comparison shows that Q(infinity)(H) carries the clear signature of non-Markovian effects for H not equal 1/2. The preasymptotic behavior of <l(max)(t)> is also discussed.
Collapse
|
8
|
Soriano J, Braslavsky I, Xu D, Krichevsky O, Stavans J. Universality of persistence exponents in two-dimensional Ostwald ripening. PHYSICAL REVIEW LETTERS 2009; 103:226101. [PMID: 20366113 DOI: 10.1103/physrevlett.103.226101] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2009] [Indexed: 05/29/2023]
Abstract
We measured persistence exponents theta(phi) of Ostwald ripening in two dimensions, as a function of the area fraction phi occupied by coarsening domains. The values of theta(phi) in two systems, succinonitrile and brine, quenched to their liquid-solid coexistence region, compare well with one another, providing compelling evidence for the universality of the one-parameter family of exponents. For small phi, theta(phi) approximately = 0.39phi, as predicted by a model that assumes no correlations between evolving domains. These constitute the first measurements of persistence exponents in the case of phase transitions with a conserved order parameter.
Collapse
Affiliation(s)
- Jordi Soriano
- Departament d'ECM, Facultat de Física, Universitat de Barcelona, Diagonal 647, E-08028 Barcelona, Spain
| | | | | | | | | |
Collapse
|