alphaFeedback

WorksMathematics

Narrow escape problems beyond spherical holes: slits and tubes

AI contributions
  • Claude Opus 5.5
View details

Abstract

The narrow escape problem consists in characterizing the exit event of a Brownian motion trapped inside a domain with reflecting boundary and small absorbing holes. Using recent advances in the quasi-stationary distribution approach, we first give an explicit set of hypotheses that allows one to derive the asymptotics of the exit time and of the exit hole distribution. We then apply it in dimension three to several hole shapes: spherical holes, thin spheroidal holes shrinking to curves on the boundary (slits), and holes located at the end of thin tubes attached to the domain. We illustrate the accuracy of the results with numerical simulations, and extend them to higher dimensions.

Introduction

The narrow escape problem was initially introduced in biology, in the context of the diffusion of molecules in a cell [27], in particular to model the diffusion of neurotransmitters in a synapse. This physical phenomenon can be modeled by a reflected Brownian motion constrained to a domain, with a small absorbing part on its boundary, modeling the exit of the molecule, or its reaction with a receptor. Under the assumption that the absorbing part is small, the goal is to compute the mean exit time of the molecule and, when there are several absorbing parts, the probability of exiting through each of them. This problem has attracted interest both in the physics community [57, 24, 49] [5] and in the mathematical community [39, 38, 30] [1, 15]. In this work, we extend the contributions of [38] [13, 14] to study different shapes of the absorbing part, within a similar framework: the quasi-stationary approach. This introduction is organized as follows. First, we define the basic setting of the narrow escape problem and the quasi-stationary distribution in § 1.1; then we describe the various hole shapes studied in this work and their connection to the literature in § 1.2; finally, the contributions and the outline of the paper are given in § 1.3.

Introduction to the narrow escape problem and the quasi-stationary distribution

In this section, we briefly introduce the mathematical fundamentals of the narrow escape problem and justify the quasi-stationary approach we use in the rest of the paper. Let Ω⊂R3\Omega\subset\mathbb{R}^{3} be a smooth bounded domain with boundary ∂Ω\partial\Omega. Denote by (Xt)t≥0(X_t)_{t\geq0} the stochastic process modeling the diffusion of a particle trapped inside Ω\Omega with reflecting boundary conditions on ∂Ω\partial\Omega. The goal is to study the exit of the particle through a small absorbing part ΓDε⊂∂Ω\Gamma_D^\varepsilon\subset\partial\Omega, where ε>0\varepsilon> 0 controls the size of ΓDε\Gamma_D^\varepsilon. Explicitly, this exit can be characterized by a pair of random variables: the exit time

τε=inf⁡{t>0∣Xt∈ΓDε‾},(1)\tau_\varepsilon= \inf\{t > 0 \mid X_t \in\overline{\Gamma_D^\varepsilon}\}, \tag*{(1)}

and the exit position

Xτε∈ΓDε‾.(2)X_{\tau_\varepsilon} \in\overline{\Gamma_D^\varepsilon}. \tag*{(2)}

In general, both of these random variables depend on the initial condition X0X_0 of the process (Xt)t≥0(X_t)_{t\geq0}. To overcome this difficulty, in this work we consider the initial condition to be distributed according to the quasi-stationary distribution νε\nu_{\varepsilon}. The quasi-stationary distribution νε\nu_{\varepsilon} is a probability measure on Ω\Omega that is stationary under the dynamics conditioned on the event that the particle has not escaped yet [16, 17]. It has been shown in [13, Theorem 1.1, Remark 1.2] that νε\nu_{\varepsilon} exists in a wide class of narrow escape problems and that, for any initial condition X0∈ΩX_0 \in\Omega,

Law⁡(Xt∣t<τε)⟶t→+∞νε.\operatorname{Law}(X_t \mid t < \tau_{\varepsilon}) \underset{t \to+\infty}{\longrightarrow} \nu_{\varepsilon}.

Hence, at large times, the law of the process (Xt)t≥0(X_t)_{t \ge0}, conditioned on the fact that it has not escaped yet, tends to the quasi-stationary distribution (see Figure 1 for an illustration). This relaxation towards νε\nu_{\varepsilon} occurs on a time scale that stays bounded as ε→0\varepsilon\to0, and thus is negligible compared to the exit time τε\tau_{\varepsilon}, which diverges in this limit. Additionally, starting from νε\nu_{\varepsilon}, [13, Proposition 1.3] shows that both τε\tau_{\varepsilon} and XτεX_{\tau_{\varepsilon}} can be derived from νε\nu_{\varepsilon} (see their explicit relation in § 2.3). Consequently, the goal of this paper is to study the quasi-stationary distribution νε\nu_{\varepsilon}, for different shapes of the absorbing part ΓDε\Gamma_{\mathrm{D}}^{\varepsilon}, and to derive the asymptotic behavior of the exit time τε\tau_{\varepsilon} and exit position XτεX_{\tau_{\varepsilon}} in the limit ε→0\varepsilon\to0.

Histogram of the exit time of 100 000 particles

Figure 1. Histogram of the exit time of 100 000 particles in a unit ball of dimension 3 with 2 absorbing spherical caps of angle 0.2 radians on its boundary. Notice the exponential tail of the distribution, drawn in black, typical of the exit time under the quasi-stationary distribution. The vertical red line represents the time from which one can heuristically consider the particles to be distributed according to the quasi-stationary distribution. In this case, less than 20% of the particles have escaped before this time. This proportion shrinks to 0 in the limit ε→0\varepsilon\to0 (see [13, Section 4.2] for an illustration of this phenomenon).

Hole shapes and bibliography review

Let us now describe the different hole shapes studied in this work, which model various physical phenomena; see the global illustration in Figure 2.

Illustration of the different hole shapes considered in this work

Figure 2. Illustration of the different hole shapes considered in this work. The gray plane represents the boundary of the domain ∂Ω\partial\Omega, while the colored areas are the holes ΓDε\Gamma_{\mathrm{D}}^{\varepsilon}. The parameter ε\varepsilon illustrates how the holes shrink.

The narrow escape problem is the standard setting introduced by [27], where the holes are geodesic disks on the boundary of the domain that shrink to a point as ε→0\varepsilon\to0. It is represented in blue in Figure 2 (B). Some notable contributions in the literature on this version of the problem are [64, 56, 57, 45] [5, 1, 15] which all focus on the mean exit time or on the exit time from a fixed starting point. The quasi-stationary distribution approach has been used in [14] to study the exit time and exit position in dimension 2. Notice two other contributions [38] [13], which use the quasi-stationary distribution but for slightly modified holes that are half-spheres inside the domain; see Figure 2 (A). One goal of this paper is to understand the difference between the two settings.

The narrow slit problem is a variant of the narrow escape problem where one considers a hole that shrinks to a curve of finite length on the boundary of the domain instead of a point. In practice, this is built from a prolate spheroid that shrinks to a curve segment (see Figure 2 (C)). Notice that this problem is fundamentally different, as the hole shrinks to an object of dimension 1. Heuristically, one could expect the mean exit time to diverge logarithmically, as in the two-dimensional narrow escape problem, pointing to the importance of the codimension of the hole instead of the dimension of the domain.

The geometry behind the narrow slit problem is not new and has been studied in a wide range of applications, from electromagnetism [9] and mechanics [10] to hydrodynamics [35] and acoustics [43]. Yet, to the best of our knowledge, this regime has not been studied in the context of the narrow escape problem. However, using the results on ellipsoidal holes from [58, 45], we can gain insight into the behavior of the narrow slit problem. Formally, taking the ellipticity of the ellipsoid to scale as ε−1\varepsilon^{-1} in [45 Theorem 1.3] leads to a logarithmic divergence of the mean exit time, which confirms the heuristic prediction above. Finally, note that this problem can be tackled with capacity methods to obtain the scaling in ε\varepsilon of the mean exit time (see [21]).

The narrow transition problem is the variant illustrated in purple in Figure 2 (D), where one considers a tube of fixed length with a shrinking cross-section. The exit hole is at the end of the tube. The narrow escape community has considered this problem since [28]; see a list of recent contributions on this topic in [39, 31, 30, 47] [6]. Let us also point out the strong connection between this problem and Helmholtz resonators in acoustics/electrostatics [32, 29, 33] [4].

Another interpretation of this problem arises from statistical physics, in the context of transitions from one state to another (see [65] for a discussion of several problems involving Brownian motions trapped in channels, and [22] for entropy-driven phenomena). Indeed, if at the end of the tube there is another domain Ω′\Omega', then the escape time through the tube is also the transition time from Ω\Omega to Ω′\Omega'. Transition times are fundamental objects in the study of metastable systems [36, 25, 62]. Usually, the metastability is induced by an energetic barrier between Ω\Omega and Ω′\Omega', hence the transition time is trivially linked by a factor 2 to the escape time from Ω\Omega (see [37 Remark 1]). In this work, where Ω\Omega and Ω′\Omega' are linked by a narrow tube, we see that this is not the case, and the transition time can be much larger than the escape time.

The analogy between the narrow transition problem and the energetic barrier framework can be pushed further. The transition time for the latter is given by the Eyring–Kramers law [36, 25, 37],

E[τ]∼2π∣λ1(z)∣∣det⁡∇2V(z)∣det⁡∇2V(x)exp⁡(V(z)−V(x)h),\mathbb{E}[\tau] \sim\frac{2\pi}{\lvert\lambda_1(z)\rvert}\sqrt{\frac{\lvert\det\nabla^2 V(z)\rvert}{\det\nabla^2 V(x)}}\exp\left(\frac{V(z)-V(x)}{h}\right),

where hh is proportional to the temperature and plays the role of the small parameter, VV is the energetic potential, λ1\lambda_1 is the negative eigenvalue of the Hessian of VV at the saddle point zz between the two wells, and xx is the local minimum of VV in the well where the process starts. Formally, the ratio of the two Hessians can be seen as a measure of the volume of the well over the volume of the transition region ∣Ω∣ε2L\frac{\lvert\Omega\rvert}{\varepsilon^2 L}, where LL is the length of the tube and ε2\varepsilon^2 its cross-sectional area. The factor ∣λ1∣−1\lvert\lambda_1\rvert^{-1} is a measure of the rate to traverse the transition region, here L2L^2. Taking the exponent to 0 as there are no energetic barriers in the narrow transition problem, the formula for the mean exit time can be understood as

E[τ]∼L∣Ω∣ε2.\mathbb{E}[\tau] \sim L\frac{\lvert\Omega\rvert}{\varepsilon^2}.

In § 1.3, we see that these heuristics lead to the correct scaling behavior for the mean exit time in the narrow transition problem, making a link between metastability arising from an energetic constraint and from an entropic constraint.

Contribution and outline of the paper

The contribution of this work is to derive the exit time distribution and the exit hole distribution for the narrow escape problem (in two variants, A and B), the slit problem (C), and the transition problem (D), within a single quasi-stationary distribution framework. Let us emphasize that the novelty of this work is to demonstrate that the framework introduced in [38] can be applied to a wide range of hole shapes, and to derive asymptotically the exit time distribution and the exit hole distribution from the quasi-stationary distribution, in the limit ε→0\varepsilon\to0. To do this, new quasi-modes have been built to approximate the quasi-stationary distribution for the different hole shapes. This technique is explained in § 2.3. The results of this paper are gathered in Table 1, which summarizes the different hole shapes and the asymptotics of their mean exit time and the exit hole distribution in a reference setting (one hole for the mean exit time, two holes for the exit hole) in the limit ε→0\varepsilon\to0.

Hole typeHole descriptionScaling of the surface areaMean exit timeExit holeReference
ASpherical holes inside the domainε2\varepsilon^2∣Ω∣2πε\frac{|\Omega|}{2\pi\varepsilon}13\frac{1}{3}Proposition 3.1
BSpherical holes on the boundaryε2\varepsilon^2∣Ω∣4ε\frac{|\Omega|}{4\varepsilon}13\frac{1}{3}Proposition 3.3
CProlate spheroid around a curved segment of length 2ℓ2\ell on ∂Ω\partial\Omega, half-width ε\varepsilonεℓ\varepsilon\ell∣Ω∣2πℓ∣log⁡(ε)∣\frac{|\Omega|}{2\pi\ell|\log(\varepsilon)|}12\frac{1}{2}Proposition 4.2
DEnd of a tube of fixed length LL with shrinking cross-section of area πε2\pi\varepsilon^2ε2\varepsilon^2∣Ω∣Lπε2\frac{|\Omega|L}{\pi\varepsilon^2}15\frac{1}{5}Proposition 5.3

Table 1. Summary of the different hole shapes considered in this work and their properties.

For the sake of comparison, we also include a column stating how the surface area of each hole scales with ε\varepsilon. Notice that the results in Table 1 are consistent with the heuristics on the mean exit time developed in § 1.2. Furthermore, we make explicit the difference between the spherical holes on the boundary studied in [14] and the spherical holes studied in [38] [13], which are half-spheres inside the domain. The scaling of the mean exit time differs by a factor π/2\pi/2 between the two settings, giving a rigorous explanation of [13, Remark 1.5].

The notation A, B, C and D refers to the illustrations in Figure 2. See (Type A), (Type B), (Type C) and (Type D) in § 2.2 for the precise definitions of the holes. The mean exit time is Eνε[τε]\mathbb{E}_{\nu_\varepsilon}[\tau_\varepsilon], the expectation of (1.1) under the quasi-stationary distribution νε\nu_\varepsilon. The exit hole column refers to the probability of exiting through a hole of parameter ε\varepsilon in a system with another hole of parameter 2ε2\varepsilon (with the same ℓ\ell for C and the same LL for D). The reference column gives the proposition in the paper where the result is stated and proved. Each result is given as its asymptotic equivalent in the limit ε→0\varepsilon\to0.

The results presented in Table 1 are believed to be derivable in higher dimensions, at the cost of lengthy computations. This is one goal of the numerical experiments in § 6, which aim to illustrate the results of this work and to explore the higher-dimensional case.

The paper is organized as follows. In § 2, we introduce the main notation and explain the quasi-stationary distribution framework. Then, in § 3, we present our results on the spherical holes A and B. In § 4, we construct a quasi-mode for the narrow slit problem and derive the asymptotics of the exit time and exit hole distributions. The same is done for the narrow transition problem in § 5. Finally, we present numerical experiments in § 6 to illustrate the results of this work and explore the higher-dimensional case. Additionally, in Section A, we discuss the Poincaré–Wirtinger inequality uniform in ε\varepsilon in the several geometries of this work. We then provide the proof of Theorem 2.2 in Section B. Section C gathers the fundamental properties of the Green function considered in this work, and Section D is a technical section on the computation of several integrals.

Notation and detailed method

In this section we introduce the mathematical notation, first in general in § 2.1, then for the different hole shapes in § 2.2. Finally, in § 2.3, we explain how to use the quasi-stationary distribution to derive the exit time and exit hole distributions, using a quasi-mode.

Notation

2.1. Notation. In this section, we define the notation and the mathematical concepts used throughout this work. We first introduce geometric considerations, before introducing the operator associated with the process (Xt)t≥0(X_t)_{t\geq0} defined in (3). Note that in this paper, four geometric settings are considered. For simplicity, we use the same notation in all four settings, with the understanding that the underlying geometric setting is clear from the context.

We recall that Ω\Omega is a smooth bounded domain of R3\mathbb{R}^3 with boundary ∂Ω\partial\Omega. Let (Xt)t≥0(X_t)_{t\geq0} be the following stochastic process evolving in Ω‾\overline{\Omega}:

dXt=2 dBt−n(Xt) dLt,(3)\mathrm{d}X_t = \sqrt{2}\,\mathrm{d}B_t - n(X_t)\,\mathrm{d}L_t, \tag*{(3)}

where (Bt)t≥0(B_t)_{t\geq0} is a standard Brownian motion in R3\mathbb{R}^3, (Lt)t≥0(L_t)_{t\geq0} is the local time of (Xt)t≥0(X_t)_{t\geq0} on ∂Ω\partial\Omega (see [60, 63]), and nn is the outward unit normal vector to ∂Ω\partial\Omega. This process is a reflected Brownian motion in Ω\Omega.

For N≥1N\geq1 and ε>0\varepsilon>0, we consider a family of NN (relatively) open connected sets Γεk⊂Ω‾\Gamma_\varepsilon^k\subset\overline{\Omega}, k∈{1,…,N}k\in\{1,\ldots,N\}, called the holes, with a small parameter rε(k)>0r_\varepsilon^{(k)}>0 that controls their size (for the (Type C) holes, wε(k)w_\varepsilon^{(k)} plays a similar role). We write ΓDε=⋃k=1NΓεk\Gamma_\mathrm{D}^\varepsilon=\bigcup_{k=1}^N\Gamma_\varepsilon^k for their union. The small hole asymptotics is defined by the limit ε→0\varepsilon\to0, where rε(k)→0r_\varepsilon^{(k)}\to0 and ∣Γεk∣→0|\Gamma_\varepsilon^k|\to0. In this work, we emphasize that we only consider the limit where the holes shrink but do not move. Hence, we always assume that the holes remain separated by a distance bounded below uniformly in ε\varepsilon.

The generator of (Xt)t≥0(X_t)_{t\geq0}. The generator of the process (Xt)t≥0(X_t)_{t\geq0}, defined in (3), killed on ΓDε\Gamma_\mathrm{D}^\varepsilon is the Laplacian Δ\Delta on Ωε\Omega_\varepsilon with Dirichlet boundary conditions on ΓDε\Gamma_\mathrm{D}^\varepsilon and Neumann boundary conditions on ΓNε\Gamma_\mathrm{N}^\varepsilon. In the following, we work with the nonnegative operator Lε=−Δ\mathcal{L}_\varepsilon=-\Delta, whose domain is

D(Lε)={u∈H1(Ωε)|Δu∈L2(Ωε), u=0 on ΓDε, ∂nu=0 on ΓNε}.(4)\mathcal{D}(\mathcal{L}_\varepsilon)=\left\{u\in H^1(\Omega_\varepsilon)\mathrel{}\middle|\mathrel{} \Delta u\in L^2(\Omega_\varepsilon),\ u=0\ \text{on }\Gamma_\mathrm{D}^\varepsilon,\ \partial_nu=0\ \text{on }\Gamma_\mathrm{N}^\varepsilon\right\}. \tag*{(4)}

where the normal derivative ∂nu\partial_nu is defined by duality on H−1/2(∂Ωε)H^{-1/2}(\partial\Omega_\varepsilon) (see Section B for an explicit definition). The operator Lε\mathcal{L}_\varepsilon can be defined as the self-adjoint operator associated with the closed form

References

  1. [22]D. Frenkel. Entropy-driven phase transitions. Physica A: Statistical Mechanics and its Applications, 263(1-4):26–38, 1999. DOI: 10.1016/S0378-4371(98)00501-9. Proceedings of the 20th IUPAP International Conference on Statistical Physics.DOI
  2. [23]D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second order. Classics in Mathematics. Springer-Verlag, Berlin, 2001. DOI: 10.1007/978-3-642-61798-0. Reprint of the 1998 edition.DOI
  3. [24]D. Gomez and A. F. Cheviakov. Asymptotic analysis of narrow escape problems in nonspherical three-dimensional domains. Phys. Rev. E, 91(1):012137, 2015. DOI: 10.1103/PhysRevE.91.012137.
  4. [25]P. Hänggi, P. Talkner, and M. Borkovec. Reaction-rate theory: fifty years after Kramers. Reviews of Modern Physics, 62(2):251–341, Apr. 1990. DOI: 10.1103/RevModPhys.62.251.
  5. [26]D. Hoff. Concentration properties and pointwise bounds for the Green’s functions for the Laplace operator. J. Differential Equations, 427:601–640, 2025. DOI: 10.1016/j.jde.2025.01.083.
  6. [27]D. Holcman and Z. Schuss. Escape through a small opening: receptor trafficking in a synaptic membrane. Journal of Statistical Physics, 117(5-6):975–1014, 2004. DOI: 10.1007/s10955-004-5712-8.
  7. [28]D. Holcman and Z. Schuss. The narrow escape problem. SIAM Review, 56(2):213–257, 2014. DOI: 10.1137/120889395.DOI
  8. [29]M. S. Howe. On the Helmholtz resonator. Journal of Sound and Vibration, 45(3):427–440, 1976. DOI: 10.1016/0022-460X(76)90397-7.DOI
  9. [30]W.-T. Hsu. Asymptotic analysis on narrow tubes: narrow escape problems and diffusion processes, 2025. arXiv: 2508.15060 [math.AP].
  10. [31]W.-T. Hsu. Metastability of diffusion processes in narrow tubes, 2025. arXiv: 2511.21548 [math.PR].
  11. [32]U. Ingard. On the theory and design of acoustic resonators. The Journal of the Acoustical Society of America, 25(6):1037–1061, Nov. 1953. DOI: 10.1121/1.1907235.
  12. [33]S. Jimbo. The singularly perturbed domain and the characterization for the eigenfunctions with Neumann boundary condition. Journal of Differential Equations, 77(2):322–350, 1989. DOI: 10.1016/0022-0396(89)90147-2.
  13. [34]S. Jimbo and S. Kosugi. Spectra of domains with partial degeneration. J. Math. Sci. Univ. Tokyo, 16(3):269–414, 2009.
  14. [35]L. Koons and E. Lauga. Slender-ribbon theory. Physics of Fluids, 28(1):013101, Jan. 2016. DOI: 10.1063/1.4938566.
  15. [36]H. A. Kramers. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica, 7(4):284–304, 1940. DOI: 10.1016/S0031-8914(40)90098-2.
  16. [37]T. Lelièvre, D. Le Peutrec, and B. Nectoux. Eyring-Kramers exit rates for the overdamped Langevin dynamics: the case with saddle points on the boundary. J. Éc. polytech. Math., 12:881–982, 2025. DOI: 10.5802/jep.303. arXiv: 2207.09284 [math.PR].
  17. [38]T. Lelièvre, M. Rachid, and G. Stoltz. A spectral approach to the narrow escape problem in the disk, 2024. arXiv: 2401.06903 [math.AP].
  18. [39]X. Li. Matched asymptotic analysis to solve the narrow escape problem in a domain with a long neck. J. Phys. A, 47(50):505202, 2014. DOI: 10.1088/1751-8113/47/50/505202.
  19. [40]A. E. Lindsay, A. J. Bernoff, T. Goodwill, and J. G. Hoskins. The three-dimensional Neumann Green’s function for general surfaces: singular asymptotics and boundary integral methods, 2026. arXiv: 2606.12714 [math.NA].
  20. [41]M. Mascagni and C.-O. Wang. ε-shell error analysis for “Walk on spheres” algorithms. Math. Comput. Simulation, 63(2):93–104, 2003. DOI: 10.1016/S0378-4754(03)00038-7.
  21. [42]B. Miller, R. Sawhney, K. Crane, and I. Gkioulekas. Walkin’ Robin: walk on stars with Robin boundary conditions. ACM Trans. Graph., 43(4):Art. 41, 2024. DOI: 10.1145/3658153.
  22. [43]P. M. Morse and P. J. Rubenstein. The diffraction of waves by ribbons and by slits. Physical Review, 54(11):895–898, Dec. 1938. DOI: 10.1103/physrev.54.895.
  23. [44]M. E. Muller. Some continuous Monte Carlo methods for the Dirichlet problem. Ann. Math. Statist., 27(3):569–589, 1956. DOI: 10.1214/aoms/1177728169.
  24. [45]M. Nursultanov, J. C. Tzou, and L. Tzou. On the mean first arrival time of Brownian particles on Riemannian manifolds. J. Math. Pures Appl. (9), 150:202–240, 2021. DOI: 10.1016/j.matpur.2021.04.006.DOI
  25. [46]W. Prager and J. L. Synge. Approximations in elasticity based on the concept of function space. Quart. Appl. Math., 5(3):241–269, 1947. DOI: 10.1090/qam/25902.DOI
  26. [47]V. Richardson, Y. H. Ling, and S. D. Lawley. Yet another look at narrow escape through a tube. SIAM Journal on Applied Mathematics, 86(4):1632–1653, 2026. DOI: 10.1137/25M1824197.DOI
  27. [48]D. Ruiz. On the uniformity of the constant in the Poincaré inequality. Adv. Nonlinear Stud., 12(4):889–903, 2012. DOI: 10.1515/ans-2012-0413.DOI
  28. [49]J.-F. Rupprecht, O. Bénichou, D. S. Grebenkov, and R. Voituriez. Exit time distribution in spherically symmetric two-dimensional domains. Journal of Statistical Physics, 158(1):192–230, 2015. DOI: 10.1007/s10955-014-1116-6.DOI
  29. [50]K. K. Sabelfeld. Monte Carlo Methods in Boundary Value Problems. Springer Series in Computational Physics. Springer-Verlag, Berlin, 1991. DOI: 10.1007/978-3-642-75977-2.
  30. [51]S. Salsa. Partial Differential Equations in Action: From Modelling to Theory. Universitext. Springer-Verlag Italia, Milan, 2008. DOI: 10.1007/978-88-470-0752-9.DOI
  31. [52]G. Savaré. Regularity and perturbation results for mixed second order elliptic problems. Comm. Partial Differential Equations, 22(5-6):869–899, 1997. DOI: 10.1080/03605309708821287.DOI
  32. [53]R. Sawhney and K. Crane. Monte Carlo geometry processing: a grid-free approach to PDE-based methods on volumetric domains. ACM Trans. Graph., 39(4):Art. 123, 2020. DOI: 10.1145/3386569.3392374.DOI
  33. [54]R. Sawhney, B. Miller, I. Gkioulekas, and K. Crane. Walk on stars: a grid-free Monte Carlo method for PDEs with Neumann boundary conditions. ACM Trans. Graph., 42(4):Art. 80, 2023. DOI: 10.1145/3592398.DOI
  34. [55]R. Sawhney, D. Seyb, W. Jarosz, and K. Crane. Grid-free Monte Carlo for PDEs with spatially varying coefficients. ACM Trans. Graph., 41(4):Art. 53, 2022. DOI: 10.1145/3528223.3530134.DOI
  35. [56]A. Singer, Z. Schuss, and D. Holcman. Narrow escape and leakage of Brownian particles. Phys. Rev. E, 78(5):051111, Nov. 2008. DOI: 10.1103/PhysRevE.78.051111.DOI
  36. [57]A. Singer, Z. Schuss, and D. Holcman. Narrow escape, Part III: non-smooth domains and Riemann surfaces. J. Stat. Phys., 122(3):491–509, 2006. DOI: 10.1007/s10955-005-8028-4.DOI
  37. [58]A. Singer, Z. Schuss, D. Holcman, and R. S. Eisenberg. Narrow escape, Part I. Journal of Statistical Physics, 122(3):437–463, 2006. DOI: 10.1007/s10955-005-8026-6.DOI
  38. [59]J. Stone and N. Matsakis. Rayon: a data parallelism library for Rust. https://github.com/rayon-rs/rayon, 2024.
  39. [60]H. Tanaka. Stochastic differential equations with reflecting boundary condition in convex regions. Hiroshima Math. J., 9(1):163–177, 1979. DOI: 10.32917/hmj/1206135203.DOI
  40. [61]The Rust Team. The Rust programming language. https://rust-lang.org, 2020.
  41. [62]A. F. Voter. Introduction to the kinetic Monte Carlo method. In K. E. Sickafus, E. A. Kotomin, and B. P. Uberuaga, editors, Radiation Effects in Solids. Volume 235, NATO Science Series, pages 1–23. Springer, Dordrecht, 2007. DOI: 10.1007/978-1-4020-5295-8_1.DOI
  42. [63]Y. Wang and J. Li. From reflecting Brownian motion to reflected stochastic differential equations: a systematic survey and complementary study. SSRN Electronic Journal, Oct. 2020. DOI: 10.2139/ssrn.3688563.DOI
  43. [64]M. J. Ward and J. B. Keller. Strong localized perturbations of eigenvalue problems. SIAM J. Appl. Math., 53(3):770–798, 1993. DOI: 10.1137/0153038.DOI
  44. [65]R. Zwanzig. Diffusion past an entropy barrier. The Journal of Physical Chemistry, 96(10):3926–3930, May 1992. DOI: 10.1021/j100189a004.DOI

Paper details

vlnc:W_pr7tdytk
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 4

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Claude Opus 5.5 · View statement →

Looking for related papers…
Contents