Researcher profile

D. S. Grebenkov

D. S. Grebenkov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
8topics
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

9 published item(s)

preprint2022arXiv

First-encounter time of two diffusing particles in two- and three-dimensional confinement

The statistics of the first-encounter time of diffusing particles changes drastically when they are placed under confinement. In the present work, we make use of Monte Carlo simulations to study the behavior of a two-particle system in two- and three-dimensional domains with reflecting boundaries. Based on the outcome of the simulations, we give a comprehensive overview of the behavior of the survival probability $S(t)$ and the associated first-encounter time probability density $H(t)$ over a broad time range spanning several decades. In addition, we provide numerical estimates and empirical formulas for the mean first-encounter time $\langle \cal{T}\rangle $, as well as for the decay time $T$ characterizing the monoexponential long-time decay of the survival probability. Based on the distance between the boundary and the center of mass of two particles, we obtain an empirical lower bound $t_B$ for the time at which $S(t)$ starts to significantly deviate from its counterpart for the no boundary case. Surprisingly, for small-sized particles, the dominant contribution to $T$ depends only on the total diffusivity $D=D_1+D_2$, in sharp contrast to the one-dimensional case. This contribution can be related to the Wiener sausage generated by a fictitious Brownian particle with diffusivity $D$. In two dimensions, the first subleading contribution to $T$ is found to depend weakly on the ratio $D_1/D_2$. We also investigate the slow-diffusion limit when $D_2 \ll D_1$ and discuss the transition to the limit when one particle is a fixed target. Finally, we give some indications to anticipate when $T$ can be expected to be a good approximation for $\langle \cal{T}\rangle$.

preprint2020arXiv

First-encounter time of two diffusing particles in confinement

We investigate how confinement may drastically change both the probability density of the first-encounter time and the related survival probability in the case of two diffusing particles. To obtain analytical insights into this problem, we focus on two one-dimensional settings: a half-line and an interval. We first consider the case with equal particle diffusivities, for which exact results can be obtained for the survival probability and the associated first-encounter time density over the full time domain. We also evaluate the moments of the first-encounter time when they exist. We then turn to the case when the diffusivities are not equal, and focus on the long-time behavior of the survival probability. Our results highlight the great impact of boundary effects in diffusion-controlled kinetics even for simple one-dimensional settings, as well as the difficulty of obtaining analytic results as soon as translational invariance of such systems is broken.

preprint2020arXiv

Universal spectral features of different classes of random diffusivity processes

.Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic $1/f^2$-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.

preprint2019arXiv

Full distribution of first exit times in the narrow escape problem

In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small "escape window" in the otherwise impermeable boundary, once it arrives to this window and over-passes an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.

preprint2017arXiv

On spectral properties of the Bloch-Torrey operator in two dimensions

We investigate a two-dimensional Schrödinger operator, $-h^2 Δ+iV(x)$, with a purely complex potential $iV(x)$. A rigorous definition of this non-selfadjoint operator is provided for bounded and unbounded domains with common boundary conditions (Dirichlet, Neumann, Robin and transmission). We propose a general perturbative approach to construct its quasimodes in the semi-classical limit. An alternative WKB construction is also discussed. These approaches are local and thus valid for both bounded and unbounded domains, allowing one to compute the approximate eigenvalues to any order in the small $h$ limit. The general results are further illustrated on the particular case of the Bloch-Torrey operator, $-h^2Δ+ ix_1$, for which a four-term asymptotics is explicitly computed. Its high accuracy is confirmed by a numerical computation of the eigenvalues and eigenfunctions of this operator for a disk and circular annuli. The localization of eigenfunctions near the specific boundary points is revealed. Some applications in the field of diffusion nuclear magnetic resonance are discussed.

preprint2016arXiv

The complex Airy operator with a semi-permeable barrier

We consider a suitable extension of the complex Airy operator, $-d^2/dx^2 + ix$, on the real line with a transmission boundary condition at the origin. We provide a rigorous definition of this operator and study its spectral properties. In particular, we show that the spectrum is discrete, the space generated by the generalized eigenfunctions is dense in $L^2$ (completeness), and we analyze the decay of the associated semi-group. We also present explicit formulas for the integral kernel of the resolvent in terms of Airy functions, investigate its poles, and derive the resolvent estimates.

preprint2016arXiv

The escape problem for mortal walkers

We introduce and investigate the escape problem for random walkers that may eventually die, decay, bleach, or lose activity during their diffusion towards an escape or reactive region on the boundary of a confining domain. In the case of a first-order kinetics (i.e., exponentially distributed lifetimes), we study the effect of the associated death rate onto the survival probability, the exit probability, and the mean first passage time. We derive the upper and lower bounds and some approximations for these quantities. We reveal three asymptotic regimes of small, intermediate and large death rates. General estimates and asymptotics are compared to several explicit solutions for simple domains, and to numerical simulations. These results allow one to account for stochastic photobleaching of fluorescent tracers in bio-imaging, degradation of mRNA molecules in genetic translation mechanisms, or high mortality rates of spermatozoa in the fertilization process. This is also a mathematical ground for optimizing storage containers and materials to reduce the risk of leakage of dangerous chemicals or nuclear wastes.

preprint2013arXiv

Following a Trend with an Exponential Moving Average: Analytical Results for a Gaussian Model

We investigate how price variations of a stock are transformed into profits and losses (P&Ls) of a trend following strategy. In the frame of a Gaussian model, we derive the probability distribution of P&Ls and analyze its moments (mean, variance, skewness and kurtosis) and asymptotic behavior (quantiles). We show that the asymmetry of the distribution (with often small losses and less frequent but significant profits) is reminiscent to trend following strategies and less dependent on peculiarities of price variations. At short times, trend following strategies admit larger losses than one may anticipate from standard Gaussian estimates, while smaller losses are ensured at longer times. Simple explicit formulas characterizing the distribution of P&Ls illustrate the basic mechanisms of momentum trading, while general matrix representations can be applied to arbitrary Gaussian models. We also compute explicitly annualized risk adjusted P&L and strategy turnover to account for transaction costs. We deduce the trend following optimal timescale and its dependence on both auto-correlation level and transaction costs. Theoretical results are illustrated on the Dow Jones index.

preprint2013arXiv

Optimal and sub-optimal quadratic forms for non-centered Gaussian processes

Individual random trajectories of stochastic processes are often analyzed by using quadratic forms such as time averaged (TA) mean square displacement (MSD) or velocity auto-correlation function (VACF). The appropriate quadratic form is expected to have a narrow probability distribution in order to reduce statistical uncertainty of a single measurement. We consider the problem of finding the optimal quadratic form that minimizes a chosen cumulant moment (e.g., the variance) of the probability distribution, under the constraint of fixed mean value. For discrete non-centered Gaussian processes, we construct the optimal quadratic form by using the spectral representation of the cumulant moments. Moreover, we obtain a simple explicit formula for the smallest achievable cumulant moment that may serve as a quality benchmark for other quadratic forms. We illustrate the optimality issues by comparing the optimal variance with the variances of the TA MSD and TA VACF of fractional Brownian motion superimposed with a constant drift and independent Gaussian noise.