Researcher profile

Denis S. Grebenkov

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

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

14 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.

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.

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.

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.