Minimizing capture time with many small traps in heterogeneous media
Abstract
We study the problem of optimally placing a large number of small absorbing traps to minimize the mean first-passage time (MFPT) of particles diffusing in a heterogeneous medium with space-dependent diffusivity D(x) and prescribed initial particle distribution ω(x). In two and three dimensions we identify two distinct regimes, depending on how strongly the trap ensemble depletes the particles. In the weak trapping regime (fewer traps) the optimal trap density ρ is the arithmetic average ρ= 1/2(μ+ω) of the normalized inverse diffusivity μ∝ 1/D and the initial particle distribution ω. In the strong trapping regime (more traps), it is proportional to their geometric average, ρ∝ √(μω). We revisit the classical Green's function approach and show that it applies to weak trapping only; in 2D this requires the trap size to be exponentially small in the number of traps, so that most applications of interest fall in the strong trapping regime instead. We therefore develop a homogenization approach that captures both regimes. A special choice of initial particle distribution is ω∝ 1/D, which corresponds to the equilibrium particle distribution in the absence of traps according to the Itô interpretation. In this case we show that the optimal trap distribution is the stationary particle distribution itself, regardless of the trapping strength, and this makes the MFPT constant throughout the domain; the same trap distribution also minimizes the worst-case capture time. Finally, in one dimension we find instead ρ∝(μω) 1/3, and a similar analysis covers space-dependent drift and thin domains of variable cross-section. Direct numerical optimization confirms the analytical results.
Related
- Microlocal defect functionals in VMO: Geometric localisation and applications to highly heterogeneous media
- An accurate parallel finite difference method for the variable-order fractional Laplacian in heterogeneous media
- GeONet: Geometry-Encoded Operator Network for Predicting Flow in Heterogeneous Porous Media
- On a conjecture of Almgren: area-minimizing submanifolds with fractal singular sets
- On Simon's area-minimizing hypersurfaces with fractal singular sets
BibTeX
@misc{grebenkov2026minimizing,
title = {{Minimizing capture time with many small traps in heterogeneous media}},
author = {Denis S. Grebenkov and Theodore Kolokolnikov},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.39990},
eprint = {2609.39990},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.39990},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Minimizing capture time with many small traps in heterogeneous media
AU - Denis S. Grebenkov
AU - Theodore Kolokolnikov
PY - 2026
DA - 2026/09/30
PB - arXiv
DO - 10.48550/arxiv.2609.39990
UR - https://arxiv.org/abs/2609.39990
AB - We study the problem of optimally placing a large number of small absorbing traps to minimize the mean first-passage time (MFPT) of particles diffusing in a heterogeneous medium with space-dependent diffusivity D(x) and prescribed initial particle distribution ω(x). In two and three dimensions we identify two distinct regimes, depending on how strongly the trap ensemble depletes the particles. In the weak trapping regime (fewer traps) the optimal trap density ρ is the arithmetic average ρ= 1/2(μ+ω) of the normalized inverse diffusivity μ∝ 1/D and the initial particle distribution ω. In the strong trapping regime (more traps), it is proportional to their geometric average, ρ∝ √(μω). We revisit the classical Green's function approach and show that it applies to weak trapping only; in 2D this requires the trap size to be exponentially small in the number of traps, so that most applications of interest fall in the strong trapping regime instead. We therefore develop a homogenization approach that captures both regimes. A special choice of initial particle distribution is ω∝ 1/D, which corresponds to the equilibrium particle distribution in the absence of traps according to the Itô interpretation. In this case we show that the optimal trap distribution is the stationary particle distribution itself, regardless of the trapping strength, and this makes the MFPT constant throughout the domain; the same trap distribution also minimizes the worst-case capture time. Finally, in one dimension we find instead ρ∝(μω) ^1/3, and a similar analysis covers space-dependent drift and thin domains of variable cross-section. Direct numerical optimization confirms the analytical results.
ER -