Alan E. Lindsay,Andrew J. Bernoff
Abstract:Cellular scale decision making is modulated by the dynamics of signalling molecules and their diffusive trajectories from a source to small absorbing sites on the cellular surface. Diffusive capture problems are computationally challenging due to the complex geometry and the applied boundary conditions together with intrinsically long transients that occur before a particle is captured. This paper reports on a particle-based Kinetic Monte Carlo (KMC) method that provides rapid accurate simulation of arrival statistics for (i) a half-space bounded by a surface with a finite collection of absorbing traps and (ii) the domain exterior to a convex cell again with absorbing traps. We validate our method by replicating classical results and in addition, newly developed boundary homogenization theories and matched asymptotic expansions on capture rates. In the case of non-spherical domains, we describe a new shielding effect in which geometry can play a role in sharpening cellular estimates on the directionality of diffusive sources.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how particles diffuse and are trapped onto small absorption sites on the cell surface in three - dimensional space under complex geometric structures. Specifically, the paper focuses on calculating the time distribution of particles arriving at these absorption sites, especially the diffusion - trapping problems outside the half - space, the plane, and convex polyhedra. These problems are very common in biological phenomena and engineering systems, such as cell - signaling, cargo transportation, etc.
### Main problems
1. **Calculating the arrival - time distribution**: For particles released from position \(x_0\), the key problem is to calculate the dynamic flux on each absorption site and its dependence on the number and spatial configuration of absorption sites.
2. **Complex geometries and mixed boundary conditions**: Due to the complexity of geometric structures and mixed boundary conditions (partially reflective, partially absorptive), such problems are challenging in numerical simulations.
3. **Long - transient processes**: Before the particles are trapped, long - transient processes usually occur, which makes the simulation more difficult.
### Mathematical model
The diffusion problem in the paper can be described by the following partial differential equation:
\[
\frac{\partial p}{\partial t}=D \Delta p, \quad x \in \Omega, \quad t > 0
\]
Initial condition:
\[
p(x, 0)=\delta(x - x_0), \quad x \in \Omega
\]
Boundary conditions:
\[
p = 0\quad \text{on}\quad x \in \Gamma_a
\]
\[
D \nabla p\cdot\hat{n}=0\quad \text{on}\quad x \in \Gamma_r
\]
where \(D\) is the diffusion coefficient of the particles, \(\Omega\) is a three - dimensional region, and its boundary \(\partial\Omega\) is divided into an absorptive part \(\Gamma_a\) and a reflective part \(\Gamma_r\).
### Solutions
The paper proposes a particle - based Kinetic Monte Carlo (KMC) method to quickly and accurately simulate the arrival statistics of particles by sampling Brownian motion. The specific steps include:
1. **Projection operators**: Combine the heat - equation projection operators of three exact solutions to handle projections from a point to a plane, from a point to a sphere, and from the center of a hemisphere to a sphere respectively.
2. **Mixed boundary - value problems**: The KMC method can bypass the obstacles of traditional methods (such as the finite - element method, the finite - difference method, and the boundary - integral method) when dealing with surface - flux singularities.
3. **Convergence**: Although the accuracy of the KMC method converges slowly, the precision can be improved by parallel - computing a large number of particles.
### Applications and verification
- **Verification of classical results**: The paper verifies the effectiveness of the method by replicating classical results and verifying the new boundary - homogenization theory and matched asymptotic expansions.
- **New effects in non - spherical domains**: In non - spherical domains, the geometric shape can affect the directional estimation of the diffusion source, and the paper describes this new shielding effect.
In conclusion, this paper aims to develop an efficient numerical method for quickly and accurately calculating the diffusion - trapping problems of particles in three - dimensional convex geometric structures, especially in cell - scale decision - making processes.