BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Satya N. Majumdar

Satya N. Majumdar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechmath-phmath.MPmath.PRcond-mat.softcond-mat.quant-gasPopulations and EvolutionComputationcond-matcond-mat.dis-nnDigital Librarieshep-thMachine Learningmath.COphysics.soc-phquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
53works
0followers
16topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2026arXiv

The distribution of the maximum of independent resetting Brownian motions

The probability distribution of the maximum $M_t$ of a single resetting Brownian motion (RBM) of duration $t$ and resetting rate $r$, properly centred and scaled, is known to converge to the standard Gumbel distribution of the classical extreme value theory. This Gumbel law describes the typical fluctuations of $M_t$ around its average $\sim \ln (r t)$ for large $t$ on a scale of $O(1)$. Here we compute the large-deviation tails of this distribution when $M_t = O(t)$ and show that the large-deviation function has a singularity where the second derivative is discontinuous, signalling a dynamical phase transition. Then we consider a collection of independent RBMs with initial (and resetting) positions uniformly distributed with a density $ρ$ over the negative half-line. We show that the fluctuations in the initial positions of the particles modify the distribution of $M_t$. The average over the initial conditions can be performed in two different ways, in analogy with disordered systems: (i) the annealed case where one averages over all possible initial conditions and (ii) the quenched case where one considers only the contributions coming from typical initial configurations. We show that in the annealed case, the limiting distribution of the maximum is characterized by a new scaling function, different from the Gumbel law but the large-deviation function remains the same as in the single particle case. In contrast, for the quenched case, the limiting (typical) distribution remains Gumbel but the large-deviation behaviors are new and nontrivial. Our analytical results, both for the typical as well as for the large-deviation regime of $M_t$, are verified numerically with extremely high precision, down to $10^{-250}$ for the probability density of $M_t$.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Exact extreme, order and sum statistics in a class of strongly correlated system

Even though strongly correlated systems are abundant, only a few exceptional cases admit analytical solutions. In this paper we present a large class of solvable systems with strong correlations.. We consider a set of $N$ independent and identically distributed (i.i.d) random variables $\{X_1,\, X_2,\ldots, X_N\}$ whose common distribution has a parameter $Y$ (or a set of parameters) which itself is random with its own distribution. For a fixed value of this parameter $Y$, the $X_i$ variables are independent and we call them conditionally independent and identically distributed (c.i.i.d). However, once integrated over the distribution of the parameter $Y$, the $X_i$ variables get strongly correlated, yet retaining a solvable structure for various observables, such as for the sum and the extremes of $X_i$'s. This provides a simple procedure to generate a class of solvable strongly correlated systems. We illustrate how this procedure works via three physical examples where $N$ particles on a line perform independent (i) Brownian motions, (ii) ballistic motions with random initial velocities, and (iii) Lévy flights, but they get strongly correlated via {\it simultaneous resetting} to the origin. Our results are verified in numerical simulations. This procedure can be used to generate an endless variety of solvable strongly correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Metropolis Monte Carlo sampling: convergence, localization transition and optimality

Among random sampling methods, Markov Chain Monte Carlo algorithms are foremost. Using a combination of analytical and numerical approaches, we study their convergence properties towards the steady state, within a random walk Metropolis scheme. Analysing the relaxation properties of some model algorithms sufficiently simple to enable analytic progress, we show that the deviations from the target steady-state distribution can feature a localization transition as a function of the characteristic length of the attempted jumps defining the random walk. While the iteration of the Monte Carlo algorithm converges to equilibrium for all choices of jump parameters, the localization transition changes drastically the asymptotic shape of the difference between the probability distribution reached after a finite number of steps of the algorithm and the target equilibrium distribution. We argue that the relaxation before and after the localisation transition is respectively limited by diffusion and rejection rates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An exactly solvable predator prey model with resetting

We study a simple model of a diffusing particle (the prey) that on encounter with one of a swarm of diffusing predators can either perish or be reset to its original position at the origin. We show that the survival probability of the prey up to time $t$ decays algebraically as $\sim t^{-θ(p, γ)}$ where the exponent $θ$ depends continuously on two parameters of the model, with $p$ denoting the probability that a prey survives upon encounter with a predator and $γ= D_A/(D_A+D_B)$ where $D_A$ and $D_B$ are the diffusion constants of the prey and the predator respectively. We also compute exactly the probability distribution $P(N|t_c)$ of the total number of encounters till the capture time $t_c$ and show that it exhibits an anomalous large deviation form $P(N|t_c)\sim t_c^{- Φ\left(\frac{N}{\ln t_c}=z\right)}$ for large $t_c$. The rate function $Φ(z)$ is computed explicitly. Numerical simulations are in excellent agreement with our analytical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Condensation transition in large deviations of self-similar Gaussian processes with stochastic resetting

We study the fluctuations of the area $A(t)= \int_0^t x(τ)\, dτ$ under a self-similar Gaussian process (SGP) $x(τ)$ with Hurst exponent $H>0$ (e.g., standard or fractional Brownian motion, or the random acceleration process) that stochastically resets to the origin at rate $r$. Typical fluctuations of $A(t)$ scale as $\sim \sqrt{t}$ for large $t$ and on this scale the distribution is Gaussian, as one would expect from the central limit theorem. Here our main focus is on atypically large fluctuations of $A(t)$. In the long-time limit $t\to\infty$, we find that the full distribution of the area takes the form $P_{r}\left(A|t\right)\sim\exp\left[-t^αΦ\left(A/t^β\right)\right]$ with anomalous exponents $α=1/(2H+2)$ and $β= (2H+3)/(4H+4)$ in the regime of moderately large fluctuations, and a different anomalous scaling form $P_{r}\left(A|t\right)\sim\exp\left[-tΨ\left(A/t^{\left(2H+3\right)/2}\right)\right]$ in the regime of very large fluctuations. The associated rate functions $Φ(y)$ and $Ψ(w)$ depend on $H$ and are found exactly. Remarkably, $Φ(y)$ has a singularity that we interpret as a first-order dynamical condensation transition, while $Ψ(w)$ exhibits a second-order dynamical phase transition above which the number of resetting events ceases to be extensive. The parabolic behavior of $Φ(y)$ around the origin $y=0$ correctly describes the typical, Gaussian fluctuations of $A(t)$. Despite these anomalous scalings, we find that all of the cumulants of the distribution $P_{r}\left(A|t\right)$ grow linearly in time, $\langle A^n\rangle_c\approx c_n \, t$, in the long-time limit. For the case of reset Brownian motion (corresponding to $H=1/2$), we develop a recursive scheme to calculate the coefficients $c_n$ exactly and use it to calculate the first 6 nonvanishing cumulants.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Counting statistics for non-interacting fermions in a rotating trap

We study the ground state of $N \gg 1$ noninteracting fermions in a two-dimensional harmonic trap rotating at angular frequency $Ω>0$. The support of the density of the Fermi gas is a disk of radius $R_e$. We calculate the variance of the number of fermions ${\cal N}_R$ inside a disk of radius $R$ centered at the origin for $R$ in the bulk of the Fermi gas. We find rich and interesting behaviours in two different scaling regimes: (i) $Ω/ ω<1 $ and (ii) $1 - Ω/ ω= O(1/N)$, where $ω$ is the angular frequency of the oscillator. In the first regime (i) we find that ${\rm Var}\,{\cal N}_{R}\simeq\left(A\log N+B\right)\sqrt{N}$ and we calculate $A$ and $B$ as functions of $R/R_e$, $Ω$ and $ω$. We also predict the higher cumulants of ${\cal N}_{R}$ and the bipartite entanglement entropy of the disk with the rest of the system. In the second regime (ii), the mean fermion density exhibits a staircase form, with discrete plateaus corresponding to filling $k$ successive Landau levels, as found in previous studies. Here, we show that ${\rm Var}\,{\cal N}_{R}$ is a discontinuous piecewise linear function of $\sim (R/R_e) \sqrt{N}$ within each plateau, with coefficients that we calculate exactly, and with steps whose precise shape we obtain for any $k$. We argue that a similar piecewise linear behavior extends to all the cumulants of ${\cal N}_{R}$ and to the entanglement entropy. We show that these results match smoothly at large $k$ with the above results for $Ω/ω=O(1)$. These findings are nicely confirmed by numerical simulations. Finally, we uncover a universal behavior of ${\rm Var}\,{\cal N}_{R}$ near the fermionic edge. We extend our results to a three-dimensional geometry, where an additional confining potential is applied in the $z$ direction.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Freezing transitions of Brownian particles in confining potentials

We study the mean first passage time (MFPT) to an absorbing target of a one-dimensional Brownian particle subject to an external potential $v(x)$ in a finite domain. We focus on the cases in which the external potential is confining, of the form $v(x)=k|x-x_0|^n/n$, and where the particle&#39;s initial position coincides with $x_0$. We first consider a particle between an absorbing target at $x=0$ and a reflective wall at $x=c$. At fixed $x_0$, we show that when the target distance $c$ exceeds a critical value, there exists a nonzero optimal stiffness $k_{\rm opt}$ that minimizes the MFPT to the target. However, when $c$ lies below the critical value, the optimal stiffness $k_{\rm opt}$ vanishes. Hence, for any value of $n$, the optimal potential stiffness undergoes a continuous &#34;freezing&#34; transition as the domain size is varied. On the other hand, when the reflective wall is replaced by a second absorbing target, the freezing transition in $k_{\rm opt}$ becomes discontinuous. The phase diagram in the $(x_0,n)$-plane then exhibits three dynamical phases and metastability, with a &#34;triple&#34; point at $(x_0/c\simeq 0.17185$, $n\simeq 0.39539)$. For harmonic or higher order potentials $(n\ge 2)$, the MFPT always increases with $k$ at small $k$, for any $x_0$ or domain size. These results are contrasted with problems of diffusion under optimal resetting in bounded domains.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gap probability and full counting statistics in the one dimensional one-component plasma

We consider the $1d$ one-component plasma (OCP) in thermal equilibrium, consisting of $N$ equally charged particles on a line, with pairwise Coulomb repulsion and confined by an external harmonic potential. We study two observables: (i) the distribution of the gap between two consecutive particles in the bulk and (ii) the distribution of the number of particles $N_I$ in a fixed interval $I=[-L,+L]$ inside the bulk, the so-called full-counting-statistics (FCS). For both observables, we compute, for large $N$, the distribution of the typical as well as atypical large fluctuations. We show that the distribution of the typical fluctuations of the gap are described by the scaling form ${\cal P}_{\rm gap, bulk}(g,N) \sim N H_α(g\,N)$, where $α$ is the interaction coupling and the scaling function $H_α(z)$ is computed explicitly. It has a faster than Gaussian tail for large $z$: $H_α(z) \sim e^{-z^3/(96 α)}$ as $z \to \infty$. Similarly, for the FCS, we show that the distribution of the typical fluctuations of $N_I$ is described by the scaling form ${\cal P}_{\rm FCS}(N_I,N) \sim 2α\, U_α[2 α(N_I - \bar{N}_I)]$, where $\bar{N}_I = L\,N/(2 α)$ is the average value of $N_I$ and the scaling function $U_α(z)$ is obtained explicitly. For both observables, we show that the probability of large fluctuations are described by large deviations forms with respective rate functions that we compute explicitly. Our numerical Monte-Carlo simulations are in good agreement with our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Integer partitions and exclusion statistics: Limit shapes and the largest part of Young diagrams

We compute the limit shapes of the Young diagrams of the minimal difference $p$ partitions and provide a simple physical interpretation for the limit shapes. We also calculate the asymptotic distribution of the largest part of the Young diagram and show that the scaled distribution has a Gumbel form for all $p$. This Gumbel statistics for the largest part remains unchanged even for general partitions of the form $E=\sum_i n_i i^{1/ν}$ with $ν>0$ where $n_i$ is the number of times the part $i$ appears.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Number of distinct sites visited by a resetting random walker

We investigate the number $V_p(n)$ of distinct sites visited by an $n$-step resetting random walker on a $d$-dimensional hypercubic lattice with resetting probability $p$. In the case $p=0$, we recover the well-known result that the average number of distinct sites grows for large $n$ as $\langle V_0(n)\rangle\sim n^{d/2}$ for $d<2$ and as $\langle V_0(n)\rangle\sim n$ for $d>2$. For $p>0$, we show that $\langle V_p(n)\rangle$ grows extremely slowly as $\sim \left[\log(n)\right]^d$. We observe that the recurrence-transience transition at $d=2$ for standard random walks (without resetting) disappears in the presence of resetting. In the limit $p\to 0$, we compute the exact crossover scaling function between the two regimes. In the one-dimensional case, we derive analytically the full distribution of $V_p(n)$ in the limit of large $n$. Moreover, for a one-dimensional random walker, we introduce a new observable, which we call imbalance, that measures how much the visited region is symmetric around the starting position. We analytically compute the full distribution of the imbalance both for $p=0$ and for $p>0$. Our theoretical results are verified by extensive numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal Resetting Brownian Bridges

We introduce a resetting Brownian bridge as a simple model to study search processes where the total search time $t_f$ is finite and the searcher returns to its starting point at $t_f$. This is simply a Brownian motion with a Poissonian resetting rate $r$ to the origin which is constrained to start and end at the origin at time $t_f$. We first provide a rejection-free algorithm to generate such resetting bridges in all dimensions by deriving an effective Langevin equation with an explicit space-time dependent drift $\tilde μ({\bf x},t)$ and resetting rate $\tilde r({\bf x}, t)$. We also study the efficiency of the search process in one-dimension by computing exactly various observables such as the mean-square displacement, the hitting probability of a fixed target and the expected maximum. Surprisingly, we find that there exists an optimal resetting rate $r^*$ that maximizes the search efficiency, even in the presence of a bridge constraint. We show however that the physical mechanism responsible for this optimal resetting rate for bridges is entirely different from resetting Brownian motions without the bridge constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reducing mean first passage times with intermittent confining potentials: a realization of resetting processes

During a random search, resetting the searcher&#39;s position from time to time to the starting point often reduces the mean completion time of the process. Although many different resetting models have been studied over the past ten years, only a few can be physically implemented. Here we study theoretically a protocol that can be realised experimentally and which exhibits unusual optimization properties. A Brownian particle is subject to an arbitrary confining potential $v(x)$ which is switched on and off intermittently at fixed rates. Motion is constrained between an absorbing wall located at the origin and a reflective wall. When the walls are sufficiently far apart, the interplay between free diffusion during the &#34;off&#34; phases and attraction toward the potential minimum during the &#34;on&#34; phases gives rise to rich behaviours, not observed in ideal resetting models. For potentials of the form $v(x)=k|x-x_0|^n/n$, with $n>0$, the switch-on and switch-off rates that minimise the mean first passage time (MFPT) to the origin undergo a continuous phase transition as the potential stiffness $k$ is varied. When $k$ is above a critical value $k_c$, potential intermittency enhances target encounter: the minimal MFPT is lower than the Kramer&#39;s time and is attained for a non-vanishing pair of switching rates. We focus on the harmonic case $n=2$, extending previous results for the piecewise linear potential ($n=1$) in unbounded domains. We also study the non-equilibrium stationary states emerging in this process.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Run-and-tumble particle in one-dimensional confining potential: Steady state, relaxation and first passage properties

We study the dynamics of a one-dimensional run and tumble particle subjected to confining potentials of the type $V(x) = α\, |x|^p$, with $p>0$. The noise that drives the particle dynamics is telegraphic and alternates between $\pm 1$ values. We show that the stationary probability density $P(x)$ has a rich behavior in the $(p, α)$-plane. For $p>1$, the distribution has a finite support in $[x_-,x_+]$ and there is a critical line $α_c(p)$ that separates an active-like phase for $α> α_c(p)$ where $P(x)$ diverges at $x_\pm$, from a passive-like phase for $α< α_c(p)$ where $P(x)$ vanishes at $x_\pm$. For $p<1$, the stationary density $P(x)$ collapses to a delta function at the origin, $P(x) = δ(x)$. In the marginal case $p=1$, we show that, for $α< α_c$, the stationary density $P(x)$ is a symmetric exponential, while for $α> α_c$, it again is a delta function $P(x) = δ(x)$. For the special cases $p=2$ and $p=1$, we obtain exactly the full time-dependent distribution $P(x,t)$, that allows us to study how the system relaxes to its stationary state. In addition, in these two cases, we also study analytically the full distribution of the first-passage time to the origin. Numerical simulations are in complete agreement with our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Edge fluctuations and third-order phase transition in harmonically confined long-range systems

We study the distribution of the position of the rightmost particle $x_{\max}$ in a $N$-particle Riesz gas in one dimension confined in a harmonic trap. The particles interact via long-range repulsive potential, of the form $r^{-k}$ with $-2<k<\infty$ where $r$ is the inter-particle distance. In equilibrium at temperature $O(1)$, the gas settles on a finite length scale $L_N$ that depends on $N$ and $k$. We numerically observe that the typical fluctuation of $y_{\max} = x_{\max}/L_N$ around its mean is of $O(N^{-η_k})$. Over this length scale, the distribution of the typical fluctuations has a $N$ independent scaling form. We show that the exponent $η_k$ obtained from the Hessian theory predicts the scale of typical fluctuations remarkably well. The distribution of atypical fluctuations to the left and right of the mean $\langle y_{\max} \rangle$ are governed by the left and right large deviation functions, respectively. We compute these large deviation functions explicitly $\forall k>-2$. We also find that these large deviation functions describe a pulled to pushed type phase transition as observed in Dyson&#39;s log-gas ($k\to 0$) and $1d$ one component plasma ($k=-1$). Remarkably, we find that the phase transition remains $3^{\rm rd}$ order for the entire regime. Our results demonstrate the striking universality of the $3^{\rm rd}$ order transition even in models that fall outside the paradigm of Coulomb systems and the random matrix theory. We numerically verify our analytical expressions of the large deviation functions via Monte Carlo simulation using an importance sampling algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generating stochastic trajectories with global dynamical constraints

We propose a method to exactly generate Brownian paths $x_c(t)$ that are constrained to return to the origin at some future time $t_f$, with a given fixed area $A_f = \int_0^{t_f}dt\, x_c(t)$ under their trajectory. We derive an exact effective Langevin equation with an effective force that accounts for the constraint. In addition, we develop the corresponding approach for discrete-time random walks, with arbitrary jump distributions including Lévy flights, for which we obtain an effective jump distribution that encodes the constraint. Finally, we generalise our method to other types of dynamical constraints such as a fixed occupation time on the positive axis $T_f=\int_0^{t_f}dt\, Θ\left[x_c(t)\right]$ or a fixed generalised quadratic area $\mathcal{A}_f=\int_0^{t_f}dt \,x_c^2(t)$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Localization transition in the Discrete Non-Linear Schrödinger Equation: ensembles inequivalence and negative temperatures

We present a detailed account of a first-order localization transition in the Discrete Nonlinear Schrödinger Equation, where the localized phase is associated to the high energy region in parameter space. We show that, due to ensemble inequivalence, this phase is thermodynamically stable only in the microcanonical ensemble. In particular, we obtain an explicit expression of the microcanonical entropy close to the transition line, located at infinite temperature. This task is accomplished making use of large-deviation techniques, that allow us to compute, in the limit of large system size, also the subleading corrections to the microcanonical entropy. These subleading terms are crucial ingredients to account for the first-order mechanism of the transition, to compute its order parameter and to predict the existence of negative temperatures in the localized phase. All of these features can be viewed as signatures of a thermodynamic phase where the translational symmetry is broken spontaneously due to a condensation mechanism yielding energy fluctuations far away from equipartition: actually they prefer to participate in the formation of nonlinear localized excitations (breathers), typically containing a macroscopic fraction of the total energy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mean area of the convex hull of a run and tumble particle in two dimensions

We investigate the statistics of the convex hull for a single run-and-tumble particle in two dimensions. Run-and-tumble particle, also known as persistent random walker, has gained significant interest in the recent years due to its biological application in modelling the motion of bacteria. We consider two different statistical ensembles depending on whether (i) the total number of tumbles $n$ or (ii) the total observation time $t$ is kept fixed. Benchmarking the results on perimeter, we study the statistical properties of the area of the convex hull for RTP. Exploiting the connections to extreme value statistics, we obtain exact analytical expressions for the mean area for both ensembles. For fixed-$t$ ensemble, we show that the mean possesses a scaling form in $γt$ (with $γ$ being the tumbling rate) and the corresponding scaling function is exactly computed. Interestingly, we find that it exhibits crossover from $\sim t^3$ scaling at small times $\left( t \ll γ^{-1} \right)$ to $\sim t$ scaling at large times $\left( t \gg γ^{-1} \right)$. On the other hand, for fixed-$n$ ensemble, the mean expectedly grows linearly with $n$ for $n \gg 1$. All our analytical findings are supported with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimal mean first-passage time of a Brownian searcher with resetting in one and two dimensions: Experiments, theory and numerical tests

We study experimentally, numerically and theoretically the optimal mean time needed by a Brownian particle, freely diffusing either in one or two dimensions, to reach, within a tolerance radius $R_{\text tol}$, a target at a distance $L$ from an initial position in the presence of resetting. The reset position is Gaussian distributed with width $σ$. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different target distances (values of the ratios $b=L/σ$) and target size ($a=R_\text{tol}/L$). We find an interesting phase transition at a critical value of $b$, both in one and two dimensions. The details of the calculations as well as experimental setup and limitations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Record statistics for random walks and Lévy flights with resetting

We compute exactly the mean number of records $\langle R_N \rangle$ for a time-series of size $N$ whose entries represent the positions of a discrete time random walker on the line. At each time step, the walker jumps by a length $η$ drawn independently from a symmetric and continuous distribution $f(η)$ with probability $1-r$ (with $0\leq r < 1$) and with the complementary probability $r$ it resets to its starting point $x=0$. This is an exactly solvable example of a weakly correlated time-series that interpolates between a strongly correlated random walk series (for $r=0$) and an uncorrelated time-series (for $(1-r) \ll 1$). Remarkably, we found that for every fixed $r \in [0,1[$ and any $N$, the mean number of records $\langle R_N \rangle$ is completely universal, i.e., independent of the jump distribution $f(η)$. In particular, for large $N$, we show that $\langle R_N \rangle$ grows very slowly with increasing $N$ as $\langle R_N \rangle \approx (1/\sqrt{r})\, \ln N$ for $0<r <1$. We also computed the exact universal crossover scaling functions for $\langle R_N \rangle$ in the two limits $r \to 0$ and $r \to 1$. Our analytical predictions are in excellent agreement with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stability of large complex systems with heterogeneous relaxation dynamics

We study the probability of stability of a large complex system of size $N$ within the framework of a generalized May model, which assumes a linear dynamics of each population size $n_i$ (with respect to its equilibrium value): $ \frac{\mathrm{d}\, n_i}{\mathrm{d}t} = - a_i n_i - \sqrt{T} \sum_{j} J_{ij} n_j $. The $a_i>0$&#39;s are the intrinsic decay rates, $J_{ij}$ is a real symmetric $(N\times N)$ Gaussian random matrix and $\sqrt{T}$ measures the strength of pairwise interaction between different species. Unlike in May&#39;s original homogeneous model, each species has now an intrinsic damping $a_i$ that may differ from one another. As the interaction strength $T$ increases, the system undergoes a phase transition from a stable phase to an unstable phase at a critical value $T=T_c$. We reinterpret the probability of stability in terms of the hitting time of the level $b=0$ of an associated Dyson Brownian Motion (DBM), starting at the initial position $a_i$ and evolving in `time&#39; $T$. In the large $N \to \infty$ limit, using this DBM picture, we are able to completely characterize $T_c$ for arbitrary density $μ(a)$ of the $a_i$&#39;s. For a specific flat configuration $a_i = 1 + σ\frac{i-1}{N}$, we obtain an explicit parametric solution for the limiting (as $N\to \infty$) spectral density for arbitrary $T$ and $σ$. For finite but large $N$, we also compute the large deviation properties of the probability of stability on the stable side $T < T_c$ using a Coulomb gas representation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Statistics of first-passage Brownian functionals

We study the distribution of first-passage functionals ${\cal A}= \int_0^{t_f} x^n(t)\, dt$, where $x(t)$ is a Brownian motion (with or without drift) with diffusion constant $D$, starting at $x_0>0$, and $t_f$ is the first-passage time to the origin. In the driftless case, we compute exactly, for all $n>-2$, the probability density $P_n(A|x_0)=\text{Prob}.(\mathcal{A}=A)$. This probability density has an essential singular tail as $A\to 0$ and a power-law tail $\sim A^{-(n+3)/(n+2)}$ as $A\to \infty$. The former is reproduced by the optimal fluctuation method (OFM), which also predicts the optimal paths of the conditioned process for small $A$. For the case with a drift toward the origin, where no exact solution is known for general $n>-1$, the OFM predicts the distribution tails. For $A\to 0$ it predicts the same essential singular tail as in the driftless case. For $A\to \infty$ it predicts a stretched exponential tail $-\ln P_n(A|x_0)\sim A^{1/(n+1)}$ for all $n>0$. In the limit of large Péclet number $\text{Pe}= μx_0/(2D)\gg 1$, where $μ$ is the drift velocity, the OFM predicts a large-deviation scaling for all $A$: $-\ln P_n(A|x_0)\simeq\text{Pe}\, Φ_n\left(z= A/\bar{A}\right)$, where $\bar{A}=x_0^{n+1}/{μ(n+1)}$ is the mean value of $\mathcal{A}$. We compute the rate function $Φ_n(z)$ analytically for all $n>-1$. For $n>0$ $Φ_n(z)$ is analytic for all $z$, but for $-1<n<0$ it is non-analytic at $z=1$, implying a dynamical phase transition. The order of this transition is $2$ for $-1/2<n<0$, while for $-1<n<-1/2$ the order of transition changes continuously with $n$. Finally, we apply the OFM to the case of $μ<0$ (drift away from the origin). We show that, when the process is conditioned on reaching the origin, the distribution of $\mathcal{A}$ coincides with the distribution of $\mathcal{A}$ for $μ>0$ with the same $|μ|$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Statistics of the maximum and the convex hull of a Brownian motion in confined geometries

We consider a Brownian particle with diffusion coefficient $D$ in a $d$-dimensional ball of radius $R$ with reflecting boundaries. We study the maximum $M_x(t)$ of the trajectory of the particle along the $x$-direction at time $t$. In the long time limit, the maximum converges to the radius of the ball $M_x(t) \to R$ for $t\to \infty$. We investigate how this limit is approached and obtain an exact analytical expression for the distribution of the fluctuations $Δ(t) = [R-M_x(t)]/R$ in the limit of large $t$ in all dimensions. We find that the distribution of $Δ(t)$ exhibits a rich variety of behaviors depending on the dimension $d$. These results are obtained by establishing a connection between this problem and the narrow escape time problem. We apply our results in $d=2$ to study the convex hull of the trajectory of the particle in a disk of radius $R$ with reflecting boundaries. We find that the mean perimeter $\langle L(t)\rangle$ of the convex hull exhibits a slow convergence towards the perimeter of the circle $2πR$ with a stretched exponential decay $2πR-\langle L(t)\rangle \propto \sqrt{R}(Dt)^{1/4} \,e^{-2\sqrt{2Dt}/R}$. Finally, we generalise our results to other confining geometries, such as the ellipse with reflecting boundaries. Our results are corroborated by thorough numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Active Brownian Motion in Two Dimensions

We study the dynamics of a single active Brownian particle (ABP) in two spatial dimensions. The ABP has an intrinsic time scale $D_R^{-1}$ set by the rotational diffusion constant $D_R$. We show that, at short-times $t \ll D_R^{-1}$, the presence of `activness&#39; results in a strongly anisotropic and non-diffusive dynamics in the $(xy)$ plane. We compute exactly the marginal distributions of the $x$ and $y$ position coordinates along with the radial distribution, which are all shown to be non-Brownian. In addition, we show that, at early times, the ABP has anomalous first-passage properties, characterized by non-Brownian exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anderson-like localization transition of random walks with resetting

We study several lattice random walk models with stochastic resetting to previously visited sites which exhibit a phase transition between an anomalous diffusive regime and a localization regime where diffusion is suppressed. The localized phase settles above a critical resetting rate, or rate of memory use, and the probability density asymptotically adopts in this regime a non-equilibrium steady state similar to that of the well known problem of diffusion with resetting to the origin. The transition occurs because of the presence of a single impurity site where the resetting rate is lower than on other sites, and around which the walker spontaneously localizes. Near criticality, the localization length diverges with a critical exponent that falls in the same class as the self-consistent theory of Anderson localization of waves in random media. The critical dimensions are also the same in both problems. Our study provides analytically tractable examples of localization transitions in path-dependent, reinforced stochastic processes, which can be also useful for understanding spatial learning by living organisms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Current fluctuations in non-interacting run-and-tumble particles in one-dimension

We present a general framework to study the distribution of the flux through the origin up to time $t$, in a non-interacting one-dimensional system of particles with a step initial condition with a fixed density $ρ$ of particles to the left of the origin. We focus principally on two cases: (i) when the particles undergo diffusive dynamics (passive case) and (ii) run-and-tumble dynamics for each particle (active case). In analogy with disordered systems, we consider the flux distribution both for the annealed and the quenched initial conditions, for the passive and active particles. In the annealed case, we show that, for arbitrary particle dynamics, the flux distribution is a Poissonian with a mean $μ(t)$ that we compute exactly in terms of the Green&#39;s function of the single particle dynamics. For the quenched case, we show that, for the run-and-tumble dynamics, the quenched flux distribution takes an anomalous large deviation form at large times $P_{\rm qu}(Q,t) \sim \exp\left[-ρ\, v_0\, γ\, t^2 ψ_{\rm RTP}\left(\frac{Q}{ρv_0\,t} \right) \right]$, where $γ$ is the rate of tumbling and $v_0$ is the ballistic speed between two successive tumblings. In this paper, we compute the rate function $ψ_{\rm RTP}(q)$ and show that it is nontrivial. Our method also gives access to the probability of the rare event that, at time $t$, there is no particle to the right of the origin. For diffusive and run-and-tumble dynamics, we find that this probability decays with time as a stretched exponential, $\sim \exp(-c\, \sqrt{t})$ where the constant $c$ can be computed exactly. We verify our results for these large deviations by using an importance sampling Monte-Carlo method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distribution of the Time Between Maximum and Minimum of Random Walks

We consider a one-dimensional Brownian motion of fixed duration $T$. Using a path-integral technique, we compute exactly the probability distribution of the difference $τ=t_{\min}-t_{\max}$ between the time $t_{\min}$ of the global minimum and the time $t_{\max}$ of the global maximum. We extend this result to a Brownian bridge, i.e. a periodic Brownian motion of period $T$. In both cases, we compute analytically the first few moments of $τ$, as well as the covariance of $t_{\max}$ and $t_{\min}$, showing that these times are anti-correlated. We demonstrate that the distribution of $τ$ for Brownian motion is valid for discrete-time random walks with $n$ steps and with a finite jump variance, in the limit $n\to \infty$. In the case of Lévy flights, which have a divergent jump variance, we numerically verify that the distribution of $τ$ differs from the Brownian case. For random walks with continuous and symmetric jumps we numerically verify that the probability of the event &#34;$τ= n$&#34; is exactly $1/(2n)$ for any finite $n$, independently of the jump distribution. Our results can be also applied to describe the distance between the maximal and minimal height of $(1+1)$-dimensional stationary-state Kardar-Parisi-Zhang interfaces growing over a substrate of finite size $L$. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 123, 200201 (2019)].

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ising model with stochastic resetting

We study the stationary properties of the Ising model that, while evolving towards its equilibrium state at temperature $T$ according to the Glauber dynamics, is stochastically reset to its fixed initial configuration with magnetisation $m_0$ at a constant rate $r$. Resetting breaks detailed balance and drives the system to a non-equilibrium stationary state where the magnetisation acquires a nontrivial distribution, leading to a rich phase diagram in the $(T,r)$ plane. We establish these results exactly in one-dimension and present scaling arguments supported by numerical simulations in two-dimensions. We show that resetting gives rise to a novel &#34;pseudo-ferro&#34; phase in the $(T,r)$ plane for $r > r^*(T)$ and $T>T_c$ where $r^*(T)$ is a crossover line separating the pseudo-ferro phase from a paramagnetic phase. This pseudo-ferro phase is characterised by a non-zero typical magnetisation and a vanishing gap near $m=0$ of the magnetisation distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kernels for noninteracting fermions via a Green&#39;s function approach with applications to step potentials

The quantum correlations of $N$ noninteracting spinless fermions in their ground state can be expressed in terms of a two-point function called the kernel. Here we develop a general and compact method for computing the kernel in a general trapping potential in terms of the Green&#39;s function for the corresponding single particle Schrödinger equation. For smooth potentials the method allows a simple alternative derivation of the local density approximation for the density and of the sine kernel in the bulk part of the trap in the large $N$ limit. It also recovers the density and the kernel of the so-called {\em Airy gas} at the edge. This method allows to analyse the quantum correlations in the ground state when the potential has a singular part with a fast variation in space. For the square step barrier of height $V_0$, we derive explicit expressions for the density and for the kernel. For large Fermi energy $μ>V_0$ it describes the interpolation between two regions of different densities in a Fermi gas, each described by a different sine kernel. Of particular interest is the {\em critical point} of the square well potential when $μ=V_0$. In this critical case, while there is a macroscopic number of fermions in the lower part of the step potential, there is only a finite $O(1)$ number of fermions on the shoulder, and moreover this number is independent of $μ$. In particular, the density exhibits an algebraic decay $\sim 1/x^2$, where $x$ is the distance from the jump. Furthermore, we show that the critical behaviour around $μ= V_0$ exhibits universality with respect with the shape of the barrier. This is established (i) by an exact solution for a smooth barrier (the Woods-Saxon potential) and (ii) by establishing a general relation between the large distance behavior of the kernel and the scattering amplitudes of the single-particle wave-function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mean perimeter and area of the convex hull of a planar Brownian motion in the presence of resetting

We compute exactly the mean perimeter and the mean area of the convex hull of a $2$-d Brownian motion of duration $t$ and diffusion constant $D$, in the presence of resetting to the origin at a constant rate $r$. We show that for any $t$, the mean perimeter is given by $\langle L(t)\rangle= 2 π\sqrt{\frac{D}{r}}\, f_1(rt)$ and the mean area is given by $\langle A(t) \rangle= 2π\frac{D}{r}\, f_2(rt)$ where the scaling functions $f_1(z)$ and $f_2(z)$ are computed explicitly. For large $t\gg 1/r$, the mean perimeter grows extremely slowly as $\langle L(t)\rangle \propto \ln (rt)$ with time. Likewise, the mean area also grows slowly as $\langle A(t)\rangle \propto \ln^2(rt)$ for $t\gg 1/r$. Our exact results indicate that the convex hull, in the presence of resetting, approaches a circular shape at late times. Numerical simulations are in perfect agreement with our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Noninteracting trapped Fermions in double-well potentials: inverted parabola kernel

We study a system of $N$ noninteracting spinless fermions in a confining, double-well potential in one dimension. When the Fermi energy is close to the value of the potential at its local maximum we show that physical properties, such as the average density and the fermion position correlation functions, display a universal behavior that depends only on the local properties of the potential near its maximum. This behavior describes the merging of two Fermi gases, which are disjoint at sufficiently low Fermi energies. We describe this behavior in terms of a new correlation kernel that we compute analytically and we call it the inverted parabola kernel&#34;. As an application, we calculate the mean and variance of the number of particles in an interval of size $2L$ centered around the position of the local maximum, for sufficiently small $L$. Finally, we discuss the possibility of observing our results in experiments, as well as the extensions to nonzero temperature and to higher space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal mean first-passage time for a Brownian searcher subjected to resetting: experimental and theoretical results

We study experimentally and theoretically the optimal mean time needed by a free diffusing Brownian particle to reach a target at a distance L from an initial position in the presence of resetting. Both the initial position and the resetting position are Gaussian distributed with width $σ$. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different values of the ratio b = L/$σ$ and find an interesting phase transtion at a critical value b = bc. For bc < b < $\infty$, there is a metastable optimum time which disappears for b < bc. The intrinsic diffculties in implementing these protocols in experiments are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Position distribution in a generalised run and tumble process

We study a class of stochastic processes of the type $\frac{d^n x}{dt^n}= v_0\, σ(t)$ where $n>0$ is a positive integer and $σ(t)=\pm 1$ represents an `active&#39; telegraphic noise that flips from one state to the other with a constant rate $γ$. For $n=1$, it reduces to the standard run and tumble process for active particles in one dimension. This process can be analytically continued to any $n>0$ including non-integer values. We compute exactly the mean squared displacement at time $t$ for all $n>0$ and show that at late times while it grows as $\sim t^{2n-1}$ for $n>1/2$, it approaches a constant for $n<1/2$. In the marginal case $n=1/2$, it grows very slowly with time as $\sim \ln t$. Thus the process undergoes a {\em localisation} transition at $n=1/2$. We also show that the position distribution $p_n(x,t)$ remains time-dependent even at late times for $n\ge 1/2$, but approaches a stationary time-independent form for $n<1/2$. The tails of the position distribution at late times exhibit a large deviation form, $p_n(x,t)\sim \exp\left[-γ\, t\, Φ_n\left(\frac{x}{x^*(t)}\right)\right]$, where $x^*(t)= v_0\, t^n/Γ(n+1)$. We compute the rate function $Φ_n(z)$ analytically for all $n>0$ and also numerically using importance sampling methods, finding excellent agreement between them. For three special values $n=1$, $n=2$ and $n=1/2$ we compute the exact cumulant generating function of the position distribution at all times $t$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum dynamics with stochastic reset

We study non-equilibrium dynamics of integrable and non-integrable closed quantum systems whose unitary evolution is interrupted with stochastic resets, characterized by a reset rate $r$, that project the system to its initial state. We show that the steady state density matrix of a non-integrable system, averaged over the reset distribution, retains its off-diagonal elements for any finite $r$. Consequently a generic observable $\hat O$, whose expectation value receives contribution from these off-diagonal elements, never thermalizes under such dynamics for any finite $r$. We demonstrate this phenomenon by exact numerical studies of experimentally realizable models of ultracold bosonic atoms in a tilted optical lattice. For integrable Dirac-like fermionic models driven periodically between such resets, the reset-averaged steady state is found to be described by a family of generalized Gibbs ensembles (GGE s) characterized by $r$. We also study the spread of particle density of a non-interacting one-dimensional fermionic chain, starting from an initial state where all fermions occupy the left half of the sample, while the right half is empty. When driven by resetting dynamics, the density profile approaches at long times to a nonequilibrium stationary profile that we compute exactly. We suggest concrete experiments that can possibly test our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Statistics of the Number of Records for Random Walks and Lévy Flights on a ${1D}$ Lattice

We study the statistics of the number of records $R_n$ for a symmetric, $n$-step, discrete jump process on a $1D$ lattice. At a given step, the walker can jump by arbitrary lattice units drawn from a given symmetric probability distribution. This process includes, as a special case, the standard nearest neighbor lattice random walk. We derive explicitly the generating function of the distribution $P(R_n)$ of the number of records, valid for arbitrary discrete jump distributions. As a byproduct, we provide a relatively simple proof of the generalized Sparre Andersen theorem for the survival probability of a random walk on a line, with discrete or continuous jump distributions. For the discrete jump process, we then derive the asymptotic large $n$ behavior of $P(R_n)$ as well as of the average number of records $E(R_n)$. We show that unlike the case of random walks with symmetric and continuous jump distributions where the record statistics is strongly universal (i.e., independent of the jump distribution for all $n$), the record statistics for lattice walks depends on the jump distribution for any fixed $n$. However, in the large $n$ limit, we show that the distribution of the scaled record number $R_n/E(R_n)$ approaches a universal, half-Gaussian form for any discrete jump process. The dependence on the jump distribution enters only through the scale factor $E(R_n)$, which we also compute in the large $n$ limit for arbitrary jump distributions. We present explicit results for a few examples and provide numerical checks of our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic Resetting and Applications

In this Topical Review we consider stochastic processes under resetting, which have attracted a lot of attention in recent years. We begin with the simple example of a diffusive particle whose position is reset randomly in time with a constant rate $r$, which corresponds to Poissonian resetting, to some fixed point (e.g. its initial position). This simple system already exhibits the main features of interest induced by resetting: (i) the system reaches a nontrivial nonequilibrium stationary state (ii) the mean time for the particle to reach a target is finite and has a minimum, optimal, value as a function of the resetting rate $r$. We then generalise to an arbitrary stochastic process (e.g. Lévy flights or fractional Brownian motion) and non-Poissonian resetting (e.g. power-law waiting time distribution for intervals between resetting events). We go on to discuss multiparticle systems as well as extended systems, such as fluctuating interfaces, under resetting. We also consider resetting with memory which implies resetting the process to some randomly selected previous time. Finally we give an overview of recent developments and applications in the field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time between the maximum and the minimum of a stochastic process

We present an exact solution for the probability density function $P(τ=t_{\min}-t_{\max}|T)$ of the time-difference between the minimum and the maximum of a one-dimensional Brownian motion of duration $T$. We then generalise our results to a Brownian bridge, i.e. a periodic Brownian motion of period $T$. We demonstrate that these results can be directly applied to study the position-difference between the minimal and the maximal height of a fluctuating $(1+1)$-dimensional Kardar-Parisi-Zhang interface on a substrate of size $L$, in its stationary state. We show that the Brownian motion result is universal and, asymptotically, holds for any discrete-time random walk with a finite jump variance. We also compute this distribution numerically for Lévy flights and find that it differs from the Brownian motion result.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Toward the full short-time statistics of an active Brownian particle on the plane

We study the position distribution of a single active Brownian particle (ABP) on the plane. We show that this distribution has a compact support, the boundary of which is an expanding circle. We focus on a short-time regime and employ the optimal fluctuation method (OFM) to study large deviations of the particle position coordinates $x$ and $y$. We determine the optimal paths of the ABP, conditioned on reaching specified values of $x$ and $y$, and the large deviation functions of the marginal distributions of $x$, and of $y$. These marginal distributions match continuously with &#34;near tails&#34; of the $x$ and $y$ distributions of typical fluctuations, studied earlier. We also calculate the large deviation function of the joint $x$ and $y$ distribution $P(x,y,t)$ in a vicinity of a special &#34;zero-noise&#34; point, and show that $\ln P(x,y,t)$ has a nontrivial self-similar structure as a function of $x$, $y$ and $t$. The joint distribution vanishes extremely fast at the expanding circle, exhibiting an essential singularity there. This singularity is inherited by the marginal $x$- and $y$-distributions. We argue that this fingerprint of the short-time dynamics remains there at all times.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal gap statistics for random walks for a class of jump densities

We study the order statistics of a random walk (RW) of $n$ steps whose jumps are distributed according to symmetric Erlang densities $f_p(η)\sim |η|^p \,e^{-|η|}$, parametrized by a non-negative integer $p$. Our main focus is on the statistics of the gaps $d_{k,n}$ between two successive maxima $d_{k,n}=M_{k,n}-M_{k+1,n}$ where $M_{k,n}$ is the $k$-th maximum of the RW between step 1 and step $n$. In the limit of large $n$, we show that the probability density function of the gaps $P_{k,n}(Δ) = \Pr(d_{k,n} = Δ)$ reaches a stationary density $P_{k,n}(Δ) \to p_k(Δ)$. For large $k$, we demonstrate that the typical fluctuations of the gap, for $d_{k,n}= O(1/\sqrt{k})$ (and $n \to \infty$), are described by a non-trivial scaling function that is independent of $k$ and of the jump probability density function $f_p(η)$, thus corroborating our conjecture about the universality of the regime of typical fluctuations (see G. Schehr, S. N. Majumdar, Phys. Rev. Lett. 108, 040601 (2012)). We also investigate the large fluctuations of the gap, for $d_{k,n} = O(1)$ (and $n \to \infty$), and show that these two regimes of typical and large fluctuations of the gaps match smoothly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal survival probability for a $d$-dimensional run-and-tumble particle

We consider an active run-and-tumble particle (RTP) in $d$ dimensions and compute exactly the probability $S(t)$ that the $x$-component of the position of the RTP does not change sign up to time $t$. When the tumblings occur at a constant rate, we show that $S(t)$ is independent of $d$ for any finite time $t$ (and not just for large $t$), as a consequence of the celebrated Sparre Andersen theorem for discrete-time random walks in one dimension. Moreover, we show that this universal result holds for a much wider class of RTP models in which the speed $v$ of the particle after each tumbling is random, drawn from an arbitrary probability distribution. We further demonstrate, as a consequence, the universality of the record statistics in the RTP problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Velocity and diffusion constant of an active particle in a one dimensional force field

We consider a run an tumble particle with two velocity states $\pm v_0$, in an inhomogeneous force field $f(x)$ in one dimension. We obtain exact formulae for its velocity $V_L$ and diffusion constant $D_L$ for arbitrary periodic $f(x)$ of period $L$. They involve the &#34;active potential&#34; which allows to define a global bias. Upon varying parameters, such as an external force $F$, the dynamics undergoes transitions from non-ergodic trapped states, to various moving states, some with non analyticities in the $V_L$ versus $F$ curve. A random landscape in the presence of a bias leads, for large $L$, to anomalous diffusion $x \sim t^μ$, $μ<1$, or to a phase with a finite velocity that we calculate.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A First-Order Dynamical Transition in the displacement distribution of a Driven Run-and-Tumble Particle

We study the probability distribution $P(X_N=X,N)$ of the total displacement $X_N$ of an $N$-step run and tumble particle on a line, in presence of a constant nonzero drive $E$. While the central limit theorem predicts a standard Gaussian form for $P(X,N)$ near its peak, we show that for large positive and negative $X$, the distribution exhibits anomalous large deviation forms. For large positive $X$, the associated rate function is nonanalytic at a critical value of the scaled distance from the peak where its first derivative is discontinuous. This signals a first-order dynamical phase transition from a homogeneous `fluid&#39; phase to a `condensed&#39; phase that is dominated by a single large run. A similar first-order transition occurs for negative large fluctuations as well. Numerical simulations are in excellent agreement with our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exact stationary state of a run-and-tumble particle with three internal states in a harmonic trap

We study the motion of a one-dimensional run-and-tumble particle with three discrete internal states in the presence of a harmonic trap of stiffness $μ.$ The three internal states, corresponding to positive, negative and zero velocities respectively, evolve following a jump process with rate $γ$. We compute the stationary position distribution exactly for arbitrary values of $μ$ and $γ$ which turns out to have a finite support on the real line. We show that the distribution undergoes a shape-transition as $β=γ/μ$ is changed. For $β<1,$ the distribution has a double-concave shape and shows algebraic divergences with an exponent $(β-1)$ both at the origin and at the boundaries. For $β>1,$ the position distribution becomes convex, vanishing at the boundaries and with a single, finite, peak at the origin. We also show that for the special case $β=1,$ the distribution shows a logarithmic divergence near the origin while saturating to a constant value at the boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Mean perimeter and mean area of the convex hull over planar random walks

We investigate the geometric properties of the convex hull over $n$ successive positions of a planar random walk, with a symmetric continuous jump distribution. We derive the large $n$ asymptotic behavior of the mean perimeter. In addition, we compute the mean area for the particular case of isotropic Gaussian jumps. While the leading terms of these asymptotics are universal, the subleading (correction) terms depend on finer details of the jump distribution and describe a &#34;finite size effect&#34; of discrete-time jump processes, allowing one to accurately compute the mean perimeter and the mean area even for small $n$, as verified by Monte Carlo simulations. This is particularly valuable for applications dealing with discrete-time jumps processes and ranging from the statistical analysis of single-particle tracking experiments in microbiology to home range estimations in ecology.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A simple derivation of the Tracy-Widom distribution of the maximal eigenvalue of a Gaussian unitary random matrix

In this paper, we first briefly review some recent results on the distribution of the maximal eigenvalue of a $(N\times N)$ random matrix drawn from Gaussian ensembles. Next we focus on the Gaussian Unitary Ensemble (GUE) and by suitably adapting a method of orthogonal polynomials developed by Gross and Matytsin in the context of Yang-Mills theory in two dimensions, we provide a rather simple derivation of the Tracy-Widom law for GUE. Our derivation is based on the elementary asymptotic scaling analysis of a pair of coupled nonlinear recursion relations. As an added bonus, this method also allows us to compute the precise subleading terms describing the right large deviation tail of the maximal eigenvalue distribution. In the Yang-Mills language, these subleading terms correspond to non-perturbative (in $1/N$ expansion) corrections to the two-dimensional partition function in the so called `weak&#39; coupling regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-intersecting Brownian walkers and Yang-Mills theory on the sphere

We study a system of N non-intersecting Brownian motions on a line segment [0,L] with periodic, absorbing and reflecting boundary conditions. We show that the normalized reunion probabilities of these Brownian motions in the three models can be mapped to the partition function of two-dimensional continuum Yang-Mills theory on a sphere respectively with gauge groups U(N), Sp(2N) and SO(2N). Consequently, we show that in each of these Brownian motion models, as one varies the system size L, a third order phase transition occurs at a critical value L=L_c(N)\sim \sqrt{N} in the large N limit. Close to the critical point, the reunion probability, properly centered and scaled, is identical to the Tracy-Widom distribution describing the probability distribution of the largest eigenvalue of a random matrix. For the periodic case we obtain the Tracy-Widom distribution corresponding to the GUE random matrices, while for the absorbing and reflecting cases we get the Tracy-Widom distribution corresponding to GOE random matrices. In the absorbing case, the reunion probability is also identified as the maximal height of N non-intersecting Brownian excursions (&#34;watermelons&#34; with a wall) whose distribution in the asymptotic scaling limit is then described by GOE Tracy-Widom law. In addition, large deviation formulas for the maximum height are also computed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Area distribution and the average shape of a Lévy bridge

We consider a one dimensional Lévy bridge x_B of length n and index 0 < α< 2, i.e. a Lévy random walk constrained to start and end at the origin after n time steps, x_B(0) = x_B(n)=0. We compute the distribution P_B(A,n) of the area A = \sum_{m=1}^n x_B(m) under such a Lévy bridge and show that, for large n, it has the scaling form P_B(A,n) \sim n^{-1-1/α} F_α(A/n^{1+1/α}), with the asymptotic behavior F_α(Y) \sim Y^{-2(1+α)} for large Y. For α=1, we obtain an explicit expression of F_1(Y) in terms of elementary functions. We also compute the average profile < \tilde x_B (m) > at time m of a Lévy bridge with fixed area A. For large n and large m and A, one finds the scaling form < \tilde x_B(m) > = n^{1/α} H_α({m}/{n},{A}/{n^{1+1/α}}), where at variance with Brownian bridge, H_α(X,Y) is a non trivial function of the rescaled time m/n and rescaled area Y = A/n^{1+1/α}. Our analytical results are verified by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extreme Eigenvalues of Wishart Matrices: Application to Entangled Bipartite System

We discuss an application of the random matrix theory in the context of estimating the bipartite entanglement of a quantum system. We discuss how the Wishart ensemble (the earliest studied random matrix ensemble) appears in this quantum problem. The eigenvalues of the reduced density matrix of one of the subsystems have similar statistical properties as those of the Wishart matrices, except that their {\em trace is constrained to be unity}. We focus here on the smallest eigenvalue which serves as an important measure of entanglement between the two subsystems. In the hard edge case (when the two subsystems have equal sizes) one can fully characterize the probability distribution of the minimum eigenvalue for real, complex and quaternion matrices of all sizes. In particular, we discuss the important finite size effect due to the {\em fixed trace constraint}.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Maximum of N Independent Brownian Walkers till the First Exit From the Half Space

We consider the one-dimensional target search process that involves an immobile target located at the origin and $N$ searchers performing independent Brownian motions starting at the initial positions $\vec x = (x_1,x_2,..., x_N)$ all on the positive half space. The process stops when the target is first found by one of the searchers. We compute the probability distribution of the maximum distance $m$ visited by the searchers till the stopping time and show that it has a power law tail: $P_N(m|\vec x)\sim B_N (x_1x_2... x_N)/m^{N+1}$ for large $m$. Thus all moments of $m$ up to the order $(N-1)$ are finite, while the higher moments diverge. The prefactor $B_N$ increases with $N$ faster than exponentially. Our solution gives the exit probability of a set of $N$ particles from a box $[0,L]$ through the left boundary. Incidentally, it also provides an exact solution of the Laplace&#39;s equation in an $N$-dimensional hypercube with some prescribed boundary conditions. The analytical results are in excellent agreement with Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Switching and growth for microbial populations in catastrophic responsive environments

Phase variation, or stochastic switching between alternative states of gene expression, is common among microbes, and may be important in coping with changing environments. We use a theoretical model to assess whether such switching is a good strategy for growth in environments with occasional catastrophic events. We find that switching can be advantageous, but only when the environment is responsive to the microbial population. In our model, microbes switch randomly between two phenotypic states, with different growth rates. The environment undergoes sudden &#34;catastrophes&#34;, the probability of which depends on the composition of the population. We derive a simple analytical result for the population growth rate. For a responsive environment, two alternative strategies emerge. In the &#34;no switching&#34; strategy, the population maximises its instantaneous growth rate, regardless of catastrophes. In the &#34;switching&#34; strategy, the microbial switching rate is tuned to minimise the environmental response. Which of these strategies is most favourable depends on the parameters of the model. Previous studies have shown that microbial switching can be favourable when the environment changes in an unresponsive fashion between several states. Here, we demonstrate an alternative role for phase variation in allowing microbes to maximise their growth in catastrophic responsive environments.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hitting probability for anomalous diffusion processes

We present the universal features of the hitting probability $Q(x,L)$, the probability that a generic stochastic process starting at $x$ and evolving in a box $[0,L]$ hits the upper boundary $L$ before hitting the lower boundary at 0. For a generic self-affine process (describing, for instance, the polymer translocation through a nanopore) we show that $Q(x,L)=Q(x/L)$ and the scaling function $Q(z)\sim z^ϕ$ as $z\to 0$ with $ϕ=θ/H$ where $H$ and $θ$ are respectively the Hurst exponent and the persistence exponent of the process. This result is verified in several exact calculations including when the process represents the position of a particle diffusing in a disordered potential. We also provide numerical supports for our analytical results.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Random Convex Hulls and Extreme Value Statistics

In this paper we study the statistical properties of convex hulls of $N$ random points in a plane chosen according to a given distribution. The points may be chosen independently or they may be correlated. After a non-exhaustive survey of the somewhat sporadic literature and diverse methods used in the random convex hull problem, we present a unifying approach, based on the notion of support function of a closed curve and the associated Cauchy&#39;s formulae, that allows us to compute exactly the mean perimeter and the mean area enclosed by the convex polygon both in case of independent as well as correlated points. Our method demonstrates a beautiful link between the random convex hull problem and the subject of extreme value statistics. As an example of correlated points, we study here in detail the case when the points represent the vertices of $n$ independent random walks. In the continuum time limit this reduces to $n$ independent planar Brownian trajectories for which we compute exactly, for all $n$, the mean perimeter and the mean area of their global convex hull. Our results have relevant applications in ecology in estimating the home range of a herd of animals. Some of these results were announced recently in a short communication [Phys. Rev. Lett. {\bf 103}, 140602 (2009)].

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Annihilation of Immobile Reactants on the Bethe Lattice

Two-particle annihilation reaction, A+A -> inert, for immobile reactants on the Bethe lattice is solved exactly for the initially random distribution. The process reaches an absorbing state in which no nearest-neighbor reactants are left. The approach of the concentration to the limiting value is exponential. The solution reproduces the known one-dimensional result which is further extended to the reaction A+B -> inert.

0 citations0 reviewsNo code linkedOpen access