Large-deviation functions for nonlinear functionals of a Gaussian stationary Markov process

A study of anomalous stochastic processes via generalizing fractional calculus

Due to the very importance of fractional calculus in studying anomalous stochastic processes, we systematically investigate the existing formulation of fractional calculus and generalize it to broader applied contexts. Specifically, based on the improved Riemann-Liouville fractional calculus operators and the modified Maruyama's notation for fractional Brownian motion, we develop the fractional Ito^'s calculus and derive a generalized Fokker-Planck equation corresponding to the Maruyama's process, along with which, the stochastic realizations of trajectories, both underdamped and overdamped, have been studied in terms of the stochastic dynamics equations newly formulated. This paves a way to study the path integrals and the stochastic thermodynamics of anomalous stochastic processes. We also explicitly derive several fundamental results in fractional calculus, including the relation between fractional and normal differentiation, the Laplace transform for fractional derivatives, the analytic solution of one type of generalized diffusion equations, and the fractional integration formulas. Our results advance the existing fractional calculus and provide practical references for studying anomalous diffusion, mechanics of memory materials in engineering, and stochastic analysis in fractional orders.

Occupation time of a system of Brownian particles on the line with steplike initial condition

We consider a system of noninteracting Brownian particles on the line with steplike initial condition and study the statistics of the occupation time on the positive half-line. We demonstrate that even at large times, the behavior of the occupation time exhibits long-lasting memory effects of the initialization. Specifically, we calculate the mean and the variance of the occupation time, demonstrating that the memory effects in the variance are determined by a generalized compressibility (or Fano factor), associated with the initial condition. In the particular case of the uncorrelated uniform initial condition we conduct a detailed study of two probability distributions of the occupation time: annealed (averaged over all possible initial configurations) and quenched (for a typical configuration). We show that at large times both the annealed and the quenched distributions admit large deviation form and we compute analytically the associated rate functions. We verify our analytical predictions via numerical simulations using importance sampling Monte Carlo strategy.

Nonequilibrium steady state of trapped active particles

We consider an overdamped particle with a general physical mechanism that creates noisy active movement (e.g., a run-and-tumble particle or active Brownian particle, etc.), that is confined by an external potential. Focusing on the limit in which the correlation time τ of the active noise is small, we find the nonequilibrium steady-state distribution P_{st}(X) of the particle's position X. While typical fluctuations of X follow a Boltzmann distribution with an effective temperature that is not difficult to find, the tails of P_{st}(X) deviate from a Boltzmann behavior: In the limit τ→0, they scale as P_{st}(X)∼e^{-s(X)/τ}. We calculate the large-deviation function s(X) exactly for arbitrary trapping potential and active noise in dimension d=1, by relating it to the rate function that describes large deviations of the position of the same active particle in absence of an external potential at long times. We then extend our results to d>1 assuming rotational symmetry.

Dynamical phase transition in the occupation fraction statistics for noncrossing Brownian particles

We consider a system of N noncrossing Brownian particles in one dimension. We find the exact rate function that describes the long-time large deviation statistics of their occupation fraction in a finite interval in space. Remarkably, we find that, for any general N≥2, the system undergoes N-1 dynamical phase transitions of second order. The N-1 transitions are the boundaries of N phases that correspond to different numbers of particles which are in the vicinity of the interval throughout the dynamics. We achieve this by mapping the problem to that of finding the ground-state energy for N noninteracting spinless fermions in a square-well potential. The phases correspond to different numbers of single-body bound states for the quantum problem. We also study the process conditioned on a given occupation fraction and the large-N limiting behavior.

Nonequilibrium steady state for harmonically confined active particles

We study the full nonequilibrium steady-state distribution P_{st}(X) of the position X of a damped particle confined in a harmonic trapping potential and experiencing active noise whose correlation time τ_{c} is assumed to be very short. Typical fluctuations of X are governed by a Boltzmann distribution with an effective temperature that is found by approximating the noise as white Gaussian thermal noise. However, large deviations of X are described by a non-Boltzmann steady-state distribution. We find that, in the limit τ_{c}→0, they display the scaling behavior P_{st}(X)∼e^{-s(X)/τ_{c}}, where s(X) is the large-deviation function. We obtain an expression for s(X) for a general active noise and calculate it exactly for the particular case of telegraphic (dichotomous) noise.

Large deviations in chaotic systems: Exact results and dynamical phase transition

Large deviations in chaotic dynamics have potentially significant and dramatic consequences. We study large deviations of series of finite lengths N generated by chaotic maps. The distributions generally display an exponential decay with N, associated with large-deviation (rate) functions. We obtain the exact rate functions analytically for the doubling, tent, and logistic maps. For the latter two, the solution is given as a power series whose coefficients can be systematically calculated to any order. We also obtain the rate function for the cat map numerically, uncovering strong evidence for the existence of a remarkable singularity of it that we interpret as a second-order dynamical phase transition. Furthermore, we develop a numerical tool for efficiently simulating atypical realizations of sequences if the chaotic map is not invertible, and we apply it to the tent and logistic maps.

Stiffness of random walks with reflecting boundary conditions

We study the distribution of occupation times for a one-dimensional random walk restricted to a finite interval by reflecting boundary conditions. At short times the classical bimodal distribution due to Lévy is reproduced with walkers staying mostly either to the left or right of the initial point. With increasing time, however, the boundaries suppress large excursions from the starting point, and the distribution becomes unimodal, converging to a δ distribution in the long-time limit. An approximate spectral analysis of the underlying Fokker-Planck equation yields results in excellent agreement with numerical simulations.

Anomalous scaling and first-order dynamical phase transition in large deviations of the Ornstein-Uhlenbeck process

We study the full distribution of A=∫_{0}^{T}x^{n}(t)dt, n=1,2,⋯, where x(t) is an Ornstein-Uhlenbeck process. We find that for n>2 the long-time (T→∞) scaling form of the distribution is of the anomalous form P(A;T)∼e^{-T^{μ}f_{n}(ΔA/T^{ν})} where ΔA is the difference between A and its mean value, and the anomalous exponents are μ=2/(2n-2) and ν=n/(2n-2). The rate function f_{n}(y), which we calculate exactly, exhibits a first-order dynamical phase transition which separates between a homogeneous phase that describes the Gaussian distribution of typical fluctuations, and a "condensed" phase that describes the tails of the distribution. We also calculate the most likely realizations of A(t)=∫_{0}^{t}x^{n}(s)ds and the distribution of x(t) at an intermediate time t conditioned on a given value of A. Extensions and implications to other continuous-time systems are discussed.

Effective Langevin equations leading to large deviation function of time-averaged velocity for a nonequilibrium Rayleigh piston

We study fluctuating dynamics of a freely movable piston that separates an infinite cylinder into two regions filled with ideal gas particles at the same pressure but different temperatures. To investigate statistical properties of the time-averaged velocity of the piston in the long-time limit, we perturbatively calculate the large deviation function of the time-averaged velocity. Then, we derive an infinite number of effective Langevin equations yielding the same large deviation function as in the original model. Finally, we provide two possibilities for uniquely determining the form of the effective model.

Anomalous scaling of dynamical large deviations of stationary Gaussian processes

Employing the optimal fluctuation method, we study the large deviation function of long-time averages (1/T)∫_{-T/2}^{T/2}x^{n}(t)dt,n=1,2,⋯, of centered stationary Gaussian processes. These processes are correlated and, in general, non-Markovian. We show that the anomalous scaling with time of the large-deviation function, recently observed for n>2 for the particular case of the Ornstein-Uhlenbeck process, holds for a whole class of stationary Gaussian processes.

Anomalous Scaling of Dynamical Large Deviations

The typical values and fluctuations of time-integrated observables of nonequilibrium processes driven in steady states are known to be characterized by large deviation functions, generalizing the entropy and free energy to nonequilibrium systems. The definition of these functions involves a scaling limit, similar to the thermodynamic limit, in which the integration time τ appears linearly, unless the process considered has long-range correlations, in which case τ is generally replaced by τ^{ξ} with ξ≠1. Here, we show that such an anomalous power-law scaling in time of large deviations can also arise without long-range correlations in Markovian processes as simple as the Langevin equation. We describe the mechanism underlying this scaling using path integrals and discuss its physical consequences for more general processes.

Large deviations for Markov processes with resetting

Markov processes restarted or reset at random times to a fixed state or region in space have been actively studied recently in connection with random searches, foraging, and population dynamics. Here we study the large deviations of time-additive functions or observables of Markov processes with resetting. By deriving a renewal formula linking generating functions with and without resetting, we are able to obtain the rate function of such observables, characterizing the likelihood of their fluctuations in the long-time limit. We consider as an illustration the large deviations of the area of the Ornstein-Uhlenbeck process with resetting. Other applications involving diffusions, random walks, and jump processes with resetting or catastrophes are discussed.

Symmetry for the duration of entropy-consuming intervals

We introduce the violation fraction υ as the cumulative fraction of time that a mesoscopic system spends consuming entropy at a single trajectory in phase space. We show that the fluctuations of this quantity are described in terms of a symmetry relation reminiscent of fluctuation theorems, which involve a function Φ, which can be interpreted as an entropy associated with the fluctuations of the violation fraction. The function Φ, when evaluated for arbitrary stochastic realizations of the violation fraction, is odd upon the symmetry transformations that are relevant for the associated stochastic entropy production. This fact leads to a detailed fluctuation theorem for the probability density function of Φ. We study the steady-state limit of this symmetry in the paradigmatic case of a colloidal particle dragged by optical tweezers through an aqueous solution. Finally, we briefly discuss possible applications of our results for the estimation of free-energy differences from single-molecule experiments.

Fractional Feynman-Kac equation for non-brownian functionals

We derive backward and forward fractional Feynman-Kac equations for the distribution of functionals of the path of a particle undergoing anomalous diffusion. Fractional substantial derivatives introduced by Friedrich and co-workers [Phys. Rev. Lett. 96, 230601 (2006)10.1103/PhysRevLett.96.230601] provide the correct fractional framework for the problem. For applications, we calculate the distribution of occupation times in half space and show how the statistics of anomalous functionals is related to weak ergodicity breaking.

Statistical properties of functionals of the paths of a particle diffusing in a one-dimensional random potential

We present a formalism for obtaining the statistical properties of functionals and inverse functionals of the paths of a particle diffusing in a one-dimensional quenched random potential. We demonstrate the implementation of the formalism in two specific examples: (1) where the functional corresponds to the local time spent by the particle around the origin and (2) where the functional corresponds to the occupation time spent by the particle on the positive side of the origin, within an observation time window of size t. We compute the disorder average distributions of the local time, the inverse local time, the occupation time, and the inverse occupation time and show that in many cases disorder modifies the behavior drastically.

Persistence exponents and the statistics of crossings and occupation times for Gaussian stationary processes

We consider the persistence probability, the occupation-time distribution, and the distribution of the number of zero crossings for discrete or (equivalently) discretely sampled Gaussian stationary processes (GSPs) of zero mean. We first consider the Ornstein-Uhlenbeck process, finding expressions for the mean and variance of the number of crossings and the "partial survival" probability. We then elaborate on the correlator expansion developed in an earlier paper [G. C. M. A. Ehrhardt and A. J. Bray, Phys. Rev. Lett. 88, 070602 (2002)] to calculate discretely sampled persistence exponents of GSPs of known correlator by means of a series expansion in the correlator. We apply this method to the processes d(n)x/dt(n)=eta(t) with n>/=3, incorporating an extrapolation of the series to the limit of continuous sampling. We then extend the correlator method to calculate the occupation-time and crossing-number distributions, as well as their partial-survival distributions and the means and variances of the occupation time and number of crossings. We apply these general methods to the d(n)x/dt(n)=eta(t) processes for n=1 (random walk), n=2 (random acceleration), and larger n, and to simple diffusion from random initial conditions in one to three dimensions. The results for discrete sampling are extrapolated to the continuum limit where possible.

Exact occupation time distribution in a non-Markovian sequence and its relation to spin glass models

We compute exactly the distribution of the occupation time in a discrete non-Markovian toy sequence that appears in various physical contexts such as the diffusion processes and Ising spin glass chains. The non-Markovian property makes the results nontrivial even for this toy sequence. The distribution is shown to have non-Gaussian tails characterized by a nontrivial large deviation function that is computed explicitly. An exact mapping of this sequence to an Ising spin glass chain via a gauge transformation raises an interesting question for a generic finite sized spin glass model; at a given temperature, what is the distribution (over disorder) of the thermally averaged number of spins that are aligned to their local fields? We show that this distribution remains nontrivial even at infinite temperature and can be computed explicitly in few cases such as in the Sherrington-Kirkpatrick model with Gaussian disorder.