Source author record

Denis S. Grebenkov

Denis S. Grebenkov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

The Effective Reactivity for Capturing Brownian Motion by Partially Reactive Patches on a Spherical Surface

We analyze the trapping of diffusing ligands, modeled as Brownian particles, by a sphere that has $N$ partially reactive boundary patches, each of small area and arbitrary shape, on an otherwise reflecting boundary. For such a structured target, the partial reactivity of each boundary patch is characterized by a Robin boundary condition, with a local boundary reactivity $κ_i$ for $i=1,\ldots,N$. For any spatial arrangement of well-separated patches on the surface of the sphere, the method of matched asymptotic expansions is used to derive explicit results for the capacitance $C_{\rm T}$ of the structured target, which is valid for any $κ_i>0$. This target capacitance $C_{\rm T}$ is defined in terms of a Green's matrix, which depends on the spatial configuration of patches, the local reactive capacitance $C_i(κ_i)$ of each patch and another coefficient that depends on the local geometry near a patch. The analytical dependence of $C_{i}(κ_i)$ on $κ_i$ is uncovered via a spectral expansion over Steklov eigenfunctions. For circular patches, the latter are readily computed numerically and provide an accurate fully explicit sigmoidal approximation for $C_{i}(κ_i)$. In the homogenization limit of $N\gg 1$ identical uniformly-spaced patches with $κ_i=κ$, we derive an explicit scaling law for the effective capacitance and the effective reactivity of the structured target that is valid in the limit of small patch area fraction. From a comparison with numerical simulations, we show that this scaling law provides a highly accurate approximation over the full range $κ>0$, even when there is only a moderately large number of reactive patches.

preprint2022arXiv

Depletion of Resources by a Population of Diffusing Species

Depletion of natural and artificial resources is a fundamental problem and a potential cause of economic crises, ecological catastrophes, and death of living organisms. Understanding the depletion process is crucial for its further control and optimized replenishment of resources. In this paper, we investigate a stock depletion by a population of species that undergo an ordinary diffusion and consume resources upon each encounter with the stock. We derive the exact form of the probability density of the random depletion time, at which the stock is exhausted. The dependence of this distribution on the number of species, the initial amount of resources, and the geometric setting is analyzed. Future perspectives and related open problems are discussed.

preprint2022arXiv

Mean first-passage time to a small absorbing target in three-dimensional elongated domains

We derive an approximate formula for the mean first-passage time (MFPT) to a small absorbing target of arbitrary shape inside an elongated domain of a slowly varying axisymmetric profile. For this purpose, the original Poisson equation in three dimensions is reduced an effective one-dimensional problem on an interval with a semi-permeable semi-absorbing membrane. The approximate formula captures correctly the dependence of the MFPT on the distance to the target, the radial profile of the domain, and the size and the shape of the target. This approximation is validated by Monte Carlo simulations.

preprint2020arXiv

Paradigm shift in diffusion-mediated surface phenomena

Diffusion-mediated surface phenomena are crucial for human life and industry, with examples ranging from oxygen capture by lung alveolar surface to heterogeneous catalysis, gene regulation, membrane permeation and filtration processes. Their current description via diffusion equations with mixed boundary conditions is limited to simple surface reactions with infinite or constant reactivity. In this letter, we propose a probabilistic approach based on the concept of boundary local time to investigate the intricate dynamics of diffusing particles near a reactive surface. Reformulating surface-particle interactions in terms of stopping conditions, we obtain in a unified way major diffusion-reaction characteristics such as the propagator, the survival probability, the first-passage time distribution, and the reaction rate. This general formalism allows us to describe new surface reaction mechanisms such as for instance surface reactivity depending on the number of encounters with the diffusing particle that can model the effects of catalyst fooling or membrane degradation. The disentanglement of the geometric structure of the medium from surface reactivity opens far-reaching perspectives for modeling, optimization and control of diffusion-mediated surface phenomena.

preprint2020arXiv

Surface Hopping Propagator: An Alternative Approach to Diffusion-Influenced Reactions

Dynamics of a particle diffusing in a confinement can be seen a sequence of bulk-diffusion-mediated hops on the confinement surface. Here, we investigate the surface hopping propagator that describes the position of the diffusing particle after a prescribed number of encounters with that surface. This quantity plays the central role in diffusion-influenced reactions and determines their most common characteristics such as the propagator, the first-passage time distribution, and the reaction rate. We derive explicit formulas for the surface hopping propagator and related quantities for several Euclidean domains: half-space, circular annuli, circular cylinders, and spherical shells. These results provide the theoretical ground for studying diffusion-mediated surface phenomena. The behavior of the surface hopping propagator is investigated for both "immortal" and "mortal" particles.

preprint2019arXiv

Probability distribution of the boundary local time of reflected Brownian motion in Euclidean domains

How long does a diffusing molecule spend in a close vicinity of a confining boundary or a catalytic surface? This quantity is determined by the boundary local time, which plays thus a crucial role in the description of various surface-mediated phenomena such as heterogeneous catalysis, permeation through semi-permeable membranes, or surface relaxation in nuclear magnetic resonance. In this paper, we obtain the probability distribution of the boundary local time in terms of the spectral properties of the Dirichlet-to-Neumann operator. We investigate the short-time and long-time asymptotic behaviors of this random variable for both bounded and unbounded domains. This analysis provides complementary insights onto the dynamics of diffusing molecules near partially reactive boundaries.

preprint2017arXiv

Mean perimeter and mean area of the convex hull over planar random walks

We investigate the geometric properties of the convex hull over $n$ successive positions of a planar random walk, with a symmetric continuous jump distribution. We derive the large $n$ asymptotic behavior of the mean perimeter. In addition, we compute the mean area for the particular case of isotropic Gaussian jumps. While the leading terms of these asymptotics are universal, the subleading (correction) terms depend on finer details of the jump distribution and describe a "finite size effect" of discrete-time jump processes, allowing one to accurately compute the mean perimeter and the mean area even for small $n$, as verified by Monte Carlo simulations. This is particularly valuable for applications dealing with discrete-time jumps processes and ranging from the statistical analysis of single-particle tracking experiments in microbiology to home range estimations in ecology.

preprint2016arXiv

Universal formula for the mean first passage time in planar domains

We derive a general exact formula for the mean first passage time (MFPT) from a fixed point inside a planar domain to an escape region on its boundary. The underlying mixed Dirichlet-Neumann boundary value problem is conformally mapped onto the unit disk, solved exactly, and mapped back. The resulting formula for the MFPT is valid for an arbitrary space-dependent diffusion coefficient, while the leading logarithmic term is explicit, simple, and remarkably universal. In contrast to earlier works, we show that the natural small parameter of the problem is the harmonic measure of the escape region, not its perimeter. The conventional scaling of the MFPT with the area of the domain is altered when diffusing particles are released near the escape region. These findings change the current view of escape problems and related chemical or biochemical kinetics in complex, multiscale, porous or fractal domains, while the fundamental relation to the harmonic measure opens new ways of computing and interpreting MFPTs.

preprint2015arXiv

A Multiscale Guide to Brownian Motion

We revise the Levy's construction of Brownian motion as a simple though still rigorous approach to operate with various Gaussian processes. A Brownian path is explicitly constructed as a linear combination of wavelet-based "geometrical features" at multiple length scales with random weights. Such a wavelet representation gives a closed formula mapping of the unit interval onto the functional space of Brownian paths. This formula elucidates many classical results about Brownian motion (e.g., non-differentiability of its path), providing intuitive feeling for non-mathematicians. The illustrative character of the wavelet representation, along with the simple structure of the underlying probability space, is different from the usual presentation of most classical textbooks. Similar concepts are discussed for fractional Brownian motion, Ornstein-Uhlenbeck process, Gaussian free field, and fractional Gaussian fields. Wavelet representations and dyadic decompositions form the basis of many highly efficient numerical methods to simulate Gaussian processes and fields, including Brownian motion and other diffusive processes in confining domains.

preprint2015arXiv

Anomalous, non-Gaussian tracer diffusion in heterogeneously crowded environments

A topic of intense current investigation pursues the question how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain (nucleus), or at the outer boundary (cell membrane), exhibiting distinctly different behaviour of the observed anomalous diffusion for heterogeneous crowder distributions. Specifically we reveal an extremely asymmetric spreading of the tracer even at moderate crowding fractions. In addition to the standard mean squared displacement and the local diffusion exponent of the tracer particles we investigate the magnitude and the amplitude scatter of the time averaged mean squared displacement of individual trajectories, the non-Gaussianity parameter, and the van Hove correlation function of the particle displacements. We also quantify how the average tracer diffusivity varies with the position in the domain with heterogeneous radial distribution of the crowders and examine the behaviour of the survival probability and the dynamics of first passage events of the tracer. Finally, we discuss the relevance of our results to single particle tracking measurements in biological cells.

preprint2014arXiv

Generalized Langevin Equation with Hydrodynamic Backflow: Equilibrium Properties

We review equilibrium properties for the dynamics of a single particle evolving in a visco--elastic medium under the effect of hydrodynamic backflow which includes added mass and Basset force. Arbitrary equilibrium forces acting upon the particle are also included. We discuss the derivation of the explicit expression for the thermal noise correlation function that is consistent with the fluctuation-dissipation theorem. We rely on general time-reversal arguments that apply irrespective of the external potential acting on the particle, but also allow one to retrieve existing results derived for free particles and particles in a harmonic trap. Some consequences for the analysis and interpretation of single-particle tracking experiments are briefly discussed.

preprint2014arXiv

Optimal Allocation of Trend Following Strategies

We consider a portfolio allocation problem for trend following (TF) strategies on multiple correlated assets. Under simplifying assumptions of a Gaussian market and linear TF strategies, we derive analytical formulas for the mean and variance of the portfolio return. We construct then the optimal portfolio that maximizes risk-adjusted return by accounting for inter-asset correlations. The dynamic allocation problem for $n$ assets is shown to be equivalent to the classical static allocation problem for $n^2$ virtual assets that include lead-lag corrections in positions of TF strategies. The respective roles of asset auto-correlations and inter-asset correlations are investigated in depth for the two-asset case and a sector model. In contrast to the principle of diversification suggesting to treat uncorrelated assets, we show that inter-asset correlations allow one to estimate apparent trends more reliably and to adjust the TF positions more efficiently. If properly accounted for, inter-asset correlations are not deteriorative but beneficial for portfolio management that can open new profit opportunities for trend followers.

preprint2014arXiv

Overlap of two Brownian trajectories: exact results for scaling functions

We consider two random walkers starting at the same time $t=0$ from different points in space separated by a given distance $R$. We compute the average volume of the space visited by both walkers up to time $t$ as a function of $R$ and $t$ and dimensionality of space $d$. For $d<4$, this volume, after proper renormalization, is shown to be expressed through a scaling function of a single variable $R/\sqrt{t}$. We provide general integral formulas for scaling functions for arbitrary dimensionality $d<4$. In contrast, we show that no scaling function exists for higher dimensionalities $d \geq 4$

preprint2013arXiv

Evidence of hydrodynamic and subdiffusive motion of tracers in a viscoelastic medium

We propose a theoretical model which relies on the generalized Langevin equation and may account for various dynamical features of the thermal motion of organelles, vesicles or macromolecules in viscoelastic media such as polymer networks. In particular, we consider inertial and hydrodynamic effects at short times, subdiffusive scaling at intermediate times, and eventually optical trapping at long times. Simple analytical formulas for the mean square displacement and velocity auto-correlation function are derived. The developed theory is applied to the analysis of fifty-second long trajectories of micron-sized spherical tracers in actin gels that were acquired at one microsecond temporal resolution by using optical tweezers single-particle tracking. For the first time, both the subdiffusive scaling and hydrodynamic effects are observed within a single experiment and accurately described by a minimal phenomenological model.

preprint2013arXiv

Geometrical structure of Laplacian eigenfunctions

We summarize the properties of eigenvalues and eigenfunctions of the Laplace operator in bounded Euclidean domains with Dirichlet, Neumann or Robin boundary condition. We keep the presentation at a level accessible to scientists from various disciplines ranging from mathematics to physics and computer sciences. The main focus is put onto multiple intricate relations between the shape of a domain and the geometrical structure of eigenfunctions.

preprint2013arXiv

On the exponential decay of Laplacian eigenfunctions in planar domains with branches

We consider the eigenvalue problem for the Laplace operator in a planar domain which can be decomposed into a bounded domain of arbitrary shape and elongated \branches" of variable cross-sectional profiles. When the eigenvalue is smaller than a prescribed threshold, the corresponding eigenfunction decays exponentially along each branch. We prove this behavior for Robin boundary condition and illustrate some related results by numerically computed eigenfunctions.

preprint2011arXiv

Exponential decay of Laplacian eigenfunctions in domains with branches

The behavior of Laplacian eigenfunctions in domains with branches is investigated. If an eigenvalue is below a threshold which is determined by the shape of the branch, the associated eigenfunction is proved to exponentially decay inside the branch. The decay rate is twice the square root of the difference between the threshold and the eigenvalue. The derived exponential estimate is applicable for arbitrary domains in any spatial dimension. Numerical simulations illustrate and further extend the theoretical estimate.

preprint2011arXiv

Trapped modes in finite quantum waveguides

The Laplace operator in infinite quantum waveguides (e.g., a bent strip or a twisted tube) often has a point-like eigenvalue below the essential spectrum that corresponds to a trapped eigenmode of finite L2 norm. We revisit this statement for resonators with long but finite branches that we call "finite waveguides". Although now there is no essential spectrum and all eigenfunctions have finite L2 norm, the trapping can be understood as an exponential decay of the eigenfunction inside the branches. We describe a general variational formalism for detecting trapped modes in such resonators. For finite waveguides with general cylindrical branches, we obtain a sufficient condition which determines the minimal length of branches for getting a trapped eigenmode. Varying the branch lengths may switch certain eigenmodes from non-trapped to trapped states. These concepts are illustrated for several typical waveguides (L-shape, bent strip, crossing of two stripes, etc.). We conclude that the well-established theory of trapping in infinite waveguides may be incomplete and require further development for being applied to microscopic quantum devices.