Pith. sign in

REVIEW 6 minor 73 references

Asymptotic Analysis of the Narrow Escape and Berg-Purcell problems on general three-dimensional domains with reactive boundary patches

T0 review · 0 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper derives three-term asymptotic expansions for the capacitance and the global mean first-reaction time in the small-patch limit, valid on any smooth closed three-dimensional domain with multiple reactive patches.

desk verdict Genuinely new three-term expansions for narrow escape and Berg-Purcell on general smooth 3D domains; the formal machinery holds up, with the Green's function uniformity assumption as the main soft spot. read the letter →

arxiv 2608.08626 v1 pith:LUBTDW6Q submitted 2026-08-09 math.AP

A. E. Lindsay , A.J. Bernoff , D.S. Grebenkov , J.G. Hoskins , M.J. Ward This is my paper · ORCID
classification math.AP MSC 35B2535C2035J0535Q92
keywords narrowescapeproblemBerg-PurcellmatchedasymptoticexpansionssurfaceNeumannGreen'sfunctionreactiveboundarypatchesmeanfirst-reactiontimeprincipalcurvaturesdiffusivecapture
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives three-term asymptotic expansions for two diffusive capture problems on a general smooth closed three-dimensional domain with $N$ small reactive patches: the exterior Berg-Purcell problem, which asks how many diffusing ligands a patchy cell surface captures, and the interior narrow escape problem, which asks how quickly a Brownian particle inside a domain reaches a small reactive window. The central claim is that, to order $\varepsilon^2 \log \varepsilon$, both the capacitance $C_T$ and the global mean first-reaction time are controlled by local patch data—reactive capacitance $C_i$, two monopole coefficients $E_{i\pm}$, mean curvature $H_i$, and curvature difference $Q_i$ at each patch—together with a global surface Neumann Green's matrix that encodes the placement of patches and the shape of the domain. The analysis resolves the weak, path-dependent singularity of the surface Neumann Green's function on a curved boundary, which is the technical step that lets the asymptotics go beyond spheres and constant-curvature surfaces. The formulas reduce to earlier spherical results in the appropriate limits and are checked against a boundary-integral solution on a prolate spheroid, with the predicted $O(\varepsilon^2 \log \varepsilon)$ error.

What carries the argument

The machinery is matched asymptotic expansions in the small patch diameter $\varepsilon$, organized around boundary-fitted tangential-normal coordinates $(t_1,t_2,d)$ aligned with the principal directions at each patch center. In these coordinates the Laplacian acquires explicit curvature terms proportional to mean curvature $H_i$ and curvature difference $Q_i$, and the near-patch problem becomes a half-space Robin problem for each flattened patch $\Gamma_i$. Its solution supplies the reactive capacitance $C_i(b_i)$ and dipole vector $p_i$, while the next-order inner problems define the logarithmic monopole coefficient $E_{i+}$ and the new curvature-anisotropy monopole $E_{i-}$, combined as $E_i=H_i E_{i+}+Q_i E_{i-}$. The outer solution is represented using the surface Neumann Green's function $G_e$ or $G_s$; the load-bearing technical result is the three-term singularity expansion of $G$ near a boundary source, including the path-dependent term $e(x;x_i)$ proportional to $Q$, which makes the matching of inner and outer expansions possible to the retained order.

What would settle it

A high-accuracy boundary-integral computation of $C_T$ for an ellipsoid with one non-circular patch at a point with unequal principal curvatures, repeated for a sequence of $\varepsilon \to 0$, would settle the claim: the residual against Proposition 1 must shrink like $\varepsilon^2 \log \varepsilon$, and the $O(\varepsilon)$ coefficient must contain the $Q E_-$ contribution. If either fails, the expansion is wrong.

Watch

Extended reading notes

Core claim

On the paper's own terms, the main discovery is quantitative and explicit. For the exterior problem, the capacitance satisfies $1/C_T \sim |U_0|\varepsilon^{-1}(1+\varepsilon U_{10}/U_0 \log\varepsilon+\varepsilon U_{11}/U_0)$, with $|U_0|=2/C$, $U_{10}/U_0=-(2C)^{-1}\sum_j H_j C_j^2$, and $U_{11}/U_0=2\pi C^{-1}C^T G_e C+E/C$, where $E=\sum_j(H_j E_{j+}+Q_j E_{j-})$. For the interior problem, the global mean first-reaction time has the same structural expansion with the interior Green's matrix $G_s$ replacing $G_e$ and $U_0=|\Omega|/(2\pi C)$. The new ingredient that distinguishes general surfaces from the sphere is the curvature-anisotropy term $Q_i E_{i-}$: it is nonzero only when a patch is non-circular and sits where the two principal curvatures differ, and it encodes how the orientation of an elliptic patch relative to the principal directions changes capture.

Load-bearing premise

The central premise is that the refined local formula for the surface Neumann Green's function near each patch remains accurate all the way down to the tiny inner region, so that no extra global-shape contribution shows up at the orders kept.

Editorial extensions

If this is right

  • For $N$ identical perfectly absorbing circular patches, the expansions collapse to compact formulas in which only the sum of mean curvatures $H$ and the Green's interaction energy $p_e$ or $p_s$ enter, giving a direct analytic account of how patch arrangement and surface curvature compete.
  • On a prolate spheroid with a single circular patch, the global mean first-reaction time is minimized at the equator, where both the regular part of the Green's function and the mean curvature attain their minima, as verified numerically.
  • When the patch is elliptic and aligned with principal directions, the ratios $E_{\pm}/C^2$ have closed forms, so the mean first-reaction time becomes an explicit function of the semi-axis ratio $(a_1-a_2)/(a_1+a_2)$; non-circular patches on curved surfaces are thus analytically tractable.
  • The same three-term expansion for the global mean first-reaction time yields a three-term expansion for the principal Neumann-Robin eigenvalue $\lambda_0$ of the Laplacian, extending the spherical result to general domains.
  • Optimizing patch positions to maximize capture reduces to minimizing a finite-dimensional discrete energy, whose curvature term has weight depending on patch reactivity through $\gamma(b)$, so the optimal arrangement can shift with receptor properties.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: if the expansion remains valid while the number of patches grows, it suggests a homogenized boundary condition for dense receptor arrays whose coefficients should depend on an average of local mean curvature and on a Green's-function interaction energy; the paper states this as a conjecture rather than deriving it.
  • Editorial inference: because the path-dependent term $e(x;x_i)$ enters the inner problem only through local curvature, surface features smaller than the patch scale that alter principal curvatures should change capture rates at $O(\varepsilon \log \varepsilon)$ even when patch area is fixed, a prediction that could be tested on structured surfaces.
  • Editorial inference: the same inner/outer decomposition should apply to splitting probabilities among competing receptors and to the statistics of the first few binding events, since those quantities depend on the per-patch fluxes that the matched solution already yields; the paper suggests these directions but does not carry them out.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

0 major / 6 minor

Summary. The paper develops matched asymptotic expansions for two diffusive capture problems on smooth closed three-dimensional domains with small, well-separated, partially reactive boundary patches of arbitrary shape: the exterior Berg–Purcell problem and the interior narrow escape problem for the mean first-reaction time. The main results are Proposition 1, a three-term expansion for the inverse capacitance 1/C_T, and Proposition 2, the analogous expansion for the volume-averaged MFRT. Both expansions are expressed in terms of patch capacitances C_i, monopole coefficients E_{i+} and E_{i-}, mean curvature and curvature difference at each patch, and the interior or exterior surface Neumann Green's matrix. The paper also derives a new path-dependent singularity expansion for the surface Neumann Green's function (Appendix B), introduces the new coefficient E_{i-} for non-circular patches (Appendix D), gives explicit formulas for elliptical patches (Appendix E), and validates the asymptotics against a boundary integral solver on prolate spheroids, including a Fibonacci-spiral multi-patch example and an optimization study on nearly spherical domains.

Significance. If correct, the results materially extend the prior theory, which was largely confined to spherical domains, to general smooth geometries, and they identify precisely which local and global geometric data control capture rates at the retained order. The derivation is detailed and internally consistent: the matching through O(1), the switchback logarithms, the new Appendix B singularity analysis, and the E_{i-} lemma all cohere. The paper reduces to known sphere results as special cases and numerically confirms the predicted O(epsilon^2 log epsilon) error for several configurations on prolate spheroids. These strengths give confidence in the central claims. The main results are formal asymptotics rather than theorems with rigorous remainder estimates, which is standard for this literature and is mitigated by the numerical verification. The new E_{i-} term is not directly exercised by the numerical tests, since all boundary-integral validations use circular patches, but the analytical derivation in Appendix D is careful and explicit.

minor comments (6)
  1. [Appendix B, Eqs. (B.26)-(B.29)] The three-term singularity expansion of the surface Neumann Green's function, including the path-dependent term e(x;x_i), is assumed to hold uniformly on the O(epsilon)-scale inner region around each patch. Because this uniformity underlies the matching at O(1) and the U3 monopole condition, I suggest adding one sentence stating this assumption explicitly and, if possible, a numerical check comparing the truncated expansion with a computed Green's function on a non-spherical surface at small distances.
  2. [Lemma 2.3 and Section 4] The new coefficient E_{i-}, which is one of the paper's main novelties, is never exercised by the numerical validation because all boundary-integral tests use circular patches, for which E_{i-}=0. A validation with an elliptical patch at a location with nonzero curvature difference Q_i, or an independent check of the closed-form ratios E_{1+/-}/C_1^2 in Lemma E.1, would further strengthen confidence in this term.
  3. [Proposition 2 and Eq. (3.15)] The pointwise expansion (3.15) contains an undetermined O(epsilon^2 log epsilon) constant U2, as the paper notes. Since the paper's central claim concerns the volume-averaged MFRT (3.16), this does not affect Proposition 2, but I recommend emphasizing in the statement of Proposition 2 that only the global quantity is fully determined to the stated order.
  4. [Section 4.1 and Fig. 4.1] The text describes a boundary integral method, but the caption of Fig. 4.1(b) refers to 'FEM'. Please make the terminology consistent.
  5. [Throughout] There are several typographical errors: 'encloded' in the abstract, 'incling' in the introduction, 'fill fix' in Section 4.3, and 'capacita nce' in the caption of Fig. 4.3. These should be corrected.
  6. [References] Several references are cited as '2026', 'under revision', or 'accepted' (e.g., [14], [16], [29], [42]). If the journal requires final publication data, please update these before publication.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the three-term expansions are derived by matched asymptotics with independently characterized patch-scale inputs.

full rationale

The central derivations in Propositions 1 and 2 (Eqs. 2.46 and 3.16) are carried out in this paper via matched asymptotic expansions: the outer solution is represented in terms of the surface Neumann Green's function, the inner problems on the tangent half-space are solved with stated far-field matching conditions, and the unknown constants U10 and U11 are fixed by solvability conditions. The final capacitance and mean first-reaction time are not defined as the right-hand sides of the formulas; they are extracted from the far-field behavior of the asymptotically constructed solution. The patch-scale quantities C_i, E_{i+}, and the new E_{i-} are auxiliary inputs: C_i and E_{i+} are quoted from the authors' prior work [29,31] (Lemmas 2.1 and 2.2), but they solve well-posed half-space problems that are independent of the present target quantities and can be checked by separate numerical or analytical methods, including the Steklov-based quadratures in Sec. 4.2 and the closed-form ellipse results in Appendix E. The regular parts of the Green's function are computed with the separate boundary-integral solver of [42], and the full Berg-Purcell and narrow-escape solutions in Sec. 4.1 are computed with an independent boundary integral equation for the original problems, so the comparisons in Figs. 4.1-4.3 are genuine external validation rather than fits renamed as predictions. The path-dependent singularity expansion of the Green's function (Appendix B, Eqs. B.26-B.29) is formally derived in this paper; its equivalence to the microlocal result of [47] is noted, and the limitation that no rigorous remainder bound is supplied is a correctness/rigor concern, not a circularity. No uniqueness theorem is imported from the authors to force a choice, no ansatz is smuggled in by citation, and the reductions to known spherical and single-patch results (e.g., Eq. 3.20 and the comparisons with [20,47]) serve as independent consistency checks. The self-citations that occur are normal citations to established auxiliary results and do not carry the load of the central derivation.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The coefficients C_i, E_{i+}, E_{i-}, and the Green's matrix entries are deterministic outputs of auxiliary boundary value problems and numerical solves; none are fitted to the validation data. The heuristic gamma(b) in (4.19b) is used only in a qualitative optimization discussion and does not enter the main theorems. No new physical entities are proposed; the monopole coefficients E_{i\pm} and the path-dependent singularity term e(x;x_i) are mathematical objects derived from the PDE and surface geometry.

assumptions (5)
  • domain assumption The boundary is a smooth, closed, simply connected 3D surface, and the reactive patches are small, non-overlapping, simply connected, smooth-boundary subsets.
    Stated in Section 1 before (1.3); the matched asymptotics and Green's function analysis require smoothness and patch separation |x_i - x_j| = O(1).
  • standard math The surface Neumann Green's functions (2.18a) and (3.2) exist, are unique under their normalization, and have the claimed local singularity behavior.
    Used to represent U_1, U_2, U_3 in Sections 2 and 3; local behavior is derived in Appendix B, while existence relies on standard elliptic theory for smooth surfaces.
  • domain assumption The matched asymptotic expansions converge with the assumed ordering epsilon to 0, b_i = O(1), and with error O(epsilon^2 log epsilon).
    The paper states predicted error rates in Propositions 1 and 2 and validates them numerically; formal Stokes-type matching is assumed rather than proven.
  • standard math Patch-scale results, specifically Lemma 2.1 for C_i(b_i) and Lemma 2.2 for E_{i+}(b_i), are taken from prior works [29,31] with proof summaries in Appendices C and D.
    The leading inner problem solution and its monopole coefficient are inputs to the outer matching; the current paper restates these rather than fully re-deriving them.
  • domain assumption For the interior Green's function, the smooth 1/|Omega| term in (3.2) does not affect the singularity through the retained orders.
    Stated in the last paragraph of Appendix B; true because the term is smooth, but the assertion is not shown in detail.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Asymptotic Analysis of the Narrow Escape and Berg-Purcell problems on general three-dimensional domains with reactive boundary patches." pith.science (2026). https://pith.science/paper/LUBTDW6Q

@misc{pith2026260808626,
  author       = {Pith},
  title        = {Pith review of: Asymptotic Analysis of the Narrow Escape and Berg-Purcell problems on general three-dimensional domains with reactive boundary patches},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LUBTDW6Q}},
  note         = {Machine review of arXiv:2608.08626}
}
read the original abstract

We present an asymptotic analysis of two diffusive capture problems in general smooth closed three-dimensional geometries with multiple small reactive boundary patches of arbitrary shapes. (i) The narrow escape problem seeks to determine the escape rate of Brownian particles from an enclosed region through small boundary windows. (ii) The related Berg-Purcell (or narrow entrance) problem seeks to resolve the capture rate for signaling molecules diffusing outside the cell and entering through localized reactions at membrane-bound receptors. We obtain matched asymptotic solutions of these two problems and thus address the long-standing challenge of describing the role that curvature and local reactivities play in modulating diffusive capture rates. Our explicit expansions quantify local effects on diffusive capture through the sizes, shapes, and reactivities of the patches together with the principal curvatures of the manifold at each patch. In turn, we examine global effects on diffusive capture such as the spatial configuration of patches on the manifold, as encloded by the associated surface Neumann Green's function and its regular part. The accuracy of our asymptotic formulas is validated against a full numerical solution for an ellipsoidal domain. Overall, our results yield new insights on how geometry and stochasticity combine to shape the dynamics of various biological processes.

Figures

Figures reproduced from arXiv: 2608.08626 by the authors.

Figure 1
Figure 1. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 4
Figure 4. [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 4
Figure 4. [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p024_4.png]
Figure 4
Figure 4. Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p027_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

73 extracted references · 71 canonical work pages

  1. [1]

    As a result, we look for a solution to (2.11) in the formw(ξ 1, ξ2, η) =W(λ)

    Level sets of constantλ >0 are ellipsoids that have the limiting behaviorξ 2 1/a2 1 +ξ 2 2/a2 2 →1 andη→0 asλ→0 +, while forλ→ ∞the ellipsoid becomes a large sphereξ 2 1 +ξ 2 2 +η 2 →λof radius √ λ. As a result, we look for a solution to (2.11) in the formw(ξ 1, ξ2, η) =W(λ). From separating variables in the Laplacian written in terms of ellipsoidal coord...

    work page

  2. [2]

    Two Robin patches.As a further test of (1.5), we consider two circular patches of common radiusεand reactivityb= 1

    will generate the largest capaci- tanceC T . Two Robin patches.As a further test of (1.5), we consider two circular patches of common radiusεand reactivityb= 1. We calculate that (4.13a)C= 0.271628, E + = 0.008360, E − = 0. In this case, the formula (1.5) forC T reduces to (4.13b) 1 CT ∼ 2 N Cε 1− C ¯H 2N εlogε+ε 2πCp e(x1,x 2) + ¯HE+ N C , where ¯H=H 1 +...

    work page

  3. [3]

    are (4.16)x k = a q 1−χ 2 k cosϕ k, a q 1−χ 2 k sinϕ k, c χk , k= 0, . . . , N−1. In Fig. 4.3, we plot the rescaled capacita nce (CT U0)/εagainstεfor various values ofN. The solid lines are the asymptotic approximation (2.46) where the entries of the Green’s matrix are constructed using the methods in [42]. The corresponding rescaled capacitance determine...

    work page

  4. [4]

    Using these numerical methods, we validate the main results in Propositions 1 and 2, and demonstrate that these expansions exhibit the predicted error asε→0

    Numerical validation and optimization results.In this section, we de- scribe our methods for the boundary integral solution of the Berg-Purcell (1.4) and narrow escape problems (1.9). Using these numerical methods, we validate the main results in Propositions 1 and 2, and demonstrate that these expansions exhibit the predicted error asε→0. Through some si...

    work page

  5. [5]

    First, we considered the exterior Berg-Purcell problem of receptor-mediated absorption

    Discussion.We have developed a unified asymptotic theory for two canonical diffusive capture problems in general smooth closed three-dimensional domains whose boundary consists of small non-overlapping partially reactive boundary patches of arbitrary shape on an otherwise reflecting surface. First, we considered the exterior Berg-Purcell problem of recept...

    work page 2024

  6. [6]

    +Q(ξ 2 1 −ξ 2 2) +O(ε 3), 1 |x−x i| ∼ 1 ερ 1 + εη 2ρ2 H(ξ2 1 +ξ 2

    work page

  7. [7]

    (A.9) Finally, we observe from (A.9) that|x−x i|=ερ+O(ε 3) andρ= (ξ 2 1 +ξ 2 2)1/2 whenη= 0

    +Q(ξ 2 1 −ξ 2 2) +O(ε). (A.9) Finally, we observe from (A.9) that|x−x i|=ερ+O(ε 3) andρ= (ξ 2 1 +ξ 2 2)1/2 whenη= 0. As a result, if the orthogonal projection of the Robin patch∂Ω ε i onto the tangent is a disk of radiusεa i, then the rescaled surface patch Γ i is described by ξ2 1 +ξ 2 2 1/2 ≤a i +O(ε 2). This smallO(ε 2) difference between the the ortho...

    work page

  8. [8]

    A. J. Bernoff and A. E. Lindsay. Numerical approximation of diffusive capture rates by planar and spherical surfaces with absorbing pores.SIAM J. Appl. Math., 78(1):266–290, 2018

    work page 2018

Show all 73 references
  1. [9]

    Alonso, R

    A. Alonso, R. G. Endres, and J. B. Kirkegaard. Local clustering and global spreading of receptors for optimal spatial gradient sensing.Phys. Rev. Lett., 134:158401, 2025

    2025

  2. [10]

    Arendt and ter Elst A

    W. Arendt and ter Elst A. F. M. The Dirichlet-to-Neumann operator on exterior domains.Poten. Anal., 43:313–340, 2015

    2015

  3. [11]

    Askham, L

    T. Askham, L. Greengard, J. Hoskins, L. Lu, M. O’Neil, M. Rachh, F. Vico, and V. Rokhlin. fmm3dbie: FMM-accelerated boundary integral equation solvers

  4. [12]

    J. W. S. Baron Rayleigh.The Theory of Sound, volume 2. Dover, New York, 2nd edition, 1945

    1945

  5. [13]

    B´ enichou and R

    O. B´ enichou and R. Voituriez. From first-passage times of random walks in confinement to geometry-controlled kinetics.Phys. Rep., 539(4):225–284, 2014

    2014

  6. [14]

    H. C. Berg and E. M. Purcell. Physics of chemoreception.Biophys J., 20(2):193– 219, 1977

    1977

  7. [15]

    A. J. Bernoff, A. Jilkine, A. Navarro Hern´ andez, and A. E. Lindsay. Single- cell directional sensing from just a few receptor binding events.Biophys. J., 122:3108–3116, 2023

    2023

  8. [16]

    Chakraborty, J

    S. Chakraborty, J. G. Hoskins, and A. E. Lindsay. Fast integral methods for the Neumann Green’s function: applications to capture and signaling problems in two dimensions, 2025

    2025

  9. [17]

    A. J. Bernoff and A. E. Lindsay. Kinetic Monte Carlo methods for three- dimensional diffusive capture problems in exterior domains.Royal Soc. Open Sci., 12(2):241033, 2025. The narrow escape and Berg-Purcell problems.47

    2025

  10. [18]

    A. J. Bernoff, A. E. Lindsay, and D. D. Schmidt. Boundary homogenization and capture time distributions of semipermeable membranes with periodic patterns of reactive sites.Multiscale Model. & Simul., 16(3):1411–1447, 2018

    2018

  11. [19]

    P. C. Bressloff. Asymptotic analysis of particle cluster formation in the presence of anchoring sites.Eur. Phys. J. E, 47(5):30, 2024

    2024

  12. [20]

    P. C. Bressloff and J. M. Newby. Stochastic models of intracellular transport. Rev. Mod. Phys., 85:135–196, 2013

    2013

  13. [21]

    Bundrock, A

    L. Bundrock, A. Girouard, D. S. Grebenkov, M. Levitin, and I. Polterovich. The exterior Steklov problem for Euclidean domains.J. Spect. Th., accepted; preprint is available arXiv:2511.09490, 2025

    2025 arXiv

  14. [22]

    The narrow escape problem in arbitrary dimension, 2026

    Louis Carillo, Tony Leli` evre, Thomas Normand, Gabriel Stoltz, and Urbain Vaes. The narrow escape problem in arbitrary dimension, 2026

    2026

  15. [23]

    Cavanagh, D

    H. Cavanagh, D. Kempe, J. K. Mazalo, M. Biro, and R. G. Endres. T cell mor- phodynamics reveal periodic shape oscillations in three-dimensional migration. J. Royal Soc. Interface, 19(190):20220081, 2022

    2022

  16. [24]

    D. S. Grebenkov. Spectral theory of imperfect diffusion-controlled reactions on heterogeneous catalytic surfaces.J. Chem. Phys., 151:104108, 2019

    2019

  17. [25]

    Chen and A

    X. Chen and A. Friedman. Asymptotic analysis for the narrow escape problem. SIAM J. Math. Anal., 43(6):2542–2563, 2011

    2011

  18. [26]

    Cheviakov and M

    A. Cheviakov and M. Ward. Optimization of trap locations for narrow capture problems. In D. S. Grebenkov, R. Metzler, and G. Oshanin (Eds.), editors,Target Search Problems, pages 225–246. Springer Nature Switzerland, Cham, 2024

    2024

  19. [27]

    Cheviakov and M

    A. Cheviakov and M. J. Ward. Optimizing the principal eigenvalue of the Lapla- cian in a sphere with interior traps.Math. Comp. Model., 53(7-8):1394–1409, 2011

    2011

  20. [28]

    A. F. Cheviakov, M. J. Ward, and R. Straube. An asymptotic analysis of the mean first passage time for narrow escape problems: Part II: The sphere.Mul- tiscale Model. Simul., 8(3):836–870, 2010

    2010

  21. [29]

    Dobramysl and D

    U. Dobramysl and D. Holcman. Mixed analytical-stochastic simulation method for the recovery of a Brownian gradient source from probability fluxes to small windows.J. Comput. Phys., 355:22–36, 2018

    2018

  22. [30]

    Gomez and A

    D. Gomez and A. F. Cheviakov. Asymptotic analysis of narrow escape problems in nonspherical three-dimensional domains.Phys. Rev. E, 91(1):012137, 2015

    2015

  23. [31]

    Gonz´ alez

    ´A. Gonz´ alez. Measurement of areas on a sphere using Fibonacci and latitude– longitude lattices.Math. Geosci., 42(1):49–64, 2009

    2009

  24. [32]

    Greengard, M

    L. Greengard, M. O’Neil, M. Rachh, and F. Vico. Fast multipole methods for the evaluation of layer potentials with locally-corrected quadratures.J. Comput. Phys. X, 10:100092, 2021

    2021

  25. [33]

    D. S. Grebenkov. Paradigm shift in diffusion-mediated surface phenomena.Phys. Rev. Lett., 125:078102, 2020

    2020

  26. [34]

    D. S. Grebenkov. Diffusion-controlled reactions: An overview.Molecules, 28(22):7570, 2023

    2023

  27. [35]

    D. S. Grebenkov. Spectral properties of the Dirichlet-to-Neumann operator for spheroids.Phys. Rev. E, 109:055306, 2024

    2024

  28. [36]

    D. S. Grebenkov and R. Maurette. Reactive capacitance of flat patches of arbi- trary shape.Phys. Rev. E., 113:034112, 2026

    2026

  29. [37]

    D. S. Grebenkov and M. J. Ward. The asymptotic analysis of some PDE and Steklov eigenvalue problems with partially reactive patches in 3-D.under revi- sion, SIAM Review, ArXiv: 2509.17394 [math.AP], 2026

    2026 arXiv

  30. [38]

    D. S. Grebenkov and M. J. Ward. Competition of small targets in planar domains: 48A. E. Lindsay, A.J. Bernoff, D.S. Grebenkov, J.G. Hoskins, M.J. Ward from Dirichlet to Robin and Steklov boundary condition.Eur. J. Appl. Math., 5:1–48, 2026

    2026

  31. [39]

    D. S. Grebenkov and M. J. Ward. The effective reactivity for capturing Brownian motion by partially reactive patches on a spherical surface.Multiscale Model. Simul., 24(2):660–692, 2026

    2026

  32. [40]

    S. D. Lawley, A. E. Lindsay, and C. E. Miles. Receptor organization determines the limits of single-cell source location detection.Phys. Rev. Lett., 125(1):018102, 2020

    2020

  33. [41]

    Handy and S

    G. Handy and S. D. Lawley. Revising Berg-Purcell for finite receptor kinetics. Biophys. J., 120(11):2237–2248, 2021

    2021

  34. [42]

    Holcman and Z

    D. Holcman and Z. Schuss. The narrow escape problem.SIAM Rev., 56(2):213– 257, 2014

    2014

  35. [43]

    Kaye and L

    J. Kaye and L. Greengard. A fast solver for the narrow capture and narrow escape problems in the sphere.J. Comput. Phys. X, 5:100047, 2020

    2020

  36. [44]

    H. Kersten. Grenz- und sprungrelationen f¨ ur potentiale mit quadrat- summier- barer f1¨ achenbelegung.Resultate der Mathematik, 3:17–24, 1980

    1980

  37. [45]

    Kress.Linear Integral Equations

    R. Kress.Linear Integral Equations. Applied Mathematical Sciences, Springer, 2014

    2014

  38. [46]

    Lagache and D

    T. Lagache and D. Holcman. Extended narrow escape with many windows for analyzing viral entry into the cell nucleus.J. Stat. Phys., 166(2):244–266, 2017

    2017

  39. [47]

    In our analysis below, for ease of notation we omit the subscriptefor this Green’s function

    from microlocal analysis; here, we present an alternative approach based on a local coordinate system defined in Appendix A which is suited for the high-order asymptotic analysis of (1.4) and (1.9). In our analysis below, for ease of notation we omit the subscriptefor this Gre...

  40. [48]

    Lagerstrom and D

    P. Lagerstrom and D. Reinelt. Note on logarithmic switchback terms in regular and singular perturbation problems.SIAM J. Appl. Math., 44(3):451–462, 1984

    1984

  41. [49]

    A. E. Lindsay and A. J. Bernoff. Numerical approximation of diffusive capture rates by planar and spherical surfaces with absorbing pores.SIAM J. Appl. Math., 78(1):266–290, 2018

    2018

  42. [50]

    A. E. Lindsay, A. J. Bernoff, T. Goodwill, and J. G. Hoskins. The three dimen- sional Neumann Green’s function for general surfaces: singular asymptotics and boundary integral methods, 2026

    2026

  43. [51]

    A. E. Lindsay, A. J. Bernoff, and A. Navarro Hern´ andez. Short-time diffusive fluxes over membrane receptors yields the direction of a signalling source.Royal Soc. Open Sci., 10(4):221619, 2023

    2023

  44. [52]

    A. E. Lindsay, A. J. Bernoff, and M. J. Ward. First passage statistics for the capture of a Brownian particle by a structured spherical target with multiple surface traps.Multiscale Model. Simul., 15(1):74–109, 2017

    2017

  45. [53]

    Metzler, G

    R. Metzler, G. Oshanin, and S. Redner, editors.First-Passage Phenomena and Their Applications. World Scientific, 2014

    2014

  46. [54]

    I. M. Nemenman and A. S. Silbergleit. Explicit Green’s function of a bound- ary value problem for a sphere and trapped flux analysis in Gravity Probe B experiment.J. Appl. Phys., 86(1):614–624, 1999

    1999

  47. [55]

    Nursultanov, J

    M. Nursultanov, J. C. Tzou, and L. Tzou. On the mean first arrival time of Brownian particles on Riemannian manifolds.J. Math. Pures et Appl., 150:202– 240, 2021

    2021

  48. [56]

    P. O. Persson and G. Strang. A simple mesh generator in matlab.SIAM Rev., 46(2):329–345, 2004

    2004

  49. [57]

    C. E. Plunkett and S. D. Lawley. Boundary homogenization for partially reactive patches.Multi. Model. & Simul., 22(2):784–810, 2024

    2024

  50. [58]

    I. Y. Popov. Extension theory and localization of resonances for domains of trap The narrow escape and Berg-Purcell problems.49 type.Mathematics of the USSR-Sbornik, 71(1):209, feb 1992

    1992

  51. [59]

    Redner.A Guide to First-Passage Processes

    S. Redner.A Guide to First-Passage Processes. Cambridge University Press, 2001

    2001

  52. [60]

    Reingruber and D

    J. Reingruber and D. Holcman. Gated narrow escape time for molecular signaling. Phys. Rev. Lett., 103:148102, 2009

    2009

  53. [61]

    Schamberger, A

    B. Schamberger, A. Roschger, R. Ziege, K. Anselme, M. B. Amar, M. Bykowski, A. P. G. Castro, A. Cipitria, R. Coles, R. Dimova, M. Eder, S. Ehrig, L. M. Es- cudero, M. E. Evans, P. R. Fernandes, P. Fratzl, L. Geris, N. Gierlinger, E. Han- nezo, A. Igliˇ c, J. J. K. Kirkensgaard...

    2023

  54. [62]

    Z. Schuss. The narrow escape problem: A short review of recent results.J. Sci. Comput., 53(1):194–210, 2012

    2012

  55. [63]

    Schuss.Brownian Dynamics at Boundaries and Interfaces

    Z. Schuss.Brownian Dynamics at Boundaries and Interfaces. Springer, 2013

    2013

  56. [64]

    Schuss, A

    Z. Schuss, A. Singer, and D. Holcman. The narrow escape problem for diffusion in cellular microdomains.Proc. Nat. Acad. Sci. USA, 104(41):16098–16103, 2007

    2007

  57. [65]

    Silbergleit, I

    A. Silbergleit, I. Mandel, and I. Nemenman. Potential and field singularity at a surface point charge.J. Math. Phys., 44(10):4460–4466, 2003

    2003

  58. [66]

    Singer, Z

    A. Singer, Z. Schuss, D. Holcman, and R. S. Eisenberg. Narrow escape, Part I. J. Stat. Phys., 122(3):437–463, 2006

    2006

  59. [67]

    E. M. Stein and R. Shakarchi.Complex Analysis, volume 2. Princeton University Press, Princeton, NJ, 2003

    2003

  60. [68]

    T. L. Stepien, C. Zmurchok, J. B. Hengenius, R. M. Caja Rivera, M. R. D’Orsogna, and A. E. Lindsay. Moth mating: Modeling female pheromone call- ing and male navigational strategies to optimize reproductive success.Appl. Sci., 10(18), 2020

    2020

  61. [69]

    Strieder

    W. Strieder. Diffusion-controlled reaction on an elliptic site.J. Chem. Phys., 130(17):176103, 2009

    2009

  62. [70]

    Swinbank and R

    R. Swinbank and R. J. Purser. Fibonacci grids: A novel approach to global modelling.Quart. J. Royal Meteorol. Soc., 132(619):1769–1793, 2006

    2006

  63. [71]

    Vioreanu and V

    B. Vioreanu and V. Rokhlin. Spectra of multiplication operators as a numerical tool.SIAM J. Scient. Comput., 36(1):A267–A288, 2014

    2014

  64. [72]

    F. Wei, D. Yang, R. Straube, and J. Shuai. Brownian diffusion of ion channels in different membrane patch geometries.Phys. Rev. E, 83:021919, 2011

    2011

  65. [73]

    Xiao and Z

    H. Xiao and Z. Gimbutas. A numerical algorithm for the construction of efficient quadrature rules in two and higher dimensions.Comput. Math. Appl., 59:663– 676, 2010

    2010

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.