Fractional Brownian motion in confining potentials: non-equilibrium distribution tails and optimal fluctuations

TL;DR

This paper uses the optimal fluctuation method to analyze the large-distance tails of the steady-state distribution of a fractional Brownian particle in power-law potentials, revealing non-equilibrium features and confinement conditions.

Contribution

It provides an analytical and numerical framework to determine the tail behavior of the distribution for fractional Brownian motion in confining potentials, including explicit results for the fractional Ornstein-Uhlenbeck process.

Findings

Derived the large-|x| tail behavior of the steady-state distribution.

Established conditions for particle confinement based on H and m.

Developed a numerical algorithm for optimal path computation.

Abstract

At long times, a fractional Brownian particle in a confining external potential reaches a non-equilibrium (non-Boltzmann) steady state. Here we consider scale-invariant power-law potentials $V(x)\sim |x|^m$, where $m>0$, and employ the optimal fluctuation method (OFM) to determine the large-$|x|$ tails of the steady-state probability distribution $\mathcal{P}(x)$ of the particle position. The calculations involve finding the optimal (that is, the most likely) path of the particle, which determines these tails, via a minimization of the exact action functional for this system, which has recently become available. Exploiting dynamical scale invariance of the model in conjunction with the OFM ansatz, we establish the large-$|x|$ tails of $\ln \mathcal{P}(x)$ up to a dimensionless factor $\alpha(H,m)$, where $0<H<1$ is the Hurst exponent. We determine $\alpha(H,m)$ analytically (i) in the limits of $H\to 0$ and $H\to 1$, and (ii) for $m=2$ and arbitrary $H$, corresponding to the fractional Ornstein-Uhlenbeck (fOU) process. Our results for the fOU process are in agreement with the previously known exact $\mathcal{P}(x)$ and autocovariance. The form of the tails of $\mathcal{P}(x)$ yields exact conditions, in terms of $H$ and $m$, for the particle confinement in the potential. For $H\neq 1/2$, the tails encode the non-equilibrium character of the steady state distribution, and we observe violation of time reversibility of the system except for $m=2$. To compute the optimal paths and the factor $\alpha(H,m)$ for arbitrary permissible $H$ and $m$, one needs to solve an (in general nonlinear) integro-differential equation. To this end we develop a specialized numerical iteration algorithm which accounts analytically for an intrinsic cusp singularity of the optimal paths for $H<1/2$.