1
|
De Bruyne B, Majumdar SN, Schehr G. Generating discrete-time constrained random walks and Lévy flights. Phys Rev E 2021; 104:024117. [PMID: 34525540 DOI: 10.1103/physreve.104.024117] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/27/2021] [Indexed: 11/07/2022]
Abstract
We introduce a method to exactly generate bridge trajectories for discrete-time random walks, with arbitrary jump distributions, that are constrained to initially start at the origin and return to the origin after a fixed time. The method is based on an effective jump distribution that implicitly accounts for the bridge constraint. It is illustrated on various jump distributions and is shown to be very efficient in practice. In addition, we show how to generalize the method to other types of constrained random walks such as generalized bridges, excursions, and meanders.
Collapse
Affiliation(s)
- Benjamin De Bruyne
- LPTMS, CNRS, Université Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| | - Satya N Majumdar
- LPTMS, CNRS, Université Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| | - Grégory Schehr
- Sorbonne Université, Laboratoire de Physique Théorique et Hautes Energies, CNRS UMR 7589, 4 Place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
2
|
Aliakbari A, Manshour P, Salehi MJ. Records in fractal stochastic processes. CHAOS (WOODBURY, N.Y.) 2017; 27:033116. [PMID: 28364750 DOI: 10.1063/1.4979348] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The record statistics in stationary and non-stationary fractal time series is studied extensively. By calculating various concepts in record dynamics, we find some interesting results. In stationary fractional Gaussian noises, we observe a universal behavior for the whole range of Hurst exponents. However, for non-stationary fractional Brownian motions, the record dynamics is crucially dependent on the memory, which plays the role of a non-stationarity index, here. Indeed, the deviation from the results of the stationary case increases by increasing the Hurst exponent in fractional Brownian motions. We demonstrate that the memory governs the dynamics of the records as long as it causes non-stationarity in fractal stochastic processes; otherwise, it has no impact on the record statistics.
Collapse
Affiliation(s)
- A Aliakbari
- Department of Physics, Faculty of Sciences, Persian Gulf University, 75169 Bushehr, Iran
| | - P Manshour
- Department of Physics, Faculty of Sciences, Persian Gulf University, 75169 Bushehr, Iran
| | - M J Salehi
- Department of Physics, Faculty of Sciences, Persian Gulf University, 75169 Bushehr, Iran
| |
Collapse
|
3
|
Perret A, Comtet A, Majumdar SN, Schehr G. Near-extreme statistics of Brownian motion. PHYSICAL REVIEW LETTERS 2013; 111:240601. [PMID: 24483638 DOI: 10.1103/physrevlett.111.240601] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Revised: 09/18/2013] [Indexed: 06/03/2023]
Abstract
We study the statistics of near-extreme events of Brownian motion (BM) on the time interval [0,t]. We focus on the density of states near the maximum, ρ(r,t), which is the amount of time spent by the process at a distance r from the maximum. We develop a path integral approach to study functionals of the maximum of BM, which allows us to study the full probability density function of ρ(r,t) and obtain an explicit expression for the moments <[ρ(r,t)]k> for arbitrary integer k. We also study near extremes of constrained BM, like the Brownian bridge. Finally we also present numerical simulations to check our analytical results.
Collapse
Affiliation(s)
- Anthony Perret
- Université Paris-Sud-Paris 11, CNRS, LPTMS, 91405 Orsay Cedex, France
| | - Alain Comtet
- Université Paris-Sud-Paris 11, CNRS, LPTMS, 91405 Orsay Cedex, France and Université Pierre et Marie Curie-Paris 6, 75252 Paris Cedex 05, France
| | - Satya N Majumdar
- Université Paris-Sud-Paris 11, CNRS, LPTMS, 91405 Orsay Cedex, France
| | - Grégory Schehr
- Université Paris-Sud-Paris 11, CNRS, LPTMS, 91405 Orsay Cedex, France
| |
Collapse
|
4
|
LeBlanc M, Angheluta L, Dahmen K, Goldenfeld N. Universal fluctuations and extreme statistics of avalanches near the depinning transition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:022126. [PMID: 23496478 DOI: 10.1103/physreve.87.022126] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2012] [Indexed: 06/01/2023]
Abstract
We derive exact predictions for universal scaling exponents and scaling functions associated with the statistics of maximum velocities v(m) during avalanches described by the mean-field theory of the interface depinning transition. In particular, we find a robust power-law regime in the statistics of maximum events that can explain the observed distribution of the peak amplitudes in acoustic emission experiments of crystal plasticity. Our results are expected to be broadly applicable to a broad range of systems in the mean-field interface depinning universality class, ranging from magnets to earthquakes.
Collapse
Affiliation(s)
- Michael LeBlanc
- Department of Physics, University of Illinois at Urbana-Champaign, Loomis Laboratory of Physics, 1110 West Green Street, Urbana, Illinois 61801-3080, USA
| | | | | | | |
Collapse
|
5
|
Das J. Positive feedback produces broad distributions in maximum activation attained within a narrow time window in stochastic biochemical reactions. J Chem Phys 2013; 138:015101. [PMID: 23298061 DOI: 10.1063/1.4772583] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
How do single cell fate decisions induced by activation of key signaling proteins above threshold concentrations within a time interval are affected by stochastic fluctuations in biochemical reactions? We address this question using minimal models of stochastic chemical reactions commonly found in cell signaling and gene regulatory systems. Employing exact solutions and semi-analytical methods we calculate distributions of the maximum value (N) of activated species concentrations (P(max)(N)) and the time (t) taken to reach the maximum value (P(max)(t)) within a time interval in the minimal models. We find, the presence of positive feedback interactions make P(max)(N) more spread out with a higher "peakedness" in P(max)(t). Thus positive feedback interactions may help single cells to respond sensitively to a stimulus when cell decision processes require upregulation of activated forms of key proteins to a threshold number within a time window.
Collapse
Affiliation(s)
- Jayajit Das
- Battelle Center For Mathematical Medicine, The Research Institute at the Nationwide Childrens Hospital and Department of Pediatrics, and Biophysics Graduate Program, The Ohio State University, 700 Childrens Drive, Columbus, Ohio 43205, USA.
| |
Collapse
|
6
|
Huang C, Lin JG, Ren YY. Testing for the shape parameter of generalized extreme value distribution based on the $$L_q$$ -likelihood ratio statistic. METRIKA 2012. [DOI: 10.1007/s00184-012-0409-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
7
|
LeBlanc M, Angheluta L, Dahmen K, Goldenfeld N. Distribution of maximum velocities in avalanches near the depinning transition. PHYSICAL REVIEW LETTERS 2012; 109:105702. [PMID: 23005300 DOI: 10.1103/physrevlett.109.105702] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/14/2012] [Indexed: 06/01/2023]
Abstract
We report exact predictions for universal scaling exponents and scaling functions associated with the distribution of the maximum collective avalanche propagation velocities v(m) in the mean field theory of the interface depinning transition. We derive the extreme value distribution P(v(m)|T) for the maximum velocities in avalanches of fixed duration T and verify the results by numerical simulation near the critical point. We find that the tail of the distribution of maximum velocity for an arbitrary avalanche duration, v(m), scales as P(v(m))~v(m)(-2) for large v(m). These results account for the observed power-law distribution of the maximum amplitudes in acoustic emission experiments of crystal plasticity and are also broadly applicable to other systems in the mean-field interface depinning universality class, ranging from magnets to earthquakes.
Collapse
Affiliation(s)
- Michael LeBlanc
- Department of Physics, Loomis Laboratory of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-3080, USA
| | | | | | | |
Collapse
|
8
|
Schehr G, Majumdar SN. Universal order statistics of random walks. PHYSICAL REVIEW LETTERS 2012; 108:040601. [PMID: 22400820 DOI: 10.1103/physrevlett.108.040601] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2011] [Indexed: 05/31/2023]
Abstract
We study analytically the order statistics of a time series generated by the positions of a symmetric random walk of n steps with step lengths of finite variance σ(2). We show that the statistics of the gap d(k,n) = M(k,n)-M(k+1,n) between the kth and the (k+1)th maximum of the time series becomes stationary, i.e., independent of n as n → ∞ and exhibits a rich, universal behavior. The mean stationary gap exhibits a universal algebraic decay for large k, ~d(k,∞)-/σ 1/sqrt[2πk], independent of the details of the jump distribution. Moreover, the probability density (pdf) of the stationary gap exhibits scaling, Pr(d(k,∞) = δ) ~/= (sqrt[k]/σ)P(δsqrt[k]/σ), in the regime δ~ (d(k,∞)). The scaling function P(x) is universal and has an unexpected power law tail, P(x) ~ x(-4) for large x. For δ>> (d(k,∞)) the scaling breaks down and the pdf gets cut off in a nonuniversal way. Consequently, the moments of the gap exhibit an unusual multiscaling behavior.
Collapse
Affiliation(s)
- Grégory Schehr
- Université Paris-Sud, CNRS, LPT, 91405 Orsay Cedex, France.
| | | |
Collapse
|
9
|
Cherstvy AG. Critical polyelectrolyte adsorption under confinement: planar slit, cylindrical pore, and spherical cavity. Biopolymers 2012; 97:311-7. [PMID: 22241107 DOI: 10.1002/bip.22023] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2011] [Revised: 12/02/2011] [Accepted: 12/21/2011] [Indexed: 01/22/2023]
Abstract
We explore the properties of adsorption of flexible polyelectrolyte chains in confined spaces between the oppositely charged surfaces in three basic geometries. A method of approximate uniformly valid solutions for the Green function equation for the eigenfunctions of polymer density distributions is developed to rationalize the critical adsorption conditions. The same approach was implemented in our recent study for the "inverse" problem of polyelectrolyte adsorption onto a planar surface, and on the outer surface of rod-like and spherical obstacles. For the three adsorption geometries investigated, the theory yields simple scaling relations for the minimal surface charge density that triggers the chain adsorption, as a function of the Debye screening length and surface curvature. The encapsulation of polyelectrolytes is governed by interplay of the electrostatic attraction energy toward the adsorbing surface and entropic repulsion of the chain squeezed into a thin slit or small cavities. Under the conditions of surface-mediated confinement, substantially larger polymer linear charge densities are required to adsorb a polyelectrolyte inside a charged spherical cavity, relative to a cylindrical pore and to a planar slit (at the same interfacial surface charge density). Possible biological implications are discussed briefly in the end.
Collapse
Affiliation(s)
- A G Cherstvy
- Institute of Complex Systems, ICS-2, Forschungszentrum Jülich, 52425 Jülich, Germany.
| |
Collapse
|
10
|
Moloney NR, Ozogány K, Rácz Z. Order statistics of 1/fα signals. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:061101. [PMID: 22304034 DOI: 10.1103/physreve.84.061101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2011] [Indexed: 05/31/2023]
Abstract
Order statistics of periodic, Gaussian noise with 1/f(α) power spectrum is investigated. Using simulations and phenomenological arguments, we find three scaling regimes for the average gap d(k) = (x(k) -x(k) + 1) between the kth and (k+1)st largest values of the signal. The result d(k) k(-1), known for independent, identically distributed variables, remains valid for 0 ≤ α < 1. Nontrivial, α-dependent scaling exponents, d(k) k((α-3)/2), emerge for 1 < α < 5, and, finally, α-independent scaling, d(k) ~ k, is obtained for α > 5. The spectra of average ordered values ε(k) =(x(1) - x(k))~ k(β) is also examined. The exponent β is derived from the gap scaling as well as by relating ε(k) to the density of near-extreme states. Known results for the density of near-extreme states combined with scaling suggest that β(α = 2) = 1/2, β(4) = 3/2, and β(∞) = 2 are exact values. We also show that parallels can be drawn between ε(k) and the quantum mechanical spectra of a particle in power-law potentials.
Collapse
Affiliation(s)
- N R Moloney
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str 38, D-01187 Dresden, Germany.
| | | | | |
Collapse
|
11
|
Rambeau J, Schehr G. Distribution of the time at which N vicious walkers reach their maximal height. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:061146. [PMID: 21797341 DOI: 10.1103/physreve.83.061146] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2011] [Indexed: 05/31/2023]
Abstract
We study the extreme statistics of N nonintersecting Brownian motions (vicious walkers) over a unit time interval in one dimension. Using path-integral techniques we compute exactly the joint distribution of the maximum M and of the time τ(M) at which this maximum is reached. We focus in particular on nonintersecting Brownian bridges ("watermelons without wall") and nonintersecting Brownian excursions ("watermelons with a wall"). We discuss in detail the relationships between such vicious walkers models in watermelon configurations and stochastic growth models in curved geometry on the one hand and the directed polymer in a disordered medium (DPRM) with one free end point on the other hand. We also check our results using numerical simulations of Dyson's Brownian motion and confront them with numerical simulations of the polynuclear growth model (PNG) and of a model of DPRM on a discrete lattice. Some of the results presented here were announced in a recent letter [J. Rambeau and G. Schehr, Europhys. Lett. 91, 60006 (2010)].
Collapse
Affiliation(s)
- Joachim Rambeau
- Laboratoire de Physique Théorique d'Orsay, Université Paris Sud 11 and CNRS, Orsay, France.
| | | |
Collapse
|
12
|
Palatella L, Pennetta C. Distribution of first-return times in correlated stationary signals. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:041102. [PMID: 21599110 DOI: 10.1103/physreve.83.041102] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2011] [Indexed: 05/30/2023]
Abstract
We present an analytical expression for the first return time (FRT) probability density function of a stationary correlated signal. Precisely, we start by considering a stationary discrete-time Ornstein-Uhlenbeck (OU) process with exponential decaying correlation function. The first return time distribution for this process is derived by adopting a well-known formalism typically used in the study of the FRT statistics for nonstationary diffusive processes. Then, by a subordination approach, we treat the case of a stationary process with power-law tail correlation function and diverging correlation time. We numerically test our findings, obtaining in both cases a good agreement with the analytical results. We notice that neither in the standard OU nor in the subordinated case a simple form of waiting time statistics, like stretched-exponential or similar, can be obtained while it is apparent that long time transient may shadow the final asymptotic behavior.
Collapse
Affiliation(s)
- Luigi Palatella
- CNISM UdR of Lecce and Dipartimento di Fisica, Università del Salento, Via Arnesano, I-73100 Lecce, Italy
| | | |
Collapse
|
13
|
Moloney NR, Davidsen J. Extreme value statistics and return intervals in long-range correlated uniform deviates. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:041131. [PMID: 19518197 DOI: 10.1103/physreve.79.041131] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2008] [Indexed: 05/27/2023]
Abstract
We study extremal statistics and return intervals in stationary long-range correlated sequences for which the underlying probability density function is bounded and uniform. The extremal statistics we consider (e.g., maximum relative to minimum) are such that the reference point from which the maximum is measured is itself a random quantity. We analytically calculate the limiting distributions for independent and identically distributed random variables, and use these as a reference point for correlated cases. The distributions are different from that of the maximum itself (i.e., a Weibull distribution), reflecting the fact that the distribution of the reference point either dominates over or convolves with the distribution of the maximum. The functional form of the limiting distributions is unaffected by correlations, although the convergence is slower. We show that our findings can be directly generalized to a wide class of stochastic processes. We also analyze return interval distributions, and compare them to recent conjectures of their functional form.
Collapse
Affiliation(s)
- N R Moloney
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada AB T2N 1N4.
| | | |
Collapse
|
14
|
Schehr G, Majumdar SN, Comtet A, Randon-Furling J. Exact distribution of the maximal height of p vicious walkers. PHYSICAL REVIEW LETTERS 2008; 101:150601. [PMID: 18999584 DOI: 10.1103/physrevlett.101.150601] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2008] [Indexed: 05/27/2023]
Abstract
Using path-integral techniques, we compute exactly the distribution of the maximal height Hp of p nonintersecting Brownian walkers over a unit time interval in one dimension, both for excursions p watermelons with a wall, and bridges p watermelons without a wall, for all integer p>or=1. For large p, we show that <Hp> approximately square root 2p (excursions) whereas <Hp> approximately square root p (bridges). Our exact results prove that previous numerical experiments only measured the preasymptotic behaviors and not the correct asymptotic ones. In addition, our method establishes a physical connection between vicious walkers and random matrix theory.
Collapse
Affiliation(s)
- Grégory Schehr
- Laboratoire de Physique Théorique (UMR du CNRS 8627), Université de Paris-Sud, 91405 Orsay Cedex, France
| | | | | | | |
Collapse
|
15
|
Györgyi G, Moloney NR, Ozogány K, Rácz Z. Finite-size scaling in extreme statistics. PHYSICAL REVIEW LETTERS 2008; 100:210601. [PMID: 18518593 DOI: 10.1103/physrevlett.100.210601] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2007] [Indexed: 05/26/2023]
Abstract
We study the deviations from the limit distributions in extreme value statistics arising due to the finite size (FS) of data sets. A renormalization method is introduced for the case of independent, identically distributed (iid) variables, showing that the iid universality classes are subdivided according to the exponent of the FS convergence, which determines the leading order FS shape correction function as well. It is found that, for the correlated systems of subcritical percolation and 1/f;(alpha) stationary (alpha<1) noise, the iid shape correction compares favorably to simulations. Furthermore, for the strongly correlated regime (alpha>1) of 1/f;(alpha) noise, the shape correction is obtained in terms of the limit distribution itself.
Collapse
Affiliation(s)
- G Györgyi
- Institute for Theoretical Physics-HAS, Eötvös University, Pázmány sétány 1/a, 1117 Budapest, Hungary
| | | | | | | |
Collapse
|