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 ~ t^{2H} + s^{2H} - |t-s|^{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 unequal 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) ~ R(x/t^H)/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 + epsilon, to calculate the scaling function R(y) to first order in epsilon. We find that R(y) behaves as R(y) ~ y^phi as y -> 0 (near the absorbing boundary), while R(y) ~ y^gamma exp(-y^2/2) as y -> oo, with phi = 1 - 4 epsilon + O(epsilon^2) and gamma = 1 - 2 epsilon + O(epsilon^2). Our epsilon-expansion result confirms the scaling relation phi = (1-H)/H proposed in Ref. [28]. We verify our findings via numerical simulations for H = 2/3. The tools developed here are versatile, powerful, and adaptable to different situations.Comment: 16 pages, 8 figures; revised version 2 adds discussion on spatial small-distance cutof

II. Regnskabsoversigt for Aarene 1849—68 fra Gunderslevholm Skovdistrikt.

II. Regnskabsoversigt for Aarene 1849—68 fra Gunderslevholm Skovdistrikt

Survival of a Diffusing Particle in a Transverse Shear Flow: A First-Passage Problem with Continuously Varying Persistence Exponent

We consider a particle diffusing in the y-direction, dy/dt=\eta(t), subject to a transverse shear flow in the x-direction, dx/dt=f(y), where x \ge 0 and x=0 is an absorbing boundary. We treat the class of models defined by f(y) = \pm v_{\pm}(\pm y)^\alpha where the upper (lower) sign refers to y>0 (y<0). We show that the particle survives with probability Q(t) \sim t^{-\theta} with \theta = 1/4, independent of \alpha, if v_{+}=v_{-}. If v_{+} \ne v_{-}, however, we show that \theta depends on both \alpha and the ratio v_{+}/v_{-}, and we determine this dependence.Comment: 4 page

Comment on "Mean First Passage Time for Anomalous Diffusion"

We correct a previously erroneous calculation [Phys. Rev. E 62, 6065 (2000)] of the mean first passage time of a subdiffusive process to reach either end of a finite interval in one dimension. The mean first passage time is in fact infinite.Comment: To appear in Phys. Rev.

Constraining SIDM with halo shapes: Revisited predictions from realistic simulations of early-type galaxies

We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs (Mvir sim 10 13 , M ⊙ and M ∗ ∼ 10 11 , M ⊙), simulated in CDM and a SIDM model with constant cross-section of σT/mχ = 1 cm2g-1. We measure the 3D and projected shapes of the dark matter haloes and their baryonic content using the inertia tensor and compare our measurements to the results of three HST samples of gravitational lenses and Chandra and XMM-Newton X-ray observations. We find that the inclusion of baryons greatly reduces the differences between CDM and a SIDM, together with the ability to draw constraints based on shapes. Lensing measurements reject the predictions of CDM dark-matter-only simulations and prefer one of the hydro scenarios. When we consider the total sample of lenses, observational data prefer the CDM hydro scenario. The shapes of the X-ray emitting gas are compatible with observational results in both hydro runs, with CDM predicting higher elongations only in the very centre. Contrary to previous claims at the scale of elliptical galaxies, we conclude that both CDM and our SIDM model can still explain observed shapes once we include baryons in the simulations. Our results demonstrate that this is essential to derive realistic constraints and that new simulations are needed to confirm and extend our findings

Anomalous diffusion and generalized Sparre-Andersen scaling

We are discussing long-time, scaling limit for the anomalous diffusion composed of the subordinated L\'evy-Wiener process. The limiting anomalous diffusion is in general non-Markov, even in the regime, where ensemble averages of a mean-square displacement or quantiles representing the group spread of the distribution follow the scaling characteristic for an ordinary stochastic diffusion. To discriminate between truly memory-less process and the non-Markov one, we are analyzing deviation of the survival probability from the (standard) Sparre-Andersen scaling.Comment: 5 pages, 3 figure

Record statistics and persistence for a random walk with a drift

We study the statistics of records of a one-dimensional random walk of n steps, starting from the origin, and in presence of a constant bias c. At each time-step the walker makes a random jump of length \eta drawn from a continuous distribution f(\eta) which is symmetric around a constant drift c. We focus in particular on the case were f(\eta) is a symmetric stable law with a L\'evy index 0 < \mu \leq 2. The record statistics depends crucially on the persistence probability which, as we show here, exhibits different behaviors depending on the sign of c and the value of the parameter \mu. Hence, in the limit of a large number of steps n, the record statistics is sensitive to these parameters (c and \mu) of the jump distribution. We compute the asymptotic mean record number after n steps as well as its full distribution P(R,n). We also compute the statistics of the ages of the longest and the shortest lasting record. Our exact computations show the existence of five distinct regions in the (c, 0 < \mu \leq 2) strip where these quantities display qualitatively different behaviors. We also present numerical simulation results that verify our analytical predictions.Comment: 51 pages, 22 figures. Published version (typos have been corrected

Asymptotic behavior of self-affine processes in semi-infinite domains

We propose to model the stochastic dynamics of a polymer passing through a pore (translocation) by means of a fractional Brownian motion, and study its behavior in presence of an absorbing boundary. Based on scaling arguments and numerical simulations, we present a conjecture that provides a link between the persistence exponent $\theta$ and the Hurst exponent $H$ of the process, thus sheding light on the spatial and temporal features of translocation. Furthermore, we show that this conjecture applies more generally to a broad class of self affine processes undergoing anomalous diffusion in bounded domains, and we discuss some significant examples.Comment: 4 pages, 3 figures; to be published in Phys. Rev. Let

Record statistics for biased random walks, with an application to financial data

We consider the occurrence of record-breaking events in random walks with asymmetric jump distributions. The statistics of records in symmetric random walks was previously analyzed by Majumdar and Ziff and is well understood. Unlike the case of symmetric jump distributions, in the asymmetric case the statistics of records depends on the choice of the jump distribution. We compute the record rate $P_n(c)$, defined as the probability for the $n$th value to be larger than all previous values, for a Gaussian jump distribution with standard deviation $\sigma$ that is shifted by a constant drift $c$. For small drift, in the sense of $c/\sigma \ll n^{-1/2}$, the correction to $P_n(c)$ grows proportional to arctan$(\sqrt{n})$ and saturates at the value $\frac{c}{\sqrt{2} \sigma}$. For large $n$ the record rate approaches a constant, which is approximately given by $1-(\sigma/\sqrt{2\pi}c)\textrm{exp}(-c^2/2\sigma^2)$ for $c/\sigma \gg 1$. These asymptotic results carry over to other continuous jump distributions with finite variance. As an application, we compare our analytical results to the record statistics of 366 daily stock prices from the Standard & Poors 500 index. The biased random walk accounts quantitatively for the increase in the number of upper records due to the overall trend in the stock prices, and after detrending the number of upper records is in good agreement with the symmetric random walk. However the number of lower records in the detrended data is significantly reduced by a mechanism that remains to be identified.Comment: 16 pages, 7 figure

Record Statistics for Multiple Random Walks

We study the statistics of the number of records R_{n,N} for N identical and independent symmetric discrete-time random walks of n steps in one dimension, all starting at the origin at step 0. At each time step, each walker jumps by a random length drawn independently from a symmetric and continuous distribution. We consider two cases: (I) when the variance \sigma^2 of the jump distribution is finite and (II) when \sigma^2 is divergent as in the case of L\'evy flights with index 0 < \mu < 2. In both cases we find that the mean record number grows universally as \sim \alpha_N \sqrt{n} for large n, but with a very different behavior of the amplitude \alpha_N for N > 1 in the two cases. We find that for large N, \alpha_N \approx 2 \sqrt{\log N} independently of \sigma^2 in case I. In contrast, in case II, the amplitude approaches to an N-independent constant for large N, \alpha_N \approx 4/\sqrt{\pi}, independently of 0<\mu<2. For finite \sigma^2 we argue, and this is confirmed by our numerical simulations, that the full distribution of (R_{n,N}/\sqrt{n} - 2 \sqrt{\log N}) \sqrt{\log N} converges to a Gumbel law as n \to \infty and N \to \infty. In case II, our numerical simulations indicate that the distribution of R_{n,N}/\sqrt{n} converges, for n \to \infty and N \to \infty, to a universal nontrivial distribution, independently of \mu. We discuss the applications of our results to the study of the record statistics of 366 daily stock prices from the Standard & Poors 500 index.Comment: 25 pages, 8 figure