Perturbation theory for fractional Brownian motion in presence of absorbing boundaries.

Fractional Brownian motion is a Gaussian process x(t) with zero mean and two-time correlations (x(t(1))x(t(2)))=D(t(1)(2H)+t(2)(2H)-|t(1)-t(2)|(2H)), where H, with 0<H<1, is called the Hurst exponent. For H=1/2, x(t) is a Brownian motion, while for H≠1/2, x(t) is a non-Markovian process. Here we study x(t) in presence of an absorbing boundary at the origin and focus on the probability density P(+)(x,t) for the process to arrive at x at time t, starting near the origin at time 0, given that it has never crossed the origin. It has a scaling form P(+)(x,t)t(-H)R(+)(x/t(H)). Our objective is to compute the scaling function R(+)(y), which up to now was only known for the Markov case H=1/2. We develop a systematic perturbation theory around this limit, setting H=1/2+ε, to calculate the scaling function R(+)(y) to first order in ε. We find that R(+)(y) behaves as R(+)(y)y(ϕ) as y→0 (near the absorbing boundary), while R(+)(y)~y(γ)exp(-y(2)/2) as y→∞, with ϕ=1-4ε+O(ε(2)) and γ=1-2ε+O(ε(2)). Our ε-expansion result confirms the scaling relation ϕ=(1-H)/H proposed in Zoia, Rosso, and Majumdar [Phys. Rev. Lett. 102, 120602 (2009)]. We verify our findings via numerical simulations for H=2/3. The tools developed here are versatile, powerful, and adaptable to different situations.